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Abstract 


This  work  examines  the  angular  correlation  of  speckle  patterns  produced  by  the 
coherent  illumination  of  random  rough  surfaces.  Speckle  patterns  are  produced 
when  coherent  light  is  scattered  by  a  rough  surface  or  other  random  medium.  As 
the  angle  of  illumination  changes,  the  speckle  pattern  changes  at  a  rate  determined 
by  geometry  and  by  the  details  of  the  scattering  process.  It  is  this  angular  change 
that  is  studied  in  this  thesis  by  means  of  the  correlation  coefficient. 

Experimental  studies  are  made  of  the  angular  change  of  speckle  patterns  from  low 
and  high  sloped  Gaussian  random  rough  surfaces.  They  are  illuminated  with  a 
linearly  polarised  laser  lx  am,  and  the  correlation  coefficient  of  speckle  patterns 
observed  with  a  CCD  camera  then  being  determined  by  computational  analysis. 

The  Kirchhoff  approximation  is  applied  to  the  Ilelmholtz-Kirchhoff  integral  to  obtain 
a  two-dimensional  analytical  solution  for  the  correlation  produced  by  a  low  slope, 
single  scattering.  Gaussian  surface.  The  results  are  compared  with  the  experimental 
measurements. 

A  simple  statistical  examination  of  Gaussian  random  rough  surfaces  is  made;  this 
shows  that  the  KirchhofT  approximation  is  still  valid  even  for  surfaces  which  exhibit 
multiple  scattering.  The  Kirchhoff  approximation  is  used  to  derive  an  expression  for 
the  scattered  light  from  a  high  sloped,  multiply  scattering,  surface.  This  required 
computational  implementation  —  an  advantage  being  that  surface  shading  can  be 
accounted  for.  A  very  fast  computational  method  for  determining  shaded  areas  on 
a  surface  is  presented.  Polarisation  changes  in  the  scattered  light  are  catered  for. 

A  conceptual  model  of  the  processes  involved  in  multiple  scattering  is  developed, 
successfully  explaining  various  features  of  observed  results. 
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Chapter  1 


Introduction 

1.1  Introduction  and  Literature  Review 

Interest  in  light  scattering  from  random  media  has  undergone  a  noticeable  revival 
within  the  last  decade.  Although  some  of  this  renewed  interest  may  be  attributable 
to  the  more  recent  computer  technology,  enabling  solutions  to  various  scattering 
problems  to  be  determined  numerically,  the  subject  has  generated  interest  in  the 
medical  field.  In  this  area  the  volume  scattering  and  absorption  properties  of  a  pa¬ 
tient  may  potentially  give  information  about  internal  processes  without  the  necessity 
of  large  and  expensive  equipment  or  the  use  of  harmful  radiation. 

The  driving  forces  behind  the  study  of  scattering  from  random  rough  surfaces  arc  its 
potential  as  a  non-destructive  or  remote  probe  for  measuring  the  statistical  proper¬ 
ties  of  rough  surfaces,  and  the  inescapable  effect  that  the  inherent  randomness  in  all 
surfaces  has  on  everyday  scattering  applications.  Radar  falls  into  this  last  category, 
and  may  be  regarded  as  a  scaled  up  light  scattering  problem;  similar  results  may  be 
obtained  using  light  in  a  laboratory  as  would  be  gathered  from  a  more  expensive, 
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CHAPTER  1.  INTRODUCTION 


less  controllable,  radar  set  up. 

1.1.1  Theoretical  Models 

This  section  will  present  some  of  the  theoretical  models  which  exist  for  determining 
the  problem  of  light  scattering  from  random  rough  surfaces.  It  aims  to  give  an 
introduction  to  some  of  the  methods  not  used  within  this  work. 

Approximate  theories  exist  for  very  small  height  deviations  across  the  surface  (such 
as  perturbation  theory)  and  for  surfaces  which  have  sufficiently  low  gradients  such 
that  they  may  be  considered  as  single  scattering  (Beckmann  theory  or  single  scat¬ 
ter/first  order  Kirchhoff  approximation).  Other  methods  attempt  to  predict  the 
effects  of  multiple  scattering  and  include  the  effects  of  surface  shadowing.  Compu¬ 
tational  methods  of  solving  the  integral  equations  also  exist  which  avoid  the  necessity 
of  approximations.  Most  of  the  theories  discussed  below  may  be  found  in  more  depth 
in  the  book  by  Ogilvym. 

Perturbation  Theory 

Perturbation  theory  can  be  used  for  slightly  rough  surfaces  which  deviate  from  a 
flat  plane  by  only  a  small  amount. 

The  field  on  a  surface  ==  h{r,y)  may  be  written  as 

E  =  E,  +  Es  (1.1) 

where  the  subscripts  i  and  s  represent-  the  incident  and  scattered  fields  respectively. 
E  is  a  vector  quantity  of  a  particular  polarisation. 
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The  perturbation  theory  assumes  that  the  scattered  field  may  be  written  as  a  series: 

Ea  =  Ea  o  +  Ea  i  +  Ea  2  +  •  •  •  ( 1  >2) 

where  the  numerals  represent  the  different  orders. 

Assuming  Ei  and  Ea  to  be  functions  of  A,  taking  a  Taylor  expansion  in  terms  of  A, 
terms  of  the  same  order  in  I)  may  be  equated  after  applying  the  Dirichlet  boundary 
condition,  namely  that  the  field  at  the  surface  is  zero 

£|«*  =  0  (1.3) 


at  the  mean  plane  (c  =  0).  Applying  the  Dirichlet  boundary  condition  when  the 
surface  heights  are  much  smaller  than  a  wavelength,  the  zeroth  order  term  in  the 
perturbation  series  (the  field  scattered  from  a  flat  plane)  is  found: 

£so|i=0  —  —  Ei\z-u  (1-4) 


and  the  first,  order  perturbation  term  (representing  the  effects  of  the  deviations  in 
the  surface  from  a  smooth  plane)  is 
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(1.5) 


The  llelmholtz-Kirclihoff  integral  (discussed  in  chapter  two): 


(1.6) 


may  now  be  used,  with  equation  (1.5)  as  a  boundary  condition,  to  give  the  scattered 
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where  S  is  the  mean  surface  plane  and  xj)  is  the  half  space  Green’s  function;  all  other 
notation  is  as  previously  defined. 


A  similar  method  may  be  used  to  obtain  the  result  for  the  Neumann  boundary 
condition, 


dE 


On 


-Uo  =  0 


(1.8) 


This  may  be  used  when  the  surface  gradients  are  much  less  than  one,  and  gives 
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It  is  possible,  again  with  a  similar  method,  to  obtain  the  second  order  perturbation 
of  the  scattered  field.  The  original  Taylor  expansion  must  be  taken  out  to  second 
order  (terms  in  h2)  to  obtain  the  new  boundary  condition  to  apply  to  the  Ilelmholtz- 
Kivchhoir  integral. 

The  accuracy  of  perturbation  theory  has  been  studied  by  many  authors  (cf.  section 
3.1.3  of  reference  [1]  for  a  detailed  review),  the  generally  accepted  range  of  validity 
being  ka  <C  1  where  a  is  the  rms  height,  of  the  surface,  t  the  correlation  length  and 
k  the  wavenumber  One  recent  work  compares  perturbative  calculations  with 
“exact  numerical  results'’  (as  used  in  references  [2],  [3]  and  [■!])  and  those  obtained 
by  using  the  Kirchhoff  approxiinationl’l.  It  finds  that  for  Ict  >  1  (where  7  is  the 
correlation  length  of  the  surface),  the  fourth  order  perturbative  expansion  is  better 
than  the  second  order  one.  the  latter  being  accurate  when  the  total  incoherently 
scattered  energy  is  less  than  about  9%  of  the  incident  energy.  For  the  fourth  order 
expansion,  the  theory  was  found  to  be  accurate  when  the  total  incoherently  scattered 
energy  is  less  than  around  3 9'A  of  the  incident-  energy.  If  kr  <  1,  it  was  found  that 
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the  perturbative  series  is  net  convergent,  and  the  second  order  expansion  gives  better 
results  than  the  fourth. 

Surfaces  considered  in  this  thesis  have  an  rms  roughness,  ko  1,  hence  the  pertur¬ 
bation  approximation  will  be  of  little  value. 

Extinction  Theorem 

The  extinction  theorem  derives  its  name  from  the  use  of  the  condition  that  the  field 
within  a  scattering  surface  must  cancel  the  incident  field.  For  finite  surfaces,  this 
condition  is  approximated  to  include  all  space  below  *he  surface  plane.  .  is  also 
known  as  the  null-field  method  and  the  extended  boundary  condition  method. 

Starting  with  the  usual  Hclmholtz-KirchhoiT  integral  (ecpiation  (1.6))  for  the  scat¬ 
tered  field,  and  writing  the  incident  field  such  that  the  field  below  the  surface  is  zero 
—  i.c.  the  incident  field  cancels  out  the  scattered  field  (£,  +  Ea  =  0),  then 


Equation  (1.10)  may  be  solved  for  E ,  which  may  then  be  substituted  into  equation 
(1.6)  to  give  the  scattered  field.  In  practice  it  is  usually  necessary  to  find  an  ap¬ 
proximate  solution  to  equation  (1.10),  by  using  an  expansion  similar  to  that  used 
in  perturbation  techniques.  Numerical  simulation  avoids  this  necessity  by  using  an 
iterative  technique  to  obtain  an  exact  resulted!'1]. 

Whilst  the  extinction  theorem  provides  an  exact  solution  in  principle,  the  scattered 
intensity  cannot  be  written  as  a  single  closed  form  expansion,  and  therefore  it  is  of 
little  value  in  providing  support  for  a  physical  picture  of  light  .scattering. 
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Kirchhoff  Theory 

Kirchhoff  theory  is  used  extensively  in  later  chapters  of  this  thesis  and  will  not  be 
dealt  with  in  this  section.  The  single  scattering  Beckmann  theory  is  discussed  in 
chapter  four,  and  the  double  scattering  technique  used  by  Bruce^7^8!  and  others^9^10^11^ 
is  explained  in  chapter  five. 

One  technique  increasingly  used  starts  with  the  extinction  theorem  and,  by  iterating 
the  surface  field  derived  n  times,  gives  increasing  terms  in  a  series  of  which  the  n  ~  0 
term  is  the  single  scatter  KirchhofT  solution,  larger  n  giving  terms  representing 
double  and  higher  order  scattering  termsb2h13'.  Computational  techniques  must  be 
used  to  solve  a  large  matrix  equation. 

1.1.2  Mean  Scattered  Intensity 

Research  into  light  scattering  from  random  rough  surfaces  has  remained  fairly  active 
since  the  leporting  by  Mendez  and  O’l)onnclli1,h15l  of  ‘enhanced  backseatter’.  This 
work  was  stimulated  by  the  observation  of  a  similar  effect,  in  volume  scattering^. 
The  effect  itself  has  been  known  of  for  years  in  the  astronomical  field,  referred  to  as 
the  ‘opposition  effect’  [Montgomery!1 '1  and  refs,  therein].  The  name  came  from  the 
fact  that  it  was  noticed  in  scattering  from  visible  bodies  in  the  solar  system  when 
they  and  the  sun  were  in  opposition,  i.e.,  on  opposite  sides  of  the  earth.  The  signif¬ 
icance  of  Mendez  and  O'Donnell’s  work  was  that  well  characterised  random  rough 
surfaces  were  used,  following  the  method  of  Gray!18!  (<,ee  chapter  three),  enabling 
confident  comparisons  with  computational  results  and  the  existing  approximate  an- 
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It  was  found  that  Beckmann  theory^  agreed  well  with  experimental  results  from  a 
low  slope  surface  for  small  angles  of  incidence.  A  study  of  two  steeper  slope  surfaces 
was  undertaken,  showing  the  enhanced  backscatter  phenomenon,  and  a  possible 
explanation  for  the  backscattering  process  was  outlined  —  the  enhancement  being 
due  to  constructive  interference  of  equal  and  opposite  multiply  scattering  paths. 

Experiments  have  been  conducted  comparing  low  sloped  surfaces  to  theory^  and 
various  others  examine  the  scattering  characteristics  of  known  steep  slope  surfaces 
which  display  the  enhanced  backscattering  phenomenon.  Included  in  these  is  a  com¬ 
parison  of  dielectric  and  perfectly  conducting  (gold)  surfaces  with  ident  ical  profiles^21! 
—  epoxy  resin  copies  of  the  originally  fabricated  gold  coated  surface  were  made  via 
a  silicone  gel  mould.  A  gold  coated  reproduction  was  compared  experimentally  with 
the  master  and  the  results  confirm  that  the  technique  faithfully  copies  the  original. 
Significant  enhanced  backscattering  (noticed  in  the  gold  surface)  was  not  observed 
experimentally  in  the  dielectric  case.  Results  of  a  numerical  calculation  were  com¬ 
pared  with  the  experimental  dielectric  scattering  results  showing  broad  agreement. 

Results  from  various  one-  and  two-dimensional  surfaces  have  also  been  presented  for 
different  wavelengths  of  illumination^  and  a  comparison  of  results  from  statistically 
similar  one-  and  two-dimensional  surfaces  has  been  made^  showing  results  which 
are  similar  in  as  far  as  they  may  be  compared. 

In  reference  [23]  the  Stokes’  parameters  of  the  average  scattered  light  are  utilised 
to  present  experimental  data  of  the  polarised  and  unpolarised  components.  The 
backscatter  peak  was  found  only  in  the  unpolarised  component,  and  is  therefore 
believed  to  be  a  product-  of  multiple  scattering.  Numerical  calculation  results  were 
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compared  favourably  with  experimental  results  from  one-dimensional  surfaces. 

More  extensive  studies  have  since  been  made  of  the  polarisations  of  scattered  fields 
and  the  Stokes’  parameters.  Four  unique  elements  were  found  in  the  Mueller  matrix 
of  a  one-dimensional  surface,  all  of  which  were  shown  to  be  significant  in  scattering 
from  such  a  surface  which  produces  enhanced  backscattering  t'24l,  The  complete 
Mueller  matrix  for  a  one  dimensional  perfectly  conducting  random  rough  surface 
has  been  numerically  and  experimentally  found^’3)  for  all  observed  angles  and  for  a 
number  of  different  angles  of  incidence  —  the  results  were  generally  similar  to  each 
other. 

A  complete  Mueller  matrix  has  been  numerically  and  experimentally  calculated  for 
a  gold  and  dielectric  one  dimensional  surface^.  Agreement  between  numerical  and 
experimental  results  was  good  for  the  dielectric  case,  but  differed  slightly  in  the  case 
of  the  gold  surface.  The  difference  u'as  shown  to  be  in  the  matrix  term  (mapping 
linearly  polarised  incident-  light  to  circularly  polarised  scattered  light),  and  believed 
to  be  due  to  a  difference  in  the  reflection  coefficients  of  the  surface  and  those  used 
in  the  calculation.  The  Stokes’  parameters  are  examined  and  utilised  as  in  reference 
[23]  to  show'  the  polarised  and  unpolarised  scattering  components. 

A  presentation  was  given  of  numerical  calculation  results  of  the  angular  cont  ribution 
in  the  intensity  of  the  diffuse  component  of  p-polarised  light  scattered  from  a  random 
grating  ruled  on  penetrable  and  impenetrable  dielectric  medial.  It  is  found  that 
a  metallic  (silver)  surface  gives  similar  results  to  a  perfectly  conducting  one.  This 
result  has  a  large  impact  on  the  numerical  study  of  light  scattering  as  it  is  much 
simpler  to  implement  a  perfectly  conducting  surface  than  a  metallic  one. 
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Multiple  scattering  and  shadowing  extensions  to  the  simple  KirchhofF  theory  have 
been  implemented^'^  (see  chapter  five)  —  numerical  resuits  for  single  and  double 
scatter  components  of  the  Kichhoff  theory  are  presented,  for  gold  and  dielectric 
one  dimensional  surfaces,  and  compared  favourably  with  experiment.  Experimental 
results  have  also  been  compared  to  the  double  scatter  Kirchhoff  theory  for  a  gold 
and  a  dielect  ric  surface  which  have  identical  height  profiles  t25h  Various  wavelengths 
were  used.  Generally,  the  double  scatter  Kirchhoff  theory  was  found  to  be  in  good 
agreement  with  experimental  results  for  angles  of  incidence  of  less  than  sixty  degrees. 


The  ‘exact’  method  of  the  iterative  solution  of  the  extinction  theorem^  has  been 
used  in  a  Monte-Carlo  numerical  calculation^  in  conjunction  with  the  single  seal  ter 
Kirchhoff  theory  for  a  range1  of  surface  statistics,  to  define  the  limits  of  the  single 
scatter  Kirchhoff  approximation.  The  results  presented  do  not  define  the  limits  of 
where  the  Kirchhoff  approximation  can  be  applied,  but  where  t  he  single  scattering 
theory  based  upon  the  approximation  can  be  expected  to  succeed. 


Other  work  outside  the  scope  of  correlation  measurements,  but  not  so  closely  related 
to  the  work  to  be  presented  in  the  following  chapters,  includes  joint  experimental 
and  theoretical  studies  of  the  backscattering  from  and  transmission  through  dielec¬ 
tric  and  thin  film  random  rough  surfaces^Gh27h  Speckle  techniques  have  also  been 
developed^  —  speckle  contrast  measurement s^2;)^30^  are  amongst  the  many  different 
met  hods  discussed. 
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1.1.3  Correlation  Studies 

Correlation  studies  were  first  carried  out  in  the  early  to  mid-seventies.  However 
at  this  time  the  interest  in  multiply  scattering  surfaces  had  not  developed  and  all 
works  dealt  with  low  slope  surfaces. 

Pedersen^  discussed  speckle  correlation,  giving  a  theoretical  explanation  of  var¬ 
ious  results  reported  prior  to  his  study.  Most  of  his  interest,  however,  lay  with 
polychromatic  speckle^32!  and  his  work  is  not  relevant  to  this  thesis. 

The  work  by  Leger  and  Perrin  to  determine  surface  roughness  via  a  correlation 
technique  is  directly  relevant  to  this  work.  By  exposing  a  photographic  plate  t.o 
speckle  produced  by  a  rough  surface  at  two  angles  of  incidence,  Young's  fringes 
could  be  observed  in  the  Foiuier  plane  when  the  processed  plate  was  placed  in  a. 
converging  beam^33!.  The  contrast  of  the  fringes  was  mathematically  related  to  the 
surface  roughness,  starting  with  the  general  Beckmann  solution  for  the  scattered 
field,  but  then  working  in  the  l'ourier  domain. 

In  a  later  paper3']  a  more  general  derivation  of  the  fringe  contrast  from  angular 
correlations  was  reported.  This  was  again  based  on  the  single  scattering  Beckmann 
theory,  but  took  a  real  space  approach  to  the  problem.  It  is  the  basis  of  the  derivation 
to  follow'  in  chapter  four. 

More  recently.  Ohlidafi3']  compared  the  equivalent  of  Leger  and  Perrin’s  method 
(which  uses  the  Fraunhofer  approximation)  to  results  obtained  using  the  Fresnel 
approximation.  No  comparison  with  experimental  measurements  was  made,  but  an 
analysis  of  when  each  method  should  apply  is  presented. 
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The  most  recent  studies  in  this  are?  have  been  performed  by  Michel  et.  al.  and 
Knotts  et  al.  Michel  et  alJ36l  present  numerical  results  for  the  angular  correlations 
of  scattered  amplitudes  from  a  one  dimensional  surface  and  analyse  the  results  into 
those  attributable  to  single  and  double  scattering  Knotts  et  aU37!  offer  experimental 
results  for  angular  correlations  obtained  from  a  one-dimensional  surface  and  compare 
with  numerical  results.  A  discussion  as  to  the  cause  of  various  features  observable 
in  the  results  is  given. 

The  Memory  Effect 

The  memory  effect  was  named  by  Feng  et  alJ38^  for  situations  of  volume  scattering 
(scattering  from  a  media  in  which  the  light  may  travel,  e.g.,  paint  or  milk).  They 
show  various  elements  contributing  to  the  correlation  function,  and  show  that  certain 
of  these  are  large  when  the  “momentum  transfer  of  the  incident  beams  equals  the 
transmitted  ones",  i.o..  the  change  in  the  wave  vectors  upon  scattering  is  the  same 
in  both  cases  being  correlated,  Aq!  —  Aq^.  Aq  is  defined  as  q„  —  q;,,  tin*  difference 
in  the  (transverse)  wave  vectors  of  the  incident  and  scattered  light. 

Applying  this  to  in  plane  surface  scattering  correlations,  the  memory  effect  condition 
holds  (see  chapter  four)  when 


sin  0,  —  sin  03  =  sin  <p,  —  sin  <ps  (1.11) 

where  the  subscripts  i  and  a  represent  the  incident  and  scattered  angles  respectively, 
0  being  used  to  obtain  one  speckle  pattern  and  <p  being  used  to  obtain  the  second. 


This  condition  was  used  by  Leger  and  Perrin^  twelve  years  earlier  in  their  work 


38 


CHAPTER  1.  INTRODUCTION 


described  above,  although  it  was  expressed  in  a  different  form: 

cos  0,60  —  cos  <j),b4>  (1.12) 

Equation  (1.11)  may  be  reconciled  to  equation  (1.12)  by  noting  that  0,  and  <fis 
have  been  replaced  by  0,  +  SO  and  +  S<f>  respectively  and  by  applying  a  simple 
trigonometric  transform.  As  they  were  only  dealing  with  small  changes  in  SO  and 
6<fr,  terms  of  order  SO 2  and  S<S2  can  be  removed. 

Pedersen^  also  noted  that  it  was  desirable  to  have  the  condition  satisfied,  and 
equation  (1.11)  has  been  used  in  the  work  reported  by  Michel  et  alJ3r'l. 

Nieto  Vesperinas  and  Sanchez-Ciil^  have  studied  the  memory  effect  in  relation  to 
surface  scattering,  but  they  also  notice  an  “enhanced  long  range  correlation"  ob¬ 
servable  when 

sin  0,  +  sin  0S  =  ±(sin  <?,  +  sin  </>,)  (1.13) 

The  origin  of  this  i.  believed  to  be  the  same  as  for  the  peak  observed  in  enhanced 
backseat t oring. 


1.2  Synopsis 

The  remaining  sections  of  this  chapter  consist  of  a  statement  of  the  original  work 
undertaken  and  a  description  of  the  scattering  geometry  to  be  used  throughout  the 
work. 

The  second  chapter  contains  general  theoretical  notions,  including  statistical  con¬ 
cepts  such  as  random  variables  and  processes.  It  also  contains  a  derivation  of 
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the  Helmholtz-Kirclihoff  integral  and  some  otlier  concepts  from  electromagnetism, 
namely  descriptions  of  reflection  coefficients  and  polarisation. 


Chapter  three  contains  all  of  the  experimental  methods  nsod;  these  include  manu¬ 
facturing  and  classifying  random  rough  surfaces  and  the  procedure  used  to  calculate 
the  correlation  coefficient  from  digitised  images  of  speckle  patterns. 


An  analytical  expression  for  the  correlations  expected  from  a  low-sloped,  single 
scattering  surface  is  derived  from  the  IIclmholtz-KirchhofT  integral,  in  chapter  four, 
by  applying  the  KirchhofT  approximation  to  it.  Experimental  results  for  such  a 
surface  are  compared  wit  it  the  theoretic,  al  result  and  discussed. 


The  fifth  chapter  extends  the  KirchhofT  result  to  include  the  double  scatter  term. 
The  equation  is  no  longer  analytically  soluble  and  is  implemented  computationally. 
Polarisation  changes  are  considered,  and  an  algorithm  for  rapidly  determining  Mines 
of  sight '  on  a  surface  is  presented.  It.  was  necessary  to  develop  this  algorithm  in  order 
to  calculate  point-to-point  shading  in  the  fastest  possible  time;  without  it  a  two- 
dimensional  computational  problem  of  this  nature  would  not  be  practical.  Results 
for  a  high-slopeu,  multiple  scattering  surface  are  compared  with  the  computational 
results. 


The  last  chapter  considers,  and  attempts  to  explain,  the  results  obtained.  A  possi¬ 
ble  mechanism  to  describe  the  observed  single  and  double  components  of  multiply 


scattering  surfaces  is  presented. 
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1.3  Statement  of  Originality 

There  is  not  a  great  deal  of  original  content  within  chapters  one,  two  and  three  — 
these  are  mainly  concerned  with  introducing  some  background  information  and  the 
concepts  and  techniques  available  for  use  within  this  thesis. 

Chapters  four  and  five  contain  new  experimental  results  for  the  intensity  correlation 
of  speckle  patterns  produced  by  a  method  described  in  chapter  three.  In  chapter 
four  an  existing  analytical  expression  is  extended  to  deal  with  two-dimensional  in¬ 
tensity  correlations  and  then  numerically  averaged  to  enable  a  comparison  with  the 
experimental  results.  In  chapter  five,  a  double  scatter  numerical  simulation  based 
on  the  Kirchhoff  approximation  is  extended  to  cope  with  two  dimensional  surfaces. 
A  new  fast  technique  for  calculating  the  shadows  cast,  on  the  surface  by  its  own 
fluctuations  is  presented. 

Chapter  shx  promotes  a  simple  conceptual  model  which  can  successfully  explain  the 
existence  o'  and  to  some  extent  predict  the  positions  of,  various  features  noted  in 
the  double  scatter  correlation  results  and  in  double  scatter  average  intensity  results. 


1.4  Correlation  Geometry 

Consider  a  plane  mirror  which  can  rotate  about  an  axis  which  lies  in  the  plane  of 
the  mirror.  If  an  observer  looks  in  the  direction  of  the  mirror,  he  will  see  an  image 
of  what  is  around  the  specular  direction  of  observation.  Now  if  the  mirror  is  rotated 
about  the  axis,  and  the  observer  moves  around  the  same  axis  until  he  can  sec  the 
same  scene  as  before,  t lieu  the  he  will  have  rotated  twice  as  much  as  the  mirror. 
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Figure  1.1:  The  geometry  to  be  used  to  obtain  the  experimental  and  theoretical 
results.  Observation  is  always  around  the  specular  direction. 

Extending  this  arrangement  to  the  correlation  of  speckle  patterns  formed  by  scatter¬ 
ing  from  a  random  rough  surface:  if  the  scattering  surface  is  considered  as  random 
fluctuations  on  top  of  a  plane  mirror,  the  phase  variations  introduced  to  the  scat¬ 
tered  light  should  be  most  identical  at  different  incident  angles  if  always  observing 
in  the  specular  direction  and  illuminating  on  the  rotation  axis  (such  that  the  saint* 
random  lluctuatu  us  are  illuminated).  If  these  conditions  are  met,  only  the  change 
in  the  phase  variations  due  to  the  different  angles  at  which  light  is  incident  on  the 
surface  should  be  observed. 
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The  geometrical  arrangement  to  be  used  for  investigating  the  angular  correlation  of 
speckle  patterns  will  be  conducted  by  always  viewing  speckle  pai  erns  in  the  specular 
direction  to  the  incident  illumination  (see  figure  1.1).  The  argun  -nt  presented  above 
will  be  examined  mathematically,  for  the  case  of  a  Gaussian  random  rough  surface, 
later  in  chapter  four. 


Chapter  2 


Mathematical  Concepts 

2.1  Statistics 

la  this  section,  various  statistical  notions  and  formulae,  which  will  be  required  later, 
will  be  presented.  Unless  otherwise  stated,  all  variables  will  be  real  and  continuous. 
The  bulk  of  the  material  below  can  be  found  in  textbooks  (e.g.,  Papoulisf10!  or 
Goodman^"!). 

2.1.1  Probability 

Probability  may  be  defined  in  three  ways,  but  generally  all  of  the  definitions  are 
required  to  form  a  complete  picture,  each  definition  considered  in  isolation  lacking 
something  which  would  enable  it  to  stand  on  its  own.  In  fact,  the  most  mathe¬ 
matically  rigorous  definition  of  probability  (the  axiomatic  approach)  is  incapable  of 
assigning  a  numerical  probability  to  any  event  other  than  the  cert  ain  one,  so  it,  and 
the  reams  ol  theory  based  upon  it,  would  be  useless  in  the  real  world  without,  the 
other,  ‘dirtier',  relative  frequency  and  classical  (a  priori )  approaches. 
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The  definitions  of  probability  are: 

•  the  axiomatic  definition,  which  assigns  a  number  to  an  event,  the  number 
restricted  to  be  positive,  the  probability  of  the  certain  event  being  one,  and 
the  probability  of  two  mutually  exclusive  events  being  defined  as  the  sum  of 
the  two  individual  probabilities; 

•  the  relative  frequency  definition,  which  defines  probability  of  an  event  as  the 
fraction  of  the  number  of  times  the  event  occurs  over  the  number  of  exper¬ 
iments  performed  as  the  number  of  experiments  performed  tends  to  infinity; 
and 

•  the  classical  (a  priori)  definition,  which  counts  the  possible  outcomes  and 
defines  the  probability  as  the  fraction  of  the  number  of  outcomes  in  which  t  he 
event  is  expected  to  occur  in  a  single  experiment  over  the  total  number  of 
possible  outcomes. 

The  most  useful  single  definition  of  probability  is  a  modified  form  of  the  axiomatic 
approach.  The  probability  is  a  number  associated  with  an  outcome  of  a  random 
experiment  (an  experiment  where  the  outcome  is  not  previously  determined),  this 
number  being  a  linear  representation  of  the  likelihood  that  the  associated  outcome 
will  be  observed,  or  might  be  expected  to  be  observed.  The  probability  is  1  for  an 
outcome  which  will  definitely  happen,  0  for  an  outcome  which  could  never  happen. 

This  approach  combines  the  individual  definitions  and  incorporates  everything  nec¬ 
essary  to  develop  and  use  probability  theory.  It  can  be  easily  shown  that  since  the 
probability  is  linear,  the  probability  of  occurrence  of  any  one  of  some  mutually  ex- 
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elusive  outcomes  is  the  sum  of  their  individual  probabilities,  thus  fulfilling  the  last 
axiomatic  condition.  Moreover,  the  relative  frequency  and  classical  approaches  have 
both  been  incorporated,  such  that  something  which  happen  or  might  be  expected 
to  happen,  half  of  the  time  will  have  a  probability  of  0.5;  hence,  it  is  possible  to  as¬ 
sign  a  number  to  an  outcome  which  was  not  previously  possible  under  the  axiomatic 
definition. 

This  definition  is  presented  as  an  aid  to  the  understanding  of  probability  theory. 


2.1.2  Random  Variables 

Definition 

The  definition  of  a  random  variable  is  based  on  the  underlying  random  experiment 
and  the  “events”  or  possible  (numerical)  outcomes  of  this  experiment.  A  random 
variable  is  all  possible  outcomes  of  a  random  experiment  associated  with  their  prob¬ 
abilities.  A  discrete  example  may  clarify  this:  if  the  experiment  is  throwing  two  six 
sided  dice  and  finding  their  sum,  then  the  possible  outcomes  are  integers  between 
2  and  12  each  having  differing  probabilities  of  occurring  —  the  value  2  has  a  prob¬ 
ability  of  occurrence  of  1/36  (the  dice  showing  1  +  1),  whereas  the  value  7  has  a 
probability  of  occurrence  of  1/6  (the  dice  showing  1  +  6,2  +  5,3  +  4,4  +3,5  +  2 
or  6+1).  The  values  of  this  discrete  random  variable  will  occur  with  the  same 
frequency  as  the  outcomes  of  the  experiment. 
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Probability  Distribution  and  Density  Functions 

A  probability  distribution  function,  Fx{x),  is  a  function  associated  with  a  random 
variable,  X,  and  which  gives  the  probability  that  an  outcome  will  happen  by  x; 
i.e.,  the  probability  that  an  outcome  will  be  smaller  than  X\  is  F\'(a:i).  Note  that 
Fx{-oo)  =  0  and  Fx{o o)  —  1. 


A  probability  density  function,  Pjy(.t),  also  associated  with  a  random  variable,  X, 
assigns  a  probability  to  every  possible  outcome  x\  i.e.,  the  probability  of  the  outcome 
.Tj  occurring  is  Px{?\)-  Strictly  speaking,  the  probability  density  function  gives  the 
probability  P\{xi)(lx  that  an  outcome  is  between  X\  and  .ti  +dx.  A  discrete  random 
variable  may  be  considered  as  having  a  probability  density  function  consisting  of  a 
series  of  delta  functions. 


The  probability'  density  function  is  derived  from  the  probability  distribution  func¬ 
tion,  where  possible,  by 


P.  v  (■•»•) 


<IFx(r) 

(lx 


(2.1) 


The  inverse  relation 

Fx(x)  =  [  Px(v)<hl 

J  — oo 


(2.2) 


also  exists. 


Expected  Value 

T  iro  expected  value  (or  mean  or  statistical  average)  of  the  random  variable  A  is 
defined  as  the  sum  of  all  the  possible  outcomes  multiplied  by  their  individual  prob- 
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abilities: 

/oo 

xPx{x)dx  (2.3) 

■OO 

where  P\  is  the  probability  density  function.  When  x  is  a  possible  value  of  the 
random  variable  A',  then  Px(x)  is  its  probability.  The  mean  is  a  first  order  moment 
of  a  random  variable,  as  will  be  seen  below,  but  has  been  dealt  with  first  as  it  is  a 
concept  with  which  most  people  are  familiar  and  serves  to  introduce  the  notation. 

Moments  of  a  Random  Variable 

The  moments  of  a  random  variable  are  defined  as 

m,  =  E{ Xk}  =  r  xkPx(x)dx  (2.4) 

J  —  CO 

Here,  mQ  —  1  and  m\  is  the  mean,  77.  Generally  of  more  interest  are  the  central 
moments,  i.e.,  the  moments  of  a  random  variable  after  it  has  been  shifted  to  the 
centre  of  its  range  so  that  the  expected  value  is  zero: 

!‘k  =  E{[. V  -  =  f  {x  -  i)')h Px(x)dx  (2.5) 

J -00 

Here,  fio  =  1,  fi\  =  0  and  is  the  variance,  a2,  a  being  the  standard  deviation. 
The  more  familiar  form  of  the  variance  of  a  (discrete)  list  of  A'  numbers 


can  be  obtained  from  a  relationship  between  the  moments  and  central  moments 
(found  by  a  binomial  expansion  of  the  ( x  —  7/)A  term  in  equation  (2.5)), 

ft2  --  m2  -  7/2 


(2.7) 
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and  by  realising  that  in  this  case  of  equally  weighted  values, 


1  N 


(2.8) 


=i 


Joint  Moments  of  Two  Random  Variables 


The  joint  moments  of  two  random  variables  X  and  Y  are  given  by 

r  oo  roo 

mkl  =  E{XkY  }  =  /  /  xkylPXY(x,y)dxdy  (2.9) 

J  —OzJ -oo 

The  expected  value  of  a  function  of  two  random  variables  is  given,  more  generally, 
by 

/c*"»  roo 

I  f(*,y)PxY(*,y)dxdy  (2.10) 

-e±J — <x> 

The  order  of  the  moments  is  defined  as  k  +  /  (cf.  equation  (2.9)),  so,  for  example, 
the  first  order  joint  moments  of  two  random  variables  A'  and  V',  are  7??J0  -  yx  and 
”*oi  =  Vv 

Px y{x,y)  is  the  joint  probability  density  function  of  X  and  Y .  Generally,  this  must 
be  determined  by  finding  the  joint  distribution  function  FXy[x,y)  (the  probability 
that  X  is  smaller  than  ,r  and  V  is  smaller  than  y)  and  differentiating  with  respect, 
to  both  x  and  y.  Special  cases  of  interest  are: 


1.  Two  independent  random  variables 


Pxy(x.y)  =  Px{x)PY[y) 


(2.11) 


and 
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2.  Two  Gaussianly  distributed  random  variables,  X  and  Y,  with  means  of  i)x  and 
r/v,  variances  of  a l  and  cr^  and  a  correlation  coefficient,  Cxy- 


Px v{x,y)  = 


2lT(Tx<Tyy]\  -  C\ 


(2.12) 


exp  - 


2(1  -C\Y) 


XY 

2 


(.r  -  T]xy  2Cxy{x  -  l]x){y  -  Vv)  ( y  -  7/y)2 

'  o 


O  x0  y 


It  can  be  seen  from  the  above  that  two  uncorrelated  Gaussian  variables  are  also 
independent,  i.e.,  if  Cxy  =  0  then  equation  (2.12)  takes  on  the  form  of  equation 
(2.11). 


The  joint  central  moments  of  two  random  variables  are  defined  by 

/'«  -•  E{(X  ~  Vr)k(Y  -  //„)'}  =  /  /  (^  -  7}r)k(y  ~  i)y)'PxY (a1, y)drdy  (2.13) 

J —ocJ  —  oo 

The  second  order  central  moments  are  of  interest  here.  These  are  /< 2o  =  o2x-  /'02  =  <?l 
and  // 1 1 ,  which  is  the  covariance  of  A'  and  Y,  given  by 

/'ii  =  E{(X  -  Vx)(Y  -  7/,)}  =  E{XY}  -  E{X}E{Y)  (2.14) 


Finally,  the  correlation  coefficient  of  A"  and  Y  is  defined  as 

Gw  =  £{(A:  -  1]X){Y  -7 Jy)}  _  Hu  _  /'ll  (0  1  M 

^/Eii^^EiiY-Vyy}  ~  VJwo~2  Wy 

The  correlation  coefficient  gives  a  measure  of  the  similarity  of  the  two  random  vari¬ 
ables.  It  will  have  a  value  of  1  when  the  random  variables  are  identical  (when 
A  =  V  );  a  value  of  0  when  there  are  no  similarities  between  one  random  variable 
and  the  other;  and  a  value  of  —1  if  one  random  variable  is  t  he  opposite  of  the  other 
(when  A'  =  -1"). 
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Central  Limit  Theorem 


Another  principle  to  be  applied  later  is  the  central  limit  theorem;  this  states  that  the 
joint  probability  density  function  of  a  sum  of  independent  random  variables  tends 
towards  a  Gaussian  distribution  as  the  number  of  random  variables  tends  to  infinity. 


Following  Goodmanl41h  if  U,  —  U\ ,  U2  •  •  •  Pyv  are  N  independent  random  variables 
with  means  of  77,  and  standard  deviations  of  a,,  then  the  random  variable  Z,  defined 
as 


z.'£Siz* 
fr{  <r, 

will  have  a  probability  density  function  of  the  form 

1 


PA  =  )  =  -s-e 


~P/2 


V/27r 


(2.16) 


(2.17) 


as  N  tends  to  infinity.  Note  that  this  does  not  depend  on  the  individual  probability 
density  functions  of  U,. 


The  central  limit  theorem  is  still  applicable  with  varying  degrees  of  validity  for  finite 
N. 


2.1.3  Random  Processes 

Definition 

The  definition  of  a  random  process  is  similar  to  that  of  a  random  variable,  the  differ¬ 
ence  being  that  where  a  random  variable  represents  values  with  certain  probabilit  ies, 
a  random  process  represents  functions  with  certain  probabilities.  The  definition  is 
again  based  on  an  underlying  random  experiment  and  the  possible  outcomes  of  this 
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experiment,  but  here  the  outcomes  are  functions.  A  random  process  is  all  possi¬ 
ble  output  functions  of  a  random  experiment  associated  with  their  probabilities  of 
occurring. 

The  functions  need  not  be  previously  defined;  an  observation  of  a  random  process 
(e.g.,  noise  in  a  system  with  respect  to  time),  will  in  general  yield  one  sample 
function,  x(t).  That  sample  function  is  now  known,  but  it  may  not  be  possible 
to  find  it  again  —  another  experimental  observation  will  generally  give  a  different 
sample  function. 

It.  should  be  noted  that  a  random  process  A'(/),  as  well  as  being  a  function  x(t) 
associated  with  each  event,  is  also  a  random  variable  A"  for  each  possible  (fixed)  t. 

The  probability  density  function  of  a  random  process  is  now  also  a  function  of  the 
same  parameters  as  the  random  process;  for  example,  if  the  random  process  is  A '(/>,  ry) 
then  the  (first  order)  probability  density  function  is  Pv(.r;p,c/).  The  moments  and 
central  moments  can  be  calculated  using  this  quantity,  and  will  also  be  functions  of 
p  and  ([. 

Stationarity  and  Ergodicity 

A  random  process  is  said  to  be  stationary  in  the  strict  sense  if  its  statistics  are 
independent  of  origin;  i.e.,  if  A'(/)  and  A(/  -f  s)  have  the  same  statistics  for  any  s. 
This  implies  that  first  order  statistics  are  constant,  and  second  order  statistics  only 
depend  on  the  difference,  .s.  Reference  to  stationary  quantities  later  in  this  work 
will  mean  strictly  stationary  ones. 


52 


CHAPTER  2.  MATHEMATICAL  CONCEPTS 


A  random  process  is  crgodic  if  an  average  over  one  sample  function  (a  temporal  or 
spatial  average  in  physical  systems)  is  identical  to  the  expected  value  of  the  process 
for  fixed  parameters  (its  ensemble  average).  This  means  that  the  statistics  of  a 
random  process  X(t)  can  be  determined  by  either  one  function  x(i )  or  by  one  cf  the 
random  variables  .Y  at  fixed  t.  Any  ergodic  process  is  also  stationary. 

Moments  of  a  Random  Process 

The  mean  of  a  random  process  . V ( / )  for  fixed  i  is  the  expected  value  of  one  instance 
of  that  random  process, 

»,(/)  = /?{.Y(/)}  (2.18) 

Notice  that  it  is  a  function  of  /. 

Other  quantities  of  interest  arc  the  second  order  moment  /v{A'(/)2},  and  the  second 
order  central  moment  or  variance 

<r(<f  =  E{(X(t)-v[<))2}  (2.19) 

The  joint  moments  of  two  instances  of  a  random  process  .Y(/[)  and  A  (/ 2),  are  from 
equation  (2.9), 

-  £{A'(f,)*.Y(/2)'}  (2.20) 

nii\[tiih)  is  the  autocorrelation  function,  R.xxHiH 2),  of  A(/). 

The  joint,  central  moments  of  ‘wo  instances  of  a  random  process  A’ ( / ! )  anti  A  (t2), 
are  from  equation  (2.13), 


=  £{(A'(M  -  7, (/,))*( X{h)  -  vUz))1} 


(2.21) 
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is  the  autocovariance  of  X(t), 

Equation  (2.21)  leads,  following  equation  (2.15),  to  the  correlation  coefficient  be¬ 
tween  two  instances  of  a  random  process  A(G)  and  -V(<i): 

CxxiU  0)  =  V(ti))(X(ti)  ~  v{h))} _  __  finUiiO)  ,9  99. 

s/MxUFSC) mpfcp*?!  <’(uMh) 

2.1.4  Differentiation  of  a  Random  Process 


Papoulisl40!,  deiives  a  few  important  relations  for  differentiating  random  processes. 
First  the  expected  value  of  the  derivative  X '(t)  of  a  random  process  .V ( / )  is  equal 
to  the  derivative  of  the  expected  value  of  that  process: 

E{X\t)}  =  j(E{X(t) } 

The  autocorrelation.  R_\  \  —  mu.  of  the  derivative  of  the  process  is 

S-R.x.x 


H.x'.x'UiJi) 


(2.24) 


dtidt2 

If  A  (/)  is  stationary,  then 

j  ix  '-V'  (/ )  =  -  =  ~  K\  A-  ( 0  ( 2. 25) 

Finalh,  note  that  when  1  =  0,  the  two  random  variables  are  the  same  and  Rx x  is 
the  second  moment  of  A;  i.e.. 


/:’{[-V'(/)ni(=o  =  V?.v'A-(0)  =  ~/?X.v(0) 


(2.20) 


2.2  Gaussian  Random  Rough  Surfaces 


Gaussian  surfaces  will  be  considered  here  since  an  assumption  of  Gaussian  statis¬ 
tics  (as  will  be  seen  later)  eases  analytical  calculations,  and  the  surfaces  can  be 
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manufactured  physically  and  also  created  computationally.  Here,  a  mathematical 
description  of  a  Gaussian  surface  is  given  and  some  of  its  properties  analysed. 


The  Gaussian  random  rough  surface  shall  be  represented  by  a  stationary  random 
process  S(x,y)  with  a  zero  mean,  E{S(x,y )}  =  0.  The  probability  density  function 
of  the  process  (i.e.,  the  height  J*stribution  of  the  surface)  is  Gaussian  in  form  with 
standard  deviation 


Ps(s)  = 


y/TKOS 


t  2°s 


(2.27) 


Also,  a  Gaussian  correlation  function  shall  be  imposed  on  the  surface: 


CSs\t)  = 


(2.2S) 


i.e.,  the  surface  has  a  correlation  length  r. 

Since  the  process  has  a  zero  mean  its  autocorrelation  is,  from  equation  (2.22), 

llssV)  =  (mV)  =  4Css(t)  (2.29) 


Later,  the  curvature  distribution  of  the  surface  will  be  required.  This  is  not  difficult 
to  obtain,  the  distribution  of  the  derivative  of  a  Gaussian  process  being  Gaussian 
itself  (this  follows  from  the  fact  that  any  linear  combination  of  Gaussian  tandom 
variables  ---  such  as  differentiation  —  Is  itself  Gaussian).  All  that  is  required  is  to 
find  the  standard  deviation  of  (Ins  new  Gaussian  distribution  and  the  correlation 
function  of  the  derivative. 

Equation  (2.25)  givi  ?  the  correlation  function  of  tne  derivative  distribution  (from 
equation  (2.29))  and.  again  since  the  mean  is  zero,  equation  (2. 26),  gives  the  variance 
of  the  derivative  of  the  surface. 
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The  autocorrelation,  therefore,  of  the  first  derivative  of  the  surface  height  (the  gra¬ 
dient),  is 

/  9/2  \  9/t2,  ,3 

(2.30) 


/WM  =  (l  -  %) 


The  surface  gradient  has  a  Gaussian  probability  density  function  (cf.  equation  (2.27) 
with  zero  mean  (by  equation  2.23))  and  standard  deviation 

V2os 


&S'  —  \J  0)  — 


(2.31) 


Applying  this  again,  the  autocorrelation  of  the  second  derivative  of  the  surface  height 
(the  curvature),  is 


(  4t 2  \  l‘2cr| 

B 

(’  -  e  +  30  J  e  £  75 

<3.33,  | 

So,  the  curvature  of  a  Gaussian  random  rough  surface  has  a  zero  mean  Gaussian 
distribution  with  a  standard  deviation  of 


<*S"  =  \J  Ns"S"{  0)  = - 

T 


2  y/Z<TS 


(2.33) 


2.3  Convolution 


Of  interest  lat<*r  will  be  the  convolution  of  two  functions  </(.r)  and  /(.r),  generally 
expressed  as  h(.i)  =  /(.?•)  0  </(•>')  and  written 

h(.v)  =  /  f(u)(/(x  —  ii)du  —  [  f(x  —  u)y(u)du  (2.34) 

J-<x,  J -(x.  ’ 

Tiie  convolution  determines  the  area  of  the  product  of  the  overlapping  region  of  the 
two  functions  lor  the  various  displacements  x  of  one  with  respect  to  the  other. 


^I  '  .'.'V  .:  i: 
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2.3.1  Convolution  Theorem 

The  convolution  theorem  states  that  if  F(u)  and  G(u)  are  the  Fourier  transforms  of 
two  functions  f{x)  and  g{x),  then  the  product  of  these  two  transforms  is  the  Fourier 
transform  of  the  convolution: 


/t(.r)  =  f(x)  0  g{x)  (2.35) 

II(u)  =  I?{u)G(u)  (2.36) 

The  proof  of  this  follows  simply  by  taking  a  Fourier  transform  of  both  sides  of 
equation  (2.34). 

2.3.2  Extension  of  the  Convolution  Theorem  to  Correla¬ 
tion  Coefficients 


With  reference  to  equation  (2.22),  if  .Y(/)  has  a  zero  mean  (?/(/.)  —  0  for  all  t :),  then 
the  correlation  coefficient  may  be  written: 


£{.V(/,)..V«2)} 


f  .'x  O'-  \ 


t-->q 


Now,  if  ,Y{/)  is  stationary,  i.e.,  U  —  t  and  t2  —  /  +  and  moreover  ergodic  —  such 
that  the  expected  value  may  be  rewritten  as  a  time  average  over  a  sample  function, 
the  two  being  identical  —  then: 


Cax(*)  =  lb» 

7  —"X' 


■£rflTr{0*'{t  +  s)<lt 
Trflr  {*(t)}*dt 


(2.3S) 


The  1/27’  term  divides  through,  and  the  limit  may  be  realised. 


C.v.y(.s) 


+  ,)rlt 

JTJr(t)Y<I( 
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# 


# 


Finally,  if  x(t)  is  normalised  such  that  the  integral 

I”  [x(t)]2dt  (2.39) 

J  -OQ 

(sometimes  referred  to  as  the  total  energy  of  x(t))  is  equal  to  1,  then  the  correlation 
coefficient  becomes 

Cxx(s)=  [°°  x(t)x{t  +  s)dt  (2.40) 

Notice  that  this  is  very  similar  to  the  convolution,  equation  (2.34).  If,  as  with  equa¬ 
tion  (2.34),  Fourier  transform. ;  are  taken  of  both  sides  of  the  equation,  a  relation 
similar  to  the  convolution  theorem  should  be  expected.  In  fact,  for  two  member  func¬ 
tions  fix)  and  <y(.r)  of  a  zero  mean  ergodic  random  process  (normalised  as  required 
above),  with  Fourier  transforms  F{u)  and  G(u),  the  transform  of  the  correlation 
coefficient.  C(n)  is 

C(v)  =  Fm(u)G{u)  (2,11) 

Where  F‘{v)  is  the  complex  conjugate  of  F(u). 

Equation  2.41  will  be  used  to  calculate  the  coi relation  coefficient  between  two  ex¬ 
perimental  speckle  pattern  images  (see  chapter  three). 

2.4  Electromagnetism 

2.4.1  The  Helmholtz-Kirchhoff  Integral 

The  derivation  below  follows  that,  by  Born  and  W'olf^  and  BeckmamJlyh  The 
time  dependence  of  various  quantities  has  been  suppressed  as  this  work  deals  only 
wit  h  moiiochroiuat  ic  light  and  does  not  consider  s])ecklo  pat  terns  which  change  wit  h 


respect,  to  time. 
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The  starting  point  will  be  Green’s  second  theorem, 

IIL***  -  ^ EdV  =  JJsETn  -  (2'42) 

where  the  operator 

•S-  —  n.V  (2.43) 

an 

and  E  and  x[>  are  solutions  U  to  the  wave  equation 


(V2  +  k2)U  =  0 


which  exist  in  the  volume  V  surrounded  by  the  surface  S  which  has  an  inward 
pointing  normal  n.  These  solutions  mean  that  the  volume  integral  in  equation 
(2.42)  will  he  zero,  so 


c hi' 
On 


3F 

-  4°±dS  =  0 
an 


(2,15) 


One  of  the  solutions  to  equat  ion  (2.44),  /s,  is  a  component  of  the  electric  field,  and 
the  other,  V\  satisfies  some  conditions  for  the  geometrical  propagation  —  in  the  case 
to  be  considered  this  is  the  spherical,  or  free  space,  Green’s  function: 


_  ik.r 


(2.16) 


The  standard  technique  for  integrating  over  a  function  with  a  singularity  is  used, 
namely  that  of  surrounding  the  singularity  by  a  small  sphere  with  a  radius  tending 
to  zero,  the  surface  of  this  sphere  connecting  with  that  of  the  main  bulk  by  an 
infinitesimally  thin  ‘corridor’.  This  means  that  the  singularity  (here  for  r  =  0), 
is  excluded  from  the  volume  contained  within  the  surface.  Figure  2.1  shows  this 
process,  integration  of  equation  (2.45)  should  now  be  performed  over  S  and  S i  (and 
along  the  corridor  joining  them,  bot  h  on  the  way  to  Si  and  from  it,  but  the  integrals 
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Figure  2.1:  Integration  over  a  surface  for  a.  function  with  a  singularity  at  point  P. 


along  this  corridor  cancel  each  other  and  so  can  be  omitted).  The  integration  now 
becomes 


Ii.E?L  -  +%JS  +  JL/£  -  -  °  (2'47) 

The  second  term  of  equation  (2.47)  can  be  shown,  substituting  in  equation  (2.46 
with  r  =  R  and  since  .S'i)  is  spherical,  to  be  —4nE(P),  in  the  limit  as  /?  tends  to 
zero.  E{P)  denotes  the  electric  field  at  the  point  P.  Finally  equation  (2.47)  gives 


OE 


M’ 


OE 


=  rJI/i 


OE 


(2.18) 


This  is  the  Ilelmholtz-KirchhofT  integral.  It  gives  the  electric  field  at  a  point  P  given 
that  it,  and  its  derivative,  are  known  along  the  surface  S.  It  will  bo  the  starting 
point  for  calculations  later  on. 


If  S  is  not  a  closed  surface,  equation  (2.48)  will  still  apply,  as  a  closed  surface  S', 
which  incorporates  S  and  encloses  all  sources,  may  be  considered.  The  reflectivity 
of  the  portion  of  S'  which  does  not  coincide  with  ,S  will  be  zero,  and  hence  no 
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contribution  could  possibly  be  made  by  it  to  the  field  at  a  point  P  within  it. 

2.4.2  Polarisation 


Figure  2.2:  Polarisation  state  definitions. 


The  polarisation  state  of  radiation  is  defined  by  the  behaviour  of  the  electric  vect  or  as 
the  field  propagates  through  space.  In  the  case  where  the  electric  vector  is  fixed  with 
respect  to  position/time  (i.e..  it  does  not  rotate  about  its  direction  ol  propagation), 
the  field  is  said  to  be  linearly  polarised. 

Figure  2.2  shows  a  ray  of  linearly  polarised  light,  wave  vector  k,,  striking  a  surface. 
The  electric  vector  of  the  ray  is  defined  within  the  e~,e+  plane;  this  plane  is  normal 
to  the  wave  vector.  The  plane  of  incidence  in  the  figure  is  defined  by  k,  and  the 
normal  to  the  scattering  plane,  in  this  case  the  c-axis.  If  the  electric  vector  is 
lying  along  the  e~  direction,  perpendicular  to  the  plane  of  incidence,  'lie  light  is 
s-polarised;  if  it  is  lying  along  the  e+,  parallel  to  the  plane  of  incidence,  it  is  p- 
polarised.  Any  other  linear  polarisation  state'  may  be  expressed  as  a  superposition 
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of  these  two  orthogonal  states.  The  superscripts  —  and  -f  are  used  to  represent  a 
quantity  applicable  to  the  s-  and  p-polarisation  states  respectively. 

2.4.3  Reflection  Coefficients 


The  Fresnel  reflection  coefficients,  71,  are  defined  here  as  a  multiplicative  factor  on 
an  incident  wave,  upon  reflection  from  a  flat  plane,  to  obtain  the  reflected  wave. 
These  reflection  coefficients  are  polarisation  dependent  such  that 

£+  =  71+E+ 

e-  =  n- e: 

and  may  be  complex. 

They  are  expressed  as^d; 

cos  Oi  sin  Oi  —  cos  0 sin  <p 
cos  0,  sin  0,  -f  cos  <?sin  0 
cos  0,  sin  0  -  sin  0 ,  cos  0 

'  f  \  •  j  7  •  n  1  (  ~ .  1 J  J 

COS  U,  SI »  0  -f  SI  II  Ui  cos  0 

where  0,  and  o  are  the  angle  of  incidence  and  transmission  respectively,  as  shown  in 
figure  2. ft.  Thc\  -He  directly  applicable  at  a  dielectric  interface. 


For  a  metal  surface,  equations  (2.49)  can  still  be  applied  if  the  “transmission  angle” 
0  is  redefined.  In  general,  for  a  metal  surface  with  complex  refractive  index  ?? ( 1  +  ?/») 
and  light  incident  at  an  angle  of  0X,  the  values  of  sin  0  and  cos  0  in  equation(2. 19) 
are  given  byM: 


sin  0 


1 


IK 


71(1  +  K2) 


~  Hill  Oi 


qt  n  =■-  9 (cos  7  +  i  sin  7) 


COS  0 
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Figure  2.3:  The  geometry  used  for  defining  the  Fresnel  reflection  coefficients. 


with  q  and  7  being  gi\en  by 


4  (  (1  -  K2)  .  2  A  (  2k 


sin2  0, 


tan  27  = 


2k  sin2  0, 


7i 2 1 1  f  k2)2  —  (1  —  k7)  sin2  0, 


(2.50) 


For  a  perfect  conductor,  the  reflection  coefficients  reduce  to  a  very  simple  form: 

7v+  =  1 

n~  =  -1  (2.5i) 


2.4.4  The  Kirclihoff  Approximation 

The  KirchhofT  (or  Physical  Optics)  method  is  a  technique  for  solving  equation  (2.18), 
the  Ilclniholtz-KirchhofF  integral.  It  consists  of  approximating  the  surface  field  (the 
sum  of  both  the  incident,  and  scattered  fields  at  the  surface)  and  its  normal  derivative 
at  a  point  on  the  rough  surface  to  those  given  by  the  tangent  plane  at  that  point.. 
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This  means  that  the  scattered  field  is  simply  the  incident  field  multiplied  by  the 
reflection  coefficient.  The  approximation  may  be  written  as 


Ei,  =  (1  + 

9E -  =  (i-n}dE’ 


dn 


dn 


where  Ei  is  the  incident  field,  and  E{S  is  the  field  at  the  surface. 


Throughout  this  work,  reference  to  the  Kirchhoff  approximation  strictly  relates  to 
the  definition  above,  nol  a  single  scatter  approximation  (to  be  discussed  in  chapter 
four)  which  usually  goes  hand  in  hand  with  it. 


Validity  of  the  Kirchhoff  Approximation 


Figure  2.4:  The  tangent  to  a  surface  feature  with  (a)  large  and  (b)  small  radius  of 
curvature  compared  with  the  illuminating  wavelength. 

Obviously,  the  Kirchhofr  approximation  will  be  good  when  the  radius  of  curvature 
of  a  surface  feature  is  large  compared  with  the  illuminating  wavelength  (see  figure 
2.4),  or  identically,  the  curvature,  C,  is  small  in  comparison  with  the  reciprocal  of 
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the  wavelength: 


(2.52) 


Beckmann^  and  OgilvyM  state  two  criteria  which  are  slight  variations  upon  equa¬ 
tion  (2.52).  Both  are  derived  from  geometrical  considerations;  the  first  is 


C  « 


47T  COS  4> 

A 


(2.53) 


where  <f>  is  the  local  angle  of  incidence  on  a  tangent  plane.  Unless  <f>  is  greater  than 
85  degrees,  this  inequality  will  support  greater  curvatures  than  will  equation  (2.52). 
The  second  variation, 

C<\co**0  (2.51) 

again  weights  the  inequality  in  favour  of  larger  curvat  ures  until  the  angle  of  incidence, 
0 ,  is  greater  than  47  degrees. 


For  an  understanding  of  where  these  inequalities  come  from,  the  reader  is  referred 
to  section  4.2.2  in  reference  [1].  The  purpose  of  presenting  them  here  is  to  show  that 
equation  (2.52),  to  be  used  below,  is  sufficient,  as  limitation  on  the  applicability  of 
the  Kirchhoff  approximation. 


The  reason  for  expressing  the  inequalities  in  terms  of  C  is  that  the  curvature  dis¬ 
tribution  of  a  Gaussian  random  rough  surface  h.r--  already  been  determined;  it  is 
Gaussian,  has  a  zero  mean  and  a  standard  deviation  of 

2y/la 


r- 


(2.55) 


The  standard  deviation  will  determine  the  proportion  of  the  surface  curvatures  which 
are  large  enough  to  invalidate  equation  (2.52). 
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In  what  follows,  an  attempt  is  made  to  gain  some  idea  of  the  surface  parameters 
(standard  deviation  and  correlation  length)  which  are  allowed  within  the  Kirchhoff 
approximation.  It  is  important  to  note  that  it  is  the  approximation  of  a  surface 
point  to  its  tangent  plane  which  is  being  examined  here,  not  the  validity  of  the 
application  of  the  single  scattering  model  (which  will  be  examined  in  more  detail  in 
chapter  four).  The  dependence  of  the  inequalities  on  the  angle  of  incidence  will  be 
ignored,  as  this  only  has  a  noticeable  effect  when  large. 


The  percentage  of  curvatures  which  lie  within  equation  (2.52)  is  easily  determined. 
Equation  (2.33)  is  substituted  in  to  the  limit,  of  equation  (2.52);  the  number  of 
standard  deviations  available  within  the  Kirchhoff  approximation  for  a  particular 
surface  is  then 


N  = 


(2.56) 


2\/3crA 

which  is  easily  converted  into  the  percentage  of  curvatures  lying  within  equation 
(2.52)  by  noting  that,  the  probability  distribution  for  a  Gaussian  probability  density 
function  is  the  error  function.  The  fraction  of  curvatures,  /,  within  the  limit  of 
equation  (2.52)  is  tlierelore 


/  =  2  erf(AT) 


(2.57) 


Figure  2.5  shows  the  limits  of  equation  (2.52)  for  boundaries  of  N  —  3, 2.5, 2, 1.5 
and  1  (99.7%.  98.8%,  95.5%,  86.8%  and  68.3%.  of  curvatures  respectively  within  the 
limit). 


As  can  be  seen  from  the  figure,  within  the  light  region  (to  all  ini  ruts  and  purposes), 
the  whole  distribution  of  surface  curvatures  obeys  inequality  (2.52);  in  fact,  the 
majority  of  these  curvatures  will  be  close  to  zero  (the  mean)  and  well  within  the 
limit. 
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Figure  2.5:  Regions  of  validity  of  the  KirchhofF  approximation  for  Gaussian  surface 
parameters.  The  percentages  represent  the  proportion  of  surface  curvatures  which 
lie  within  the  bounds  of  the  limit  to  the  approximation  as  given  by  equation  (2.52) 

Although  there  is  no  clear  line  at  which  a  boundary  between  applicability  and  non¬ 
applicability  may  be  defined,  the  KirchhofF  approximation  can  be  applied  with  con¬ 
fidence  to  any  Gaussian  surface  within  the  light  region.  It  is  likely  to  meet  less 
success  with  surfaces  closer  to  the  darker  region,  fairing  worst  near  to  the  a  axis. 

The  dependence  of  applicability  on  the  angle  of  incidence  will  only  manifest  itself 
with  those  surfaces  which  lie  close  to  the  boundary  of  the  ‘good’  region  of  the  figure, 
and  then  only  when  the  angle  of  incidence  is  large. 


Chapter  3 


Experimental  Procedures 


3.1  Random  Rough  Surfaces 


All  of  the  surfaces  to  be  studied  here  will  be  two-dimensional,  i.e.,  the  roughness  of 
the  surface  will  not.  depend  upon  direction  —  it  will  vary  as  a.  function  of  both  ,r  am’ 
y.  They  will  also  have  Gaussian  height  distributions  and  correlation  functions  — -  the 
reasons  for  this  are  that  it  is  possible  to  create  a  surface  with  known  Gaussian  statis¬ 
tics;  it  is  easy  to  work  theoretically  with  such  surfaces  and  that  there  already  exists 
a  selection  of  Gaussian  surfaces  whose  parameters  are  known  and  whose  scattering 
properties  have  been  investigated^12].  One-dimensional  surfaces  arc  not  considered 
in  this  work  as  they  do  not  produce  the  required  two-dimensional  speckle  patterns. 


3.1.1  Surface  Manufacture 


In  order  to  gain  a  measure  of  control  over  the  statistical  properties  of  the  surfaces, 
a  method  developed  by  Gray]15)  was  used.  This  method  involves  coating  a  glass 
substrate  with  a  layer  of  a  photoresist  which  has  an  approximately  linear  height  re- 
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sponse  to  exposure  time.  The  coated  substrate  is  then  exposed  to  speckle  patterns 
created  by  illuminating  ground  glass.  The  surfaces  used  for  experimental  measure¬ 
ments  presented  in  this  work  were  not  manufactured  by  the  author,  but  were  chosen 
to  enable  comparisons  with  existing  published  work. 


Preparing  the  Substrate 

The  aim  here  is  to  deposit  a  flat  layer  of  photoresist  onto  the  glass  substrate  (ap¬ 
proximate  dimensions  of  50  x  50  x  '3mm),  with  sufficient  thickness  to  allow  a  surface 
with  an  rms  height  of  at  least  2 /.mi  to  remain  after  the  necessary  etching.  This  turns 
out  to  be  an  initial  coating  around  ,U fim  thick. 

The  first  stage  to  be  considered  will  be  cleaning  the  glass  substrate;  this  is  important, 
not  only  for  improving  the  adhesion  of  the  photoresist,  but  also  to  remove  any  traces 
of  dirt  or  chemical  impurities  which  may  affect  the  final  coating.  Initially,  the  glass 
plates  are  washed  in  hot  water  with  a  detergent,  to  remove  grease  and  dirt.  Next,  the 
plate  is  cleaned  with  acetone  to  remove  any  chemical  deposits  on  the  surface  —  some 
of  which  will  be  left  by  the  detergent.  Finally,  the  plate  is  wiped  with  isopropanol. 
This  last  stage  is  necessary  since  it  was  found  that  acetone  leaves  deposits  which 
adversely  affect  the  photoresist,  producing  striations. 

The  photoresist,  used,  which  was  found  to  have  a  good  linear  response  to  the  exposure 
tirru.,  is  Shipley  Sl-100-37.  'This  photoresist  is  also  sufficiently  viscous  to  proems 
a  thick  coat  on  the  substrate.  The  method  of  spin  coating  is  used  to  deposit  the 
photoresist  onto  the  substrate;  the  figures  given  below  were  found  empirically,  and 
produce  a  uniform  coat  of  photoresist  of  between  JO  and  12  /tm  thickness. 


* 

* 


m 


m 

* 
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A  cleaned  substrate  is  fastened  to  the  centre  of  a  turntable  set  up  to  spin  at  300 
revolutions  per  minute.  While  the  turntable  is  stationary,  photoresist  is  applied  to 
the  centre  of  the  substrate  with  a  syringe,  until  the  natural  spread  of  the  liquid 
reaches  the  edges  of  the  glass,  about  3 ml.  The  turntable  is  now  set  spinning  for 
four  minutes.  Again  with  the  substrate  stationary,  a  similar  dose  of  photoresist  is 
applied.  Finally  the  turntable  is  spun  for  a  further  four  minutes. 


After  allowing  twenty  four  hours  to  dry  at  room  temperature,  the  coated  substrates 
are  baked,  at  ninety  degrees  celcius  (this  is  below  the  temperature  at  which  the 
photoresist  melts)  fui  thirty  minutes,  to  drive  on  any  remaining  sol  vein  in  the  pho¬ 
toresist  layer.  It  is  now'  ready  to  have  a  surface  profile  imposed  upon  it. 


A  few  factors  need  to  be  considered  during  the  manufacture  of  tlm  coated  plates. 
In  order  to  produce  a  perfectly  flat  layer  of  photoresist,  air  bubbles  must  not  be 
allowed  in  the  photoresist;  as  long  as  care  is  taken  with  the  syringe  then  the  process 
outlined  above  will  ensure  this.  Dust  must  not  be  allowed  onto  the  plates  at  any 
stage  during  manufacture;  one  speck  of  dust  can  produce  a  relatively  large  area 
which  will  not  react  properly  to  exposure  —  this  is  especially  a  problem  during  the 
twenty  four  hour  drying  period,  any  dust  settling  on  the  surface  during  this  time 
will  stick,  and  after  exposure,  a  spike  can  remain  sticking  up  above  the  surface.  It 
is  possible  for  dust  to  settle  on  the  photoresist  near  the  edges  of  the  plate  without 
invalidating  its  use  as  a  surface.  Finally,  care  must  be  taken  not  to  scratch  the 
photoresist,  compressed  air  should  be  used  to  remove  any  dust  which  settles  after 
drying. 
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Etching  the  Surface  Profile 


The  required  surface  profile,  as  mentioned  previously,  is  Gaussianly  distributed. 
Fortunately,  a  profile  of  this  approximate  nature  can  be  achieved  by  simply  exposing 
the  photoresist  to  eight  uncorrelated  speckle  patterns,  the  number  being  found  as 
follows. 


It  can  be  shown  that  the  probability  density  of  the  intensities  in  a  speckle  pattern 
with  mean  intensity  y  obeys  negative  exponential  statistics^, 


P/(D  ~  -  exp 
V 


(3.1) 


and  that  the  (unnormalised)  probability  density  function  of  a  superposition  of  ;Y 
uncorrelated  speckle  patterns,  each  with  mean  intensity  //,  is 


PRI)  = 


jN~l  ]V N 


(N  -  l)!//v 


7 


(3.2) 


The  variance  for  this  intensity  distribution  is 


Repeated  integration  by  parts  is  necessary  to  find  this. 


Figure  3.1  shows  the  form  of  this  distribution  for  various  N.  Eight  speckle  patterns 
are  taken  as  the  optimum  number  to  achieve  a  Gaussian  profile  with  the  maximum 
possible  variance  —  the  number  at  which  a  Gaussian  curve,  with  a  variance  equal 
to  that  given  by  equation  (3.3),  is  most  similar  in  shape  to  the  probability  density 
function,  equation  (3.2).  Note  that  as  N  increases,  the  profile  becomes  more  Gaus¬ 
sian  (cf.  central  limit  theorem),  but  the  variance  of  the  distribution  reduces,  thus 
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Figure  3.1:  The  probability  density  distribution  of  the  intensity  of  the  superposition 
of  N  uncorrelatcd  speckle  patterns,  each  with  mean  r;  =  1,  for  various  N . 


reducing  the  rins  surface  height  —  a  quantity  potentially  required  to  be  as  large  as 
possible. 

These  eight  uucorrelated  speckle  patterns  are  created  by  illuminating,  in  turn,  sep¬ 
arate  plates  of  ground  glass  with  the  4.57 Nam  wavelength  of  light  from  an  argon  ion 
laser.  The  photoresist  used  is  particularly  responsive  to  wavelengths  approaching 
the  ultraviolet  region  of  the  spectrum. 

The  geometrical  arrangement  for  the  exposure  of  the  coated  substrates  is  given  in 
figure  3.2.  Collimated  light  is  incident  on  the  ground  glass  diffuser  in  an  attempt 
to  produce  a  uniform  field  at  the  surface  —  important  to  give  a  uniform  statistical 
variation  to  the  surface  height.  The  size  of  the  illuminated  area  of  the  diffuser,  and 
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the  distance  of  the  substrate  from  it.,  determine  the  speckle  size  on  the  photoresist  — 
this  determines  the  correlation  length  (the  1/e  point  of  the  autocorrelation  function) 
of  the  surface. 


Figure  3.2:  Arrangement  for  exposure  of  photoresist  coated  plates  during  surface 
creation. 


For  a  more  precise  idea  of  how  the  arrangement-  will  affect  the  correlation  length  of 
the  surface,  a  simple  statistical  study  is  necessary.  Given  that  the  beam  incident-  on 
the  ground  glass  plate  has  a  Gaussian  cross-section,  then  the  correlation  between 
two  points  of  the  speckle  pattern  isf'1! 


C[{Ar)  =  exp 


\Z2Ku>Nr 


\z 


(3-4) 


where  it  has  been  assumed  that,  the  speckle  pattern  is  statistically  stationary.  Ac  is 
the  distance  between  two  points  of  the  speckle  pattern,  and  u>  is  the  1/e  half  width 
of  the  Gaussian  beam  profile  at  the  ground  glass  plate.  A  is  the  wavelength  of  the 
light,  and  ^  is  the  distance  of  the  surface  from  the  diffuser. 


From  equation  (3.4)  if  can  be  seen  that,  the  correlation  length  of  a  speckle  at  tin 
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surface,  and  hence  the  correlation  length  of  the  surface,  will  be 


A  z 


(3.5) 


The  parameters  c  and  u>  in  the  arrangement  can  be  adjusted  until  r  is  of  the  required 
size.  Generally  7  will  be  set  such  that  it  is  just  large  enough  for  the  speckle  field 
to  appear  uniform  over  the  substrate  (about  400mm),  anything  larger  would  just 
reduce  the  amount  of  power  incident  on  the  photoresist  and  hence  increases  exposure 
time.  This  leaves  a.’  as  the  main  handle  on  r,  but  because  it  is  not  possible  to  measure 
the  1/e  half  width  of  a  Gaussian  beam  quickly  and  prec  isely,  especially  for  large  r 
when  u>  is  small  (if  the'  required  r  is  150 fan  then  tv  is  approximately  0.3mm  —  for 
the  wavelength  and  distance  specified  above),  r  can  only  be  roughly  controlled.  Fine 
tuning  is  readily  achieved  by  adjusting  c  slightly. 


The  mis  height  of  the  surface  can  be  controlled  by  exposure  time;  the  maximum  time, 
before  the  photoresist  is  removed  down  to  the  glass  substrate,  has  been  determined 
experimentally  for  substrates  prepared  using  the  procedure  described  above,  and  is 
given  in  minutes  by 


(3,6) 


where  I  is  the  average  light  intensity  in  the  plane  of  the  substrate  in  the  arbitrary 
units  given  by  the  “ Berry  meter".  The  number  275  was  found  experimentally,  and 
converts  the  intensity  readings  given  by  the  meter  into  minutes  of  exposure  time. 
The  finally  calculated  exposure  time  should  be  divided  by  eight  and  each  speckle 
pattern  should  be  projected  onto  the  photoresist  for  this  time.  This  maximum 
exposure  will  give  the  surface  an  ruts  height  of  approximately  2//m. 


After  all  exposures  are  complete,  the  surface  is  developed  in  a  mixture  of  one  part  by 
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volume  of  Shipley  AZ-303  to  five  parts  water  for  thirty  seconds.  Finally,  the  surface 
is  rinsed  under  running  water  for  about  three  minutes,  and  dried  with  compressed 
air.  A  quick  check  should  be  made  to  ensure  that  the  glass  substrate  has  not  been 
exposed  through  the  photoresist  at  any  point. 

3.1.2  Gold  Coating 

The  last  stage  in  the  fabrication  of  a  surface  is  to  deposit  a  thin  layer  of  gold,  using 
vacuum  deposition  equipment,  onto  the  photoresist.  This  layer  provides  a  highly 
conductive  front  to  the  surface,  enabling  its  scattering  properties  to  be  compared 
with  theoretical  studies  which  assume  the  surface  is  a  perfect  conductor. 

The  complex  Fresnel  reflection  coefficients  for  gold  are  given  by  the  sequence  of 
equations  in  chapter  two,  starting  with  equations  (2.49).  The  square  magnitudes 
of  the  coefFicients,  7Z+  and  7Z~ ,  the  reflectivities,  are  shown  with  respect  to  angle 
in  figure  3.3,  the  refractive  index  of  the  gold  taken  at  a  wavelength  of  0.(>33//m  as 
0.167  -f3.149d42h  It  can  be  seen  that  the  curves  are  close  to  unity  and  approximately 
flat,  indicating  that  the  reflection  coefficients  are  only  loosely  dependent  on  the  angle 
of  incidence;  a  favourable  comparison  with  a  perfect  conductor  where  |7v|2  is  1  for 
both  polarisation  states. 

3.1.3  Surface  Characterisation 

Since  only  rough  control  over  the  statistical  parameters  of  the  surface  is  attainable, 
a  method  is  required  to  gauge  more  accurately  what  these  parameters  are.  For  all 
of  the  surfaces  mentioned  in  this  work,  a  Rank  Taylor-IIobson  Talystep  mechanical 


3.1.  RANDOM  ROUGH  SURFACES 


75 


IRI  2 


l.o  r 


0.8 


0.6 


0.4 


0.2 


0.0 


-J - L 


J _ L 


0  10  20  30  40  50 


_j _ i _ i _ i 

60  70  80  90 


6  j  /degrees  j 

Figure  3.3:  Fresnel  reflection  coefficients  for  gold  at  A  =  0.633/tm  over  a  range  of 
angles. 


surface  profilometer  was  used.  It  has  been  shown  that  this  instrument  does  not 
significantly  damage  the  surface1'*'^;  nevertheless,  measurements  were  made  away 
from  the  region  on  which  light  was  to  be  incident,  and  the  assumption  that  the 
statistics  are  uniform  oxer  the  whole  surface  was  made. 

The  profilometer  operates  by  dragging  a  stylus  across  the  surface;  the  stylus  available 
was  a  diamond  wedge  of  dimensions  1.8  x  0.5/im.  A  force  is  applied  to  the  stylus 
to  maintain  its  position  in  contact  with  the  surface,  and  this  force  varies  as  the 
height  of  the  surface  changes.  Generally  the  load  applied  to  the  stylus  should  be 
around  Gmg  at  the  middle  of  its  displacement  range,  this  will  mean  that  for  a  typical 
surface,  the  force  should  vary  between  about  2 mg  and  1  Qmg  at  the  extre  nes  of  the 
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range.  The  force  is  not  constant  as  this  is  the  quantity  which  is  measured  by  the 
instrument  and  converted  into  an  electrical  signal.  This  electrical  signal  is  amplified 
and  filtered  to  remove  high  frequency  vibrations,  finally  passing  to  an  analogue- 
to-digital  converter.  The  end  result  of  this  process  is  a  digital  representation  of  a 
2 mm  line  on  the  surface;  computational  processing  can  now  be  performed  to  find 
the  various  statistical  quantities. 


Processing  tL  rofilometer  Data 

X s  it  is  not  likely  that  the  surface  could  be  positioned,  by  hand,  exactly  level  under 
the  stylus,  the  trace  must  have  any  overall  trend  in  its  gradient  removed  before 
either  of  the  surface  parameters  can  be  obtained  it.  To  achieve  this,  a  straight  line 
is  fitted  to  the  data  by  the  least  squares  method,  where  the  surface  height  at  any 
position  along  that  line  is  the  distance  to  the  trace  along  a  perpendicular  to  the  line. 
In  other  words,  the  trace  is  shifted  to  the  origin  and  rotated  such  that  the  mean 
height  and  the  mean  gradient  are  both  zero.  The  gradients  removed  by  this  method 
can  be  smaller  than  0.001. 

The  standard  deviation  of  the  trace  is  now  easily  calculated.  This  will  be  the 
standard  deviation  (mis  height)  of  the  surface  over  the  length  of  the  scan.  Figure 
3.4  shows  a  histogram  of  the  distribution  of  heights  over  one  scan,  a  Gaussian  curve  is 
also  shown  with  the  same  standard  deviation.  It  can  be  seen  that  the  data  has  a  form 
which  is  roughly  Gaussian,  the  departure  being  due  to  the  long  correlation  length  of 
the  surface  presented.  The  actual  value  of  the  standard  deviation  attributed  to  the 
surface  will  be  the  average  of  many  such  individual  trace  values;  the  error  associated 
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with  it  will  then  be  the  standard  deviation  of  the  spread  of  values. 


Figure  3.4:  Histogram  of  the  surface  heights  over  a  single  scan  of  a  surface  (surface 
#80  —  see  later),  a  Gaussian  curve  with  the  same  variance  is  also  shown. 


The  correlation  length  of  the  surface  is  obtained  by  performing  an  autocorrelation  on 
the  trace.  The  result,  after  dividing  by  the  variance,  is  the  autocorrelation  function. 
This  will  be  approximately  Gaussian  in  form  (since  the  beam  profile  used  to  create 
the  surfaces  was  approximately  Gaussian),  at  least-  for  very  low  lag  (_V)  values: 


(  Ar2 

C'.s5(  Ar)  =  exp  ( - 
A  least  squares  fit.  of  this  will  give  the  correlation  length. 

\ 


(3.7) 


T  — 


£;=i(i  ny>? 


(3.S) 


Figure  3.5  shows  the  correlation  function  calculated  for  a  surface  trace.  Equation 
(3.7)  is  also  plotted  for  a  correlation  length  calculated  by  equation  (3.8;  for  the 
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_ _ _ Lag/ pm 

Figure  3.5:  The  correlation  function  of  a  surface  scan,  a  fitted  Gaussian  correlation 
function  has  a  correlation  length  of  19.5-1  /<m. 

portion  of  the  function  which  has  a  value  greater  than  0.1;  i.e.,  n  in  equation  (3.8) 
is  determined  by  the  point  at  which  the  data  first  drops  below  0.1  with  increasing 
lag.  This  is  necessary  because  as  at  higher  lags,  the  correlation  function  will  only 
go  to  zero  for  an  infinite  trace  length.  Again  the  parameter  for  the  surface  and  its 
error  are  obtained  by  averaging  many  such  correlation  lengths. 


Stylus  Size  Considerations 

For  a  surface  where  the  radius  of  curvature  of  the  features  is  greater  than  the  width 
of  the  stylus  tip,  it  can  be  assumed  that  the  size  of  the  tip  will  not  greatly  influence 
the  measurements  obtained.  However,  once  the  surface  features  are  of  a  similar 
size  to  the  stylus  used  to  trace  them,  the  tip  dimensions  will  affect  the  measured 
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surface  heights.  The  l.S/nn  x  0.5/nn  stylus  used  here  is  suited  to  the  analysis  of 
one-dimensional  surfaces,  and  does  not  introduce  significant  error  in  these  cases^3!; 
but  with  a  two-dimensional  surface,  the  extra  width  along  one  side  will  prevent  the 
stylus  from  properly  tracing  the  surface  portion  underneath  its  centre.  The  stylus 
will  in  fact  trace  the  highest  points  under  the  whole  of  its  width,  as  can  be  seen  in 
figure  3.6. 


Figure  3.6:  Finite  stylus  tip  size  prevents  accurate  tracing  of  surface  features  of 
a  similar  size.  The  scale  has  been  exaggerated,  generally  the  effect  is  much  less 
noticeable. 

The  errors  introduced  to  the  surface  parameters  have  been  examined^,  based  on 
an  assumption  of  a  circular  tipM.  Relations  derived  between  the  measured  surface 
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parameters  and  the  actual  values,  for  Gaussian  surfaces  with  a  Gaussian  correlation 
function,  such  as  those  studied  here,  are: 


(3.9) 

(3.10) 


where  the  subscript  m  represents  the  value  as  measured  by  a  stylus  tip  with  radius 
R  and  *y  is  the  rms  curvat  ure  of  the  surface: 

2\/3<t 

1  =  (3.11) 


as  derived  in  chapter  two. 


The  actual  effect  of  this  finite  tip  size  (1.8pm  across  the  width)  wiii  introduce  an 
error  of  less  than  two  percent  to  the  measured  parameters;  although  much  larger 
than  for  one-dimensional  surfaces,  it  is  still  swamped  by  the  statistical  uncertainty, 
even  for  the  roughest  of  surfaces  st  udied  here  (see  below). 


3.1.4  The  Surfaces  Studied 


Three  surfaces  will  be  used  in  the  experimental  investigation,  all  of  which  arc  ap¬ 
proximately  Gaussian  in  both  height  distribution  and  correlation  function.  The 
parameters  of  these  surfaces  are  given  in  table  3.1.  Surface  #80  has  previously 
been  reported  to  be  a  Beckman  surface,  i.c.,  it  exhibits  only  single  scattering!12!. 
Surface  #S3  is  one  which  shows  strong  enhanced  back-scatter,  a  multiple  scatter¬ 
ing  phenomenon!12!.  Both  of  these  surfaces  were  manufactured  by  Mendez  and 
O’Donnell.  Surface  #239  was  chosen  as  an  ‘in  between’,  it  is  expected  to  border  011 
the  realm  of  multiple  scattering  —  this  was  manufactured  by  Kim!41!. 
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Surface 

<7  /  fim  t  j  um 

#80 

#239 

#S3 

2.27  ±0.2  20.7  ±2.0 

1.34  ±0.1  7.0  ±0.5 
0.9  ±0.1  2.73  ±0.08 

Table  3.1:  Kms  height-  and  correlation  length  for  the  surfaces  to  be  studied. 

The  surfaces  have  been  plotted  on  a  graph  showing  the  validity  region  of  the  Kirch¬ 
hofF  approximation  (up  to  99.7%  of  curvatures  valid)  in  figure  3.7.  It  can  be  seen 
that  all  three'  surfaces  have  gradients  varying  sufficiently  slowly  for  the  KirchhofF 
approximation  to  be  expected  to  apply. 


Figure  3.7:  Surfaces  #8 0,  #239  and  #83  plotted  on  the  KirchhofF  validity  graph 
for  a  wavelength  of  A  =  0.633/t?u.  The  large  black  areas  represent  the  errors  in  the 
measurements  of  the  surfaces. 
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3.2  Obtaining  the  Angular  Correlations 
3.2.1  Experin  ital  Equipment 

The  equipment  used  is  shown  in  figure  3.8.  The  illumination  source  is  a  Uniplia.se 
lOmVU  Helium-Neon  laser  cm'Hing  linearly  polarised  light  at  a  wavelength  of  633 nm. 
The  orientation  of  the  laser  controls  the  polarisation  of  'he  light  incident  on  the 
surface.  The  speckle  patterns  are  digitised  with  a  Piinceton  Instruments  512  x  512. 
27 /mi  pixel,  sixteen  bit  pettier  cooled  CCD  camera  (accurate  to  one  bit  —  the 
readout  noise  is  about  10  photoelcctrons),  and  transfeired  to  an  magneto-optical 
disk  for  later  off  line  processing.  A  polariser  is  positioned  in  front  of  the  camera  to 
enable  control  over  which  component  of  the  scattered  light  is  detected.  No  imaging 
optics  are  used,  a  speckle  pattern  being  formed  at  the  CCD  element  purely  by  the 
free  space  propagation  of  light  from  the  surface. 

The  camera  was  not  actually  mounted  on  the  end  of  the  arm  as  shown  in  figure  3.8 
since  its  weight  would  have  put  an  unnecessary  stress  on  the  rig’s  bearings.  Instead, 
a  counter  weight  was  removed  from  the  arm,  and  the  camera  was  mounted  in  its 
place  much  closer  to  the  bearings  —  as  shown  m  figure  3.9. 

As  the  camera  is  now  mounted  behind  the  surface,  a  mirror  must  be  positioned  at 
some  length  along  the  arm,  angled  such  that  the  light  scattered  around  the  specular 
direction  is  sent  towards  the  camera.  'FI: e  camera’s  position  and  the  angles  tin. 
camera  and  minor  were  tilted  at  were  chosen  to  give  a  path  lengt.i  from  the  surface 
to  the  camera  of  around  900m m.  The  camera  was  also  high  enough  above  the 
surface  that  its  view  was  not  obstructed  by  the  surface  mountings.  Fine  control  of 
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angles. 

the  path  length  was  obtained  by  moving  the  mirror  along  the  arm  and  illuminating 
the  surface  at  a  height  which  would  mean  the  centre  of  the  camera  was  looking 
exactly  at  the  laser  spot. 

The  angles  of  the  mirror  and  camera  mounts  were  cut  to  a  tolerance  of  less  that  0.1 
degrees.  A  measurement  of  the  path  length  after  construction  gave  a  distance  of 
000  i  5mm. 

The  size  and  number  of  speckles  formed  at  the  CCD  element  is  an  import  ant  consid¬ 
eration.  Obviously,  a  large  number  of  speckles  in  an  image  mean  that  the  statistical 
properties  of  that  image,  such  as  the  mean  and  variance,  can  be  well  defined.  How¬ 
ever,  this  desirable  cri  erion  must  be  traded  against  the  need  for  each  individual 
.speckle  to  be  well  defined,  a  limitation  imposed  by  the  pixel  size  of  the  CCD  array. 
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Figure  3.9:  The  actual  arrangement  used  to  mount  the  camera.  The  camera  is  still 
looking  in  the  specular  direction  of  scattering  from  the  surface. 


The  speckles  iri  the  image  must  each  have  sufficient  intensity  information  in  order 
for  a  correlation  to  be  meaningful.  Speckles  with  a  ten  pixel  diameter  provide  this 
information,  while  still  allowing  around  400  speckles  over  the  whole  image. 


The  number  of  speckles  on  the  CCD  element  is  controlled  by  the  size  of  the  illumi¬ 
nating  spot  on  the  surface.  Since  the  detector  is  in  the  far  field,  the  speckle  size,  u, 
will  be  set  by  the  diffraction  limit  of  the  system,  roughly  given  by: 

1.22A  / 


u>  — 


D 


(3.12) 


the  aperture  diameter,  D,  being  the  beam  size  at  the  surface  and  the  focal  length,  /, 
being  the  distance  of  the  detector  form  the  surface  (in  this  case,  900mm  assuming 
the  beam  comes  to  a  focus  at  the  CCD  element).  If  radius  of  a  speckle  is  u>,  the 
area  of  a  speckle,  7ru.'2,  is  the  14  x  14mm2  CCD  element  area  divided  by  400  (the 
number  of  speckles  required},  about  0.5mm2.  The  required  beam  diameter,  on  the 
surface  from  the  above  equation  is  therefore  roughly  3 mm. 
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Note  that  the  CCD  element  subtends  an  angular  of  around  0.0  cleg  ices  in  both 
directions  to  the  surface  (solid  angle,  0.0002  steradians);  an  individual  speckle  has 
a  width  of  about  2.5  arc  minutes  subtended  from  the  surface. 

Contiol  of  tb  *  angle  at  which  the  incident  beam  strikes  the  surface  (the  orientation 
of  the  surface),  and  the  angle  at  which  the  camera  is  positioned  (the  orientation  of 
the  rig  arm),  is  by  means  of  two  stepper  motors.  These  provide  an  angular  resolution 
of  2'  (two  arc  minutes)  per  step.  Currently,  the  geometical  situations  attainable  with 
the  rig  are  limited  to  those  of  ‘in-plane’  scattering  —  the  surface  normal,  incident- 
beam  and  direction  o*-  observation  ail  lying  in  a  plan'’.  This,  however,  is  not  a 
problem  here  as  observation  i=.  confined  to  the  specular  direction. 

A  PC,  which  performs  many  duties  in  the  laboratory,  provides  the  interface  to  the 
stepper  motors.  The  Sun  workstation  controlling  the  experiment  communicates  with 
the  PC  via  a  serial  line. 

Due  to  the  nature  of  the  measurements  sought,  special  care  must  be  taken  to  ensure 
that  the  centu'  of  the  laser  beam  is  targeted  exactly  on  the  axis  of  rotation  of 
the  surface.  Failure  to  achieve  this  will  mean  that  at  different  angles  of  incidence, 
separate  portions  of  the  surface  will  be  illuminated. 

In  order  to  gain  the  maximum  advantage  offered  b\  the  camera,  the  integration  time 
is  raised  or  lowered  such  that  the  maximum  pixel  value  in  any  one  image  is  between 
50000  and  05000  (the  maximum  value  attainable  being  05585).  This  means  that  as 
the  intensity  changes  over  the  range  of  speckle  patterns,  the  camera  will  adjust  to 
utilise  most  of  its  intensity  resolution.  The  max. mum  intensity  the  camera  can  view 
is  limited  by  the  fastest  shutter  speed  (about.  lCm.s).  As  the  correlation  process 
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described  below  removes  any  dependence  on  the  overall  intensity  of  the  speckle 
pattern,  upsetting  the  relative  intensities  between  the  patterns  will  not  affect  the 
correlation  results. 


3.2.2  The  Speckle  Patterns 


Figures  3.10,  3.12  and  3.11  show  histograms,  taken  from  a  single  frame  of  the  raw 
camera  data  of  speckle  patterns,  obtained  from  surfaces  #S0,  #239  and  #83  at 
various  angles  of  incidence,  and  for  different  polarisation  arrangements. 


Figure  3.10:  Histogram  of  the  intensity  values  of  a  speckle  pattern  formed  by  surface 
#80  at  ten  degrees  incidence  and  observation  in  the  specular  direction.  Incident  and 
detected  light  are  both  s-polarised.  The  same  data  plotted  on  a  log  scale  is  inset. 


As  can  be  seen,  surface  #S0  and  surface  # S3  produce  speckle  patterns  that  have 
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Figure  3.11:  Histogram  of  the  intensity  values  of  a  speckle  pattern  formed  by  surface 
#83  at  fifty  degrees  incidence  and  observation  in  the  specular  direction.  Incident 
light  is  p-polarised  and  the  observed  pattern  is  s-polarised.  The  same  data  plotted 
on  a  log  scale  is  inset. 


roughly  the  negative  exponential  histogram  expected!45!,  but  that  the  first  few  bins 
are  considerably  less  than  they  should  be.  This  may  be  explained  by  noting  that 
a  large  amount  of  light  is  scattered  around  the  laboratory  by  the  rough  surfaces, 
and  that  this  will  be  registered  by  the  camera  as  an  ambient  light  level,  adding  to 
every  pixel  and  shifting  the  histogram  to  the  right;  i.e.,  zero  is  no  longer  the  most 
probable  intensity. 


The  quantity  cr//(7),  which  should  be  1  for  a  true  negative  exponential,  for  the 
speckle  pattern  in  figure  3.10  is  0.72.  This  increases  to  0.91  if  the  first  two  bins  of 
the  histogram  are  removed  from  the  calculation,  giving  evidence  that  ambient  light 
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is  the  distorting  factor.  This  may  be  shown  further  by  supposing  that  the  lowest  bin 
of  the  negative  exponential  should  contain  all  the  light  which  has  appeared  in  the 
two  lower  bins  —  adding  this  in  gives  a  07/(7)  —  0.98.  The  errors  on  these  figures 
may  be  considered  as  ±5%,  being  the  statistical  error  on  a  sample  of  400  speckles. 

For  figure  3.11,  crj/(I )  is  0.54,  but  increases  to  O.SS  after  subtraction  of  the  first 
three  bins.  If  the  subtracted  amount  T  added  into  what  was  the  fourth  bin  then  the 


Figure  3.12:  Histogram  of  the  intensity  values  of  a  speckle  pattern  formed  by  surface 
#239  at  thirty  degrees  incidence  and  observation  in  the  specular  direction.  One  data 
set  has  incident  light  s  polarised  and  observed  p-polarised,  the  other  has  incident 
light  p-polarised  and  observed  s-polarised. 


As  the  mean  light  level  is  subtracted  from  each  speckle  image  before  any  further 
processing,  this  ambient  light  is  not  expected  to  greatly  affect  the  correlation  values 
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to  be  obtained  (in  terms  of  invalidating  any  comparisons  to  theoretical  results). 


The  large  shift  to  higher  values  for  surface  #239,  in  the  cross  polarised  case,  may 
again  be  attributed  to  the  ambient  light  level.  The  shift  will  be  enhanced  by  the 
low  intensity  detected  compared  with  the  ambient  light  level.  When  considering  the 
integration  time  of  over  30  seconds  and  the  large  amount  ot  light  scattered  around 
the  laboratory  which  was  not  initially  depolarised,  it  is  surprising  that  the  shift  is 
not  greater. 


The  departure  from  the  general  negative  exponential  shape  of  the  histograms  in 
figure  3.12  is  not  so  easily  explained.  It  is  possible  that  it  is  due  to  a  large  amount, 
of  co-polariscd  light  reflecting  ofT  of  the  polariser  placed  before  the  camera,  and 
striking  the  surface  again,  producing  (albeit  with  less  magnitude)  another  cross 
polarised  speckle  pattern.  The  observed  histogram  would  therefore  be  a  combination 
of  the  two  patterns,  which  would  give,  as  seen  earlier,  a  more  Claussiau  shape.  The 
specular  geometry  of  the  experiments  would  support  this  hypothesis. 


However,  it  must  be  borne  in  mind  that  the  cross  polarised  results  for  surface  #239 
were  taken  near  the  practical  limits  of  the  equipment  used:  previous  work  with  this 
surface  (studying  the  average  intensity  in  a  given  direction)  reported  no  observed 
depolarisation!'11'.  The  experimental  rig  was  not  designed  with  a  view  to  studying 
very  low  amounts  of  depolarisation.  In  future,  the  rig  could  be  modified  to  cope 
with  the  above  problem,  allowing  thorough  scattering  studies  from  surfaces  such  as 
surface  #239  which  only  evidence  slight  multiply  scattering. 
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3.2.3  Processing  the  Speckle  Images 

After  having  recorded  and  stored  a  range  of  speckle  pattern  images  (the  range  could 
be  images  at  incident  angles  separated  by  one  degree  —  one  image  representing  one 
particular  alignment  of  incidence  and  observation),  the  data  must  be  processed  to 
obtain  values  for  the  correlations. 


Correlation 


The  numerical  method  employed  to  obtain  the  correlation  coefficient  uses  the  con¬ 
volution  theorem.  Tins  method  was  chosen  over  a  simple  summation  proredme  as 
it  gives  the  location  of  the  best  correlation,  without  presuming  it  is  an  exact  overlay 
of  the  two  images.  An  outline  of  the  computational  steps  involved  is  given  below: 


1.  Zero  the  mean  and  normalise  one  image  (the  reference  image). 

2.  Perform  a  fast  Fomier  transform  (FFT)  on  the  reference  image. 

3.  Find  the  complex  conjugate  of  the  transform. 

4.  Zero  the  mean  and  normalise  another  image  (a  data  image). 

5.  FFT  the  data  image. 

6.  Perform  a  pixel  by  pixel  multiplication  of  this  transform  with  the  conjugated 
reference  transform. 

7.  Find  the  inverse  FFT  of  the  above  product. 


8.  Search  for  the  peak  in  the  correlation  map  just  produced. 
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Normalisation  is  achieved  by  setting  the  sum  of  the  value  of  all  pixels  squared  to 
one  (£(a"2)  =  1). 

This  whole  procedure  results  in  just  one  number,  the  correlation  coefficient  between 
two  speckle  patterns  at  different  angles.  This  number  is  in  the  range  0  to  1,  1 
occurring  when  the  reference  is  correlated  with  itself  (identical)  and  0  occurring 
when  the  two  images  have  no  similarities  whatsoever.  To  obtain  the  curves  presented 
later,  a  range  of  ‘data'  images  must  be  correlated  with  a  single  ‘reference*  image;  i.e.. 
the  above  procedure  must  be  repeated  from  step  1  until  suflu  ient  values  at  different 
angles  have  been  procured.  One  set  of  speckle  pattern  images  is  enough  to  produce 
correlation  curves  for  any  reference  angle  within  that  set;  only  one  experimental  run 
is  necessary. 

Figure  3.13  shows  two  speckle  patterns  and  their  correlation  map.  Figure  3.14  is 
a  trace  through  the  centre  of  the  correlation  map,  the  height  of  the  peak  is  the 
correlation  coefficient  (note  that  the  width  of  the  peak  gives  the'  speckle  si/.e  in  the 
autocorrelation  case,  i.e.  when  the  speckle  patterns  are  identical).  The  peak  will 
not  appear  at  the  centre  of  the  graph  if  the  two  speckle  patterns  have  a  relative 
shift. 

The  small  fluctuations  near  the  x-axis  in  this  graph,  indicating  that,  small  corre¬ 
lations  exist  at  relatively  large  displacements  of  one  speckle  image  with  respect  to 
the  other,  are  an  artifact  of  the  finite  size  of  the  patterns  correlated.  If  the  corre¬ 
lated  fields  we  '  infinitely  large,  these  fluctuations  would  vanish  leaving  a  smooth 
line  along  the  axis.  Notice  that  the  height  of  the  fluctuations  reaches  a  maximum 
of  about  0.05;  this  is  the  level  of  uncertainty  expected  from  images,  such  as  these. 
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Figure  3.13:  Obtaining  the  correlation  coefficient  of  two  speckle  patterns.  The  cross- 
correlation  of  two  normalised,  zero  mean  speckle  patterns  giving  the  correlation  maj  . 


containing  400  speckles. 


This  minimum  level  of  confidence  will  allect  the  correlation  coelficient.  once  it  is 
reduced  to  that  extent.  It.  will  mean  that  in  later  correlation  curves,  the  minimum 
of  the  correlation  coelficient.  will  approach  this  level  and  not  drop  to  zero  as  may  be 
expected.  Various  methods  have  been  employed  to  extract  the  correlation  coelficient 


Figure  3.14:  A  slice  along  a  horizontal  hue  tlnough  the  centre  of  the  correlation 
map  in  figure  3.13.  The  peak  of  the  curve  is  the  correlation  coefficient. 
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from  the  coi  relation  maps  created  (see  below),  but  the  level  at  which  the  correlation 
is  no  longer  meaningful  is  the  same  for  all  methods,  the  differences  between  them 
being  purely  cosmetic  —  i.e.  some  give  more  information  below  this  level  than 
others,  but  information  in  this  region  is  unreliable. 


Mean  Intensity 

As  the  relative  intensities  of  the  speckle  patterns  have  been  disturbed,  it  is  not 
possible  to  derive  an  intensity  plot  by  simply  averaging  the  pixel  values  of  each 
image.  As  any  explanation  of  the  correlations  observed  must  also  apply  to  the 
intensities,  it  was  considered  necessary  to  perform  experiments  to  obtain  the  results 
under  similar  circumstances. 

Unfortunately,  previous  intensity  results  obtained  with  (essentially)  the  equipment 
described  below,  on  the  same  surfaces  as  studied  iu  this  work,  were  not  considered 
with  observation  in  the  specular  direct  ionbfif41!. 

The  illumination  of  the  surface  is  by  means  identical  to  that  used  to  obtain  the 
correlation  results.  However,  now  a  photomultiplier  tube  (PMT)  is  used  to  take  the 
results  instead  of  a  CCD  camera,  this  has  a  lens  placed  iu  front  of  it  to  collect  the 
scattered  light.  The  PMT  is  connected  to  a  high  tension  power  supply  and  its  output 
is  fed  into  a  phase  lock-in  amplifier.  The  amplifier  is  connected  to  the  ubiquitous 
PC  which  communicates  with  it  via  an  IEEE  interface.  A  beam  chopper  placed 
about  half  way  along  the  beam  is  connected  to  the  amplifier.  The  Sun  workstation 
is  no  longer  necessary  in  the  chain. 
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3.2.4  Enhancements 

Some  experimental  limitations  on  the  accuracy  of  the  correlations  can  be  removed 
during  processing,  these  will  be  discussed  now. 


Image  Alignment 

A  problem  noticed  while  sequentially  displaying  speckle  images  at  high  speed  - 
always  a  good  idea,  to  ensure  that  the  data  looks  valid  —  was  that  as  well  as 
‘evolving’  as  expected,  the  whole  pattern  was  slowly  oscillating  from  side  to  side; 
this  was  possibly  a  problem  with  the  stepper  motors.  The  oscillation  meant  that  any 
two  speckle  images  would  not  overlap  exactly.  For  the  purposes  of  the  correlation 
this  is  generally  not  a  worry;  however,  if  any  processing  is  to  be  performed  on  the 
speckle  images,  misaligned  patterns  could  pose  a  great  problem. 

The  speckle  patterns  are  aligned  by  performing  a  correlation,  but  instead  of  finding 
the  value  of  the  peak  in  the  correlation  map,  the  position  of  the  peak  within  the 
map  is  determined.  Tins  gives  the  oftset  which  must  be  applied  to  one  image  before 
it  overlays  the  other  —  the  whole  image  set  is  aligned  in  this  way  wit  h  the  reference 
image,  and  a  zero  border  is  set  around  each  image  which  absorbs  the  areas  which 
do  not  overlap. 

After  having  employed  this  technique  it.  was  found  that  correlation  values  increased 
slightly  for  the  same  set  of  data.  Obviously  some  lion-trivial  degree  of  correlation 
must  exist  before  this  technique  could  work. 
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Pattern  Shrinkage 

If  equation  (3.5  in  section  3.1.1)  is  expanded  to  depend  on  cartesian  coordinates, 
such  that  there  exists  a  tx  and  an  orthogonal  ry,  then  it  can  be  seen  that  any 
increase  in  the  beam  width  along  the  x  direction,  u>r,  would  decrease  rx  whilst 
leaving  ry  unchanged.  This  situation  occurs  for  increasing  angles  of  incidence  in 
the  experimental  work.  As  the  surface  is  tilted  to  higher  angles,  so  the  laser  beam 
spot  elongates  across  the  surface  along  the  horizontal  direction:  this  means  that 
the  speckle  pattern  will  shrink  horizontally  (and  also  change  to  some  extent)  at 
increasing  angles  of  incidence. 

As  the  shrinkage  is  not  accounted  for  in  any  theoretical  studies  later,  two  techniques 
were  employed  in  an  attempt  to  minimize  its  effect.  The  first  involved  placing  a 
mask  on  the  surface  allowing  only  a  certain  area  to  he  illuminated,  regardless  of 
the  angle  at  which  the  light  is  incident.  The  second  used  Mali,  thin’  images  in  the 
correlation  procedure  described  above. 

Figure  3.15  shows  a  comparison  of  experimental  correlations  taken  from  the  same 
surface  (#80),  with  and  without  a  mask  over  it.  As  expected,  the  correlations  at 
low  angles:  are  unaffected,  but  as  the  angle  increases  the  correlation  values  from  the 
surface  with  the  mask  over  it  are  greatly  increased.  The  improvement  is  lost  at  even 
higher  angles  as  the  mask  introduces  enors  of  its  own  as  its  finite  thickness  casts 
shadows  onto  the  surface. 

Since  the  mask  increases  the  correlations  in  some  circumstances  (hence  removes  a 
source  of  error)  and  never  reduces  them,  a  mask  will  be  placed  on  all  surfaces  to  he 
experimentally  studied. 
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Figure  3.15:  Graph  showing  the  improvement  in  correlation  values  obtained  at  a 
given  angle  by  masking  a  surface  to  reduce  pattern  shrinking.  The  solid  line  is  the 
result  fiom  the  surface  masked,  the  dashed  one  unmasked.  Note  that  the  vertical 
distance  between  points  shows  a  large  enhancement. 

For  the  other  method,  tall,  thin  images  were  obtained  by  truncating  them  in  the  ,r 
direction.  The  reasoning  behind  doing  this  is  that  the  shrinking  effect  will  displace 
speckles  further  away  from  the  centre  of  the  image  by  greater  amounts:  removing 
the  speckles  with  large  displacements  will  reduce  a  source  of  error  in  the  correlations. 
One  drawback  with  this  method  is  that  the  number  of  speckles  being  processed  is 
reduced,  and  hence  the  reliability  of  the  results  compromised. 


Improving  the  Reliability  of  the  Lower  Value  Correlations 

As  mentioned  above,  and  seen  in  figure  3.13,  there  exists  a  limit  below  which,  the 
‘true’  correlatioi  coefficient  cannot  be  easily  distinguished  from  small  correlation 
fluctuations  due  to  the  finite  size  of  the  speckle  images.  The  level  of  this  limit, 
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however,  can  be  reduced  simply  by  knowing  that  the  true  correlation  is  near  the 
centre.  So,  if  the  area  of  the  correlation  map,  over  which  a  maximum  is  searched 
for,  is  reduced  to  a  small  section  around  the  centre,  the  maximum  found  is  more 
likely  to  be  the  true  correlation  coefficient.  It  no  longer  needs  to  be  the  largest,  peak 
in  the  map. 

A  further  improvement  would  be  to  have  an  aligned  set  of  images,  by  a  method  se  ch 
as  that  described  above,  and  look  only  at  the  value  at  the  center  of  die  correlation 
map.  However,  in  order  to  align  the  images  some  degree  of  correlation  must  exist, 
so  it  is  not  possible  to  know  how  to  position  those  images,  which,  when  correlated, 
produce  a  value  below  the  limit.  One  possible  way  to  overcome  such  a  handicap  is  to 
keep  a  running  total  of  the  shift  needed,  and  correlate  each  image  with  the  previous 
one  in  its  set;  i.e.  the  shift  needed  to  be  imposed  on  the  third  image  to  align  it 
with  the  first  is  the  sum  of  the  shifts  needed  to  align  the  second  to  the  first  and  the 
third  to  the  second.  This  method,  however,  produced  shifts  which  are  inconsistent 
with  those  produced  by  a  direct  correlation.  The  most  likely  explanation  being  that 
by  Hacking  the  shift,  using  a  cumulative  process,  any  errors  will  he  cumulative,  so 
within  just  a  few  images,  the  shift  produced  will  be  wrong.  Since  the  shills  involved 
are  only  a  few  pixels  and  must  be  integer  numbers  of  pixels,  a  fairly  large  rounding 
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Chapter  4 

Single  Scatter  Analytical  Study 

4.1  The  Single  Scatter  Approximation 

la  this  section,  the  KirchhofT  approximation  is  applied  to  the  Ilelmholtz-KirchhofF 
integral  to  obtain  a  single  scatter  approximation  of  the  far  field  amplitude  of  the 
scattered  light.  An  analysis  of  the  range  of  surface  statistics  that  give  a  single 
scattering  surface  is  also  undertaken. 

4.1.1  Derivation 

The  Kirchholf  approximation  will  be  used  to  obtain  an  expressions  for  the  single 
scatter  component  of  the  scattered  amplitude  in  the  far  field  by  applying  it  to  the 
Ilelmholtz-KirchhofF  integral  (equation  (2.48)  from  chapter  two). 

Following  Beckmann^,  an  expression  for  the  single  scatter  component,  of  the  scat¬ 
tered  light  in  the  far  field  is  derived.  The  scattering  geometry  is  shown  (one  di¬ 
mensionally)  in  figure  4.1.  Light  is  incident  on  a  point  r  with  wave  vector  k,  and 
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scatters  with  wave  vector  k,;  the  positions  of  the  source  and  observer  with  respect 
to  the  scattering  points  are  given  by  R,  and  R.,  respectively. 


Figure  4.!:  Geometry  for  single  scatter  component  of  scattered  field. 


Given  an  incident  field  on  a  surface,  the  field  after  scattering  once,  can  be 
written  from  ecpiation  ('2.48)  as 


E,  =  —  /  /  tl3- - xl>~ —clS 

47r  J  Js  On  on 


where  xp  is  the  free  space  Green’s  function: 

gikj  ,RS 

4’  =  ~rT  (L2) 

The  field  at  the  surface,  and  its  normal  derivative  in  equation  (4.1)  are  given 
by  the  KirchhofF  approximation  (cf.  chapter  two), 

=  (l+K)Ei 
0E„  ,  „ ,  OE, 

-ss-  =  (l~K,*r 
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where  Tv  is  the  Fresnel  reflection  coefficient  at  the  point  of  scattering. 


Combining  equations  (4.1  and  4.3)  gives 


and  by  taking  the  incident  field  as  a  plane  wave  of  amplitude  Eo, 

Ei  =  E0e'k-V 


the  normal  derivative  at  the  surface  becomes 


oE,  lkil. 

-  =  iE0n.ktt  K‘ 

a  n 


Now  by  taking  the  single  scattered  component  of  light  out  to  the  far  field,  Ft,  can 
be  written  in  terms  of  the  absolute  position  of  the  observer  in  the  far  field,  R0, 


Rs  =  R0  -  r 


where  now  ko  and  k.,  are  identical.  From  equation  (4.1?), 


/e»t0Ro\  , 

4-  =  V-  r 


04’ 

—  =  -m.k,(-— )t  • 


Finally,  substituting  these  into  equation  (4.4),  it  becomes 

E,  =  ~  (£7^“)  jjs  zn‘  l^(k*  -  k‘)  -  (k<  +  M  e'(k'-k*)  lc/N 


(4.10) 


representing  the  single  scatter  contribution  to  the  scattered  electric  field,  expressed 
in  terms  of  the  quantities  in  figure  4.1. 
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Beckmann^10!  proceeds  to  develop  a  more  approximate  form  of  equation  (110)  for  a 
random  rough  surface  £(.t,i /)  with  sides  of  length  2 L  (L  >>  A)  which  removes  the 
explicit  reference  to  the  surface  normal.  To  distinguish  this  from  the  original  form 
(equation  (4.10),  the  electric  field  E 5)  will  be  represented  by  the  variable  A.  The 
single  scattered  amplitude  in  the  far  field  may  then  be  written,  ignoring  the  small 
“edge”  term,  as: 


A  = 


(4.11) 


where  V  =  k,  —  k3,  F3  is  a  geometrical  factor  and  Aq  is  the  field  which  would  be 
reflected  in  the  specular  direction  by  p-polarised  incident  light  if  the  surface  were  a 
perfectly  conducting  smooth  plane. 


Expanding  r  and  V  into  cartesian  components,  and  applying  the  geometry  shown 


4.1.2  Validity  of  the  Single  Scatter  Model 

A  similar  method  to  that  employed  to  estimate  the  KirchhofF  validity  will  he  used 
here  to  obtah  an  estimate  of  which  Gaussian  surfaces  are  smooth  enough  to  he 
considered  as  single  scatterers. 

An  intuitive  assumption  will  he  made,  that  double  scattering  first  conics  into  play 
(with  changing  angle)  when  the  path  of  the  light  between  two  scattering  points  is 
in  a  plane  parallel  to  the  mean  surface  plane.  This  condition  can  he  seen  in  figure 
4.3.  The  reasoning  behind  this  assumption  being  that  a  surface  which  is  not  rough 
enough  to  have  any  gradients  complying  with  this  condition,  will  not  only  scatt  r  the 
light  away  from  itself,  but  will  have  a  comparatively  large  distance  between  any  two 
sides  cf  a  ‘valley1.  This  will  make  it  almost  impossible  for  light  which  is  propagating 
with  a  gradient  greater  than  zero  to  strike  the  surface  again.  The  technique  does 
not  cater  tor  any  diffraction  effect  which  will  occur  as  the  distance  between  the  sides 
of  a  valley  approaches  a  wavelength. 
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Following  the  assumption,  the  surface  gradient  at  the  point  of  scattering  musi  be 

m  =  tan(jr/4  —  |0|/2)  (1.11) 

where  0  is  the  angle  of  incidence  of  the  light. 


Figure  4.3:  The  condition  required  before  double  scattering  is  evident. 


As  this  is  a  limit,  it  is  possible  to  say  that  the  surface  will  be  a  single  scattering  one 
if  the  vast  majority  of  its  surface  gradients  are  less  than  m  in  magnitude.  Since  the 
surfaces  studied  here  are  Gaussian,  the  surface  gradients  themselves  have  a  Gaussian 
distribution  with  zero  mean  and  standard  deviation 


(1.15) 


(cf.  chapter  two).  As  previously  determined,  99.7%  of  a  Gaussian  quantity  lies 
within  three  standard  deviations  of  the  mean.  The  condition  may,  therefore,  be 


written 

3\/2<r 


r 


<  tan(x/4  -  |0|/2) 


(bib) 
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Figure  4.4:  Surface  parameters  for  which  the  single  scatter  KirchholF  approximation 
(Beckmann  solution)  is  valid.  Boundaries  are  placed  at  the  three  standard  deviation 
threshold.  The  surfaces  are  plotted  with  the  wavelength  at  A  =  0 XCi/nn  and  for 
normal  incidence.  Surface  #80  (not  shown)  lies  well  within  the  region  in  which  the 
Beckmann  solution  applies. 

where  00.7%  ot  the  surface  gradients  are  within  the  limit. 

Figure  4.1  is  a  graph  of  equation  (4. 1C),  for  no'tnal  incidence,  with  the  Kirchhoff 
approximation  validity  graph  overlayed.  The  limit  of  the  joint  region  of  applica¬ 
bility  (the  white  area)  is  roughly  the  same  shape  as  and  within  the  one  percent 
unitarity  error  limit  found  for  the  region  of  validity  of  the  single  scatter  KirchhofT 
theory  by  numerical  simulation^  (for  all  angles  of  incidence).  The  close  agreement 
between  the  results  derived  from  the  statistical  propositions  about,  when  the  limits 
of  the  Kirchhoff  approximation  and  single  scattering  will  occur  and  the  numerical 
results  lends  a  great  deal  of  weight  to  the  simple  mechanism.  An  explanation  of  the 
double  scatter  results  shown  in  later  chapters  will  be  developed  on  the  basis  of  this 
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argument,. 

The  graph  also  shows  that  the  Kirchholf  approximation  is  more  widely  applicable 
than  the  single  scatter  approximation;  this  fact  will  be  exploited  in  chapter  live. 
Under  an  illuminating  wavelength  of  0-633 //m,  surface  #83  is  in  the  region  of  Kircli- 
lioff  applicability,  but  not  single  scattering,  whereas  surface  #23!)  lies  just  inside  the 
single  scatter  limit.  Hearing  in  mind  that  the  single  scatter  boundary  line  was  plot¬ 
ted  for  normal  incidence,  at  more  oblique  angles  a  surface  might  cross  the  boundary 
and  become  multiply  scattering  —  in  fact,  the  limit,  of  equation  (1.16)  is  achieved, 
for  surface  #239,  at  an  angle  of  incidence  of  (3.0  ±1.0  degrees  from  normal  (the  large 
uncertainty  due  to  the  errors  in  rr  and  r). 

4.2  Correlation  Coefficient  of  Two  Speckle  Pat¬ 
terns 


This  next  derivation  follows,  to  some  extent,  the  work  performed  by  beget  and 
rerriu^h,  the  stalling  point  being  the  single  scattered  field,  equation  (1.12). 


The  amplitude  correlation  coellicient  of  two  such  fields  Ai{0i,0-z,  0j)  and  ,  Oi-  Oa) 


(>MS) 


(1.17) 
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Ap,  ■/'3|  -'h);  /'Oj 

(2/7? 


IJ  f jL  'i-^Jjj)<..iryiy,-rww)^.(r.-Ic1  -y^),l.rvlr  XnW>H 


where  Vj,  i 'v  and  V7  are  as  in  equations  (l.ld)  and  wlieie  subscript  1  indicates 
dependenre  on  angles  0 j,  02  and  0 j,  one  of  tiie  required  arrangements,  and  subscript 
2  on  angles  <,•>,,  and  <pj,  the  other  arrangenieut  necessary. 


As  mentioned  in  chapter  two,  the  surfaces  to  be  examined  here  may  be  represented 
as  a  zero  mean  Gaussian  process  and  a  Gaussian  correlation  tuuctioii.  Using  the 
result  tor  any  zero  mean  Gaussian  process  ij. 

(( ••"')  =  e-r-Vd  ( i.H)) 


then 


_  ( -5<’il'  i,  +'  -21  l  'jO--Vi  !e)]  |.j  .)(jj 

where  cr2  =  (£2)  is  the  mean  square  surface  height  (the  variance  of  and 

(•»'i  -  •i'G‘  +  (// i  -  //.■)*' 


J-2.  Hi,  in)  =  exp 


(-1.21) 


is  the  Ganssian  autocorrelation  function  for  the  surface,  correlation  length  r. 


For  large  r, 


i  1/2)  ~  1  — 


(•>’i  -  -r-2)2  +  (1/1  -  I/O2 


(1.22) 


a  parabolic  approximation  which  is  more  generally  applicable  than  the  case  presented 
here  implying  that  this  theory  may  also  apply  to  surfaces  which  have  correlation 
functions  that  are  not  Gaussian. 


Putting  equations  (-1.20  and  4.22  into  equation  4. IS)  gives 


m 
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x  JJL  e,(Vvi yi  ~ v« w)g-  -  7* ~2  (yi  -in  ? dlJl  (hj2 

By  applying  the  variable  changes  tx  —  Xi  —  x2,  ty  =  yt  —  y2,  sx  =  x,  -f  x2  and 
■Sy  —  tji  +  7/2  then  cl.Vid.V2  becomes  \dsIdtx  and  clijidy2  —  \(lsyclt y,  so, 

(AiA2)  =  (4.24) 

2  4(2L)4 

rs/2 1  '  r  -•  --  -1 


X  /  j 

J-ClL 


\J  Cl  l  e:<l;>  - T^" — ^  '^2  )t^j  dtrj 


f2L  e*v,t-vn)%  ( rriL  exp  +  ?(V  +  v-  }k 

\J-ClL  {  I  r2  S  +  + 


r2  'v  1  2  JW/V  ^ 

Notice  that  if  Xi  and  x2  are  independent,  as  is  the  case  above,  they  will  form  an 
orthogonal  coordinate  system.  Within  this  system,  any  lines  of  constant  tr  will  be 
seen  to  be  orthogonal  to  lines  of  constant  sx,  implying  that  sx  and  tx  are  independent. 
The  same  argument  may  be  used  to  show  ty  and  sy  are  independent.  Note  that  the 
variable  transformation  used  at.  this  stage  in  roference[3 1]  was  not  orthogonal. 

For  large  «,  as  is  the  case  in  equation  (4.24). 

['  e~av2c'bvdv  m  r  e-av2eibva 
J-X  J-oo 

Using  the  result  for  a  >  0, 

r  -aA  ibv.  /*  M2 

J  e  e  civ  =  y  — exp  1 - 

gives 


(  4a 


(1.2  b) 


(1.26) 


f'*L  ftu, 

J-ClL 

r“‘  e«v»-v»*  exp 

J-X2L  y\a2VZiV^ 


<?2K,  K2 

7i 
■  u 


(V4,+VJ.2)2r2 


16  £t2v;,  k2 

{Vyi+Vy^r2 
16a2  V; 


*2  J 


X 


(4.27) 
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Now  substituting  the  simple  result  for  the  remaining  integrals 

rX 


J  eiaudu  =  2A"sinc(a.Y) 


(4.2S) 


into  equation  (4.27)  yields 

(.4 ,.4;>  =  ^(14,  -V„)L  (4.29) 

.,  ■  fv^„,  „  ,,1  t!  [  ((i4,  +  i-'„)2  +  (v;,  +  i;,)!)r2i 

X  s,„c  ^_(V„  -  V„)ij  exp  | - - j 

Normalising  as  in  equation  (4.17),  the  amplitude  correlation  coefficient  of  two  speckle 
patterns  is 


CA  =  e-^a(V*i-l'*2>a  sine 
x  exp 


v/2,,. 


sme 


s/2 

w*  -v*)/- 


yW,  -VU)£ 

II  f,12±!2.  ,  _  (^  +  ^)2  +  (\;,  +  vw)2V 

K^2V  H2  “  K*  VSIV^  )' 


(4.30) 


Equation  (4.30)  is  a  general  result  for  the  two-dimensional  amplitude  correlation 
coefficient;  however,  the  quantity  sought  here  is  an  intensity  correlation,  the  ampli¬ 
tude  correlation  is  not  relevant  to  this  work.  Assuming  that  the  speckle  amplitude 
is  a  complex  Gaussian  process,  tire  following  relation  holds'15!: 

(A  h&h)  =  (4.31) 

This  leads  to  a  simple  relationship  between  the  correlation  coefficient  of  the  intensity 
fluctuatoins  and  the  amplitude  correlation  coefficient^38! 

Cl  =  ICUI2  (4.32) 


Finally,  the  two-dimensional  intensity  correlation  coefficient  is 

C/(01,02,03,0I>2,03)  =  e~°2(Vli~VrW 


(4.33) 
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xsinc2  [f(VXl  -  VX2)L ]  sine2  [^(Vy,  -  Vyi)L\ 
x  exp  f-tL  + 

x  exP  [8^  V~  v2\  +  ~  v}2  v-,i-2  J 

4.2.1  Consequences  of  the  Intensity  Correlation 


It  can  be  seen  from  equation  (4.33),  that,  the  correlation  coefficient  will  have  its 
maximum  when  14,  —  142,  Vyi  ~  I'm  an<^  Hi  =  H2,  he.,  when  the  two  arrangements 
of  observation  and  illumination  are  identical.  It  may  also  be  seen  that  there  are 
various  factors  which  vary  with  only  one  of  14,  Vy  or  lr.,  and  that  sonic  of  these 
factors  may  be  removed  purely  by  geometric  considerations. 

The  most  obvious  limitation  to  impose  upon  equation  (4.33)  is  to  only  allow  'in 
plane’  scattering  in  both  arrangements,  i.e.,  set  O3  =  (p3  =  0,  which  from  equations 
(4-13)  means  that  Vy,  -  V'w  =  0.  The  intensity  correlation  becomes: 

Ci{0u02,<Pu<Pi)  =  exp  [-a2(\4,  -  V'-2)2]  sine2  —  (H,  -14J/' 


\r}(nK.nKi  (H,  -M4,)2 

x  exp  [w  [2n+2K~~^ . “ 


(4.31) 


It  now  makes  sense  to  force  the  same  step  for  the  14 ’s  and  set  14,  =  V42  —  0.  From 
equations  (4.13)  it  can  be  seen  that  this  is  easily  achieved,  namely  be  setting  0X  =  O2 
and  4>\  —  4>2-  This  corresponds  to  observation  always  in  the  specular  direction.  If 
(f)X  —  ip2  =  Qx  -+  60,  the  intensity  correlation  coefficient  between  two  speckle  patterns 


observed  in  the  specular  direction  with  illumination  angles  0X  and  0X  +  60  is 

Cj(0i,60)  =  exp  {cosOi-cosiOr+60))2  (4.35) 

Note  that  the  correlation  no  longer  depends  on  the  correlation  length  of  the  surface, 


T, 
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A  similar  result  might  he  expected  if  VXl  =  VXJ  =  c,  where  c  is  non-zero.  From  equa¬ 
tions  (4.13).  this  means  that  there  will  be  a  maximum  in  the  correlation,  (although 
it  is  not  possible  to  separate  out  the  Vz  dependence)  when 

sin  0X  —  sin  02  —  sin  <f> x  —  sin  02  (4.36) 

This  condition  has  been  reported  by  Michel^  for  rough  surfaces,  although  it  was 
presented  earlier  in  another  form  by  Leger  and  Perriid31)  (cf.  “memory  effect"  in 
chapter  one).  It  is  always  true  for  the  specular  case. 

Figures  4.5  to  4.7  examine  a  few  of  the  possibilities  equation  (1.36)  offers;  each  figure 
consists  of  a  correlation  graph  produced  with  equation  (-1.31),  and  an  associated 
graph  showing  the  angles  0X.  02,  0 t  and  02.  In  all  cases  0X  and  are  fixed  and 
0  1  is  the  controlled  parameter;  knowing  these  three  angles,  <p2  is  deteremined  with 
equation  (4.36). 

The  first  of  the  figures,  figure  4.5,  shows  the  specular  case,  02  =  0 x  and  <p2  =  <px,  the 
reference  being  at  0X  =  10  degrees;  the  shape  of  the  curve  is  that  given  by  equation 
(4.35). 

Figure  4.6  shows  the  case  where  0X  and  02  are  fixed  at  10  degrees  and  0  degrees 
respectively,  02  is  determined,  for  a  particular  <px,  by  equation  (4.36).  As  may  be 
expected  from  consideration  of  the  VZl  —  VZ2  term, 

(cos  -f  cos02)  —  (cos  0i  +  cos  02 )  (4.37) 


the  correlation  drops  off  at  a  slower  rate  than  in  the  specular  case  shown  in  figure 


Figure  4.5:  Top:  Cj(0i,  02<  <f>t,  (f>2)  plotted  for  the  specular  case  with  the  reference  at 
10  degrees  (02  =  G\  -■  10). 

Bottom:  4*2  is  determinied  by  equation  (4.3G),  in  this  case  <f>2  ~  <P i- 


The  last  figure,  figure  4.7,  shows  the  case  where  the  reference  image  is  in  back- 
scatter,  i.e.,  02  =  —0i\  <p 2  is  given,  for  a  particular  <p\,  b)  equation  (4.3G).  For 
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Q  (01,02,41!, <>2  ) 


Figure  4.7:  Top:  Cj(0i,  02,  0i,  <fa)  plotted  for  the  reference  m  back-$catter  (0 j  =  JO 
degrees  and  <?2  =  — 10  degrees). 

Bottom:  <f>2  is  determined  by  satisfying  equation  (4.36). 


is  identical  to  the  that  for  the  specular  case,  but  now  the  variation  in  VZ2 , 


cos  d> 1  +  cos  $2 


(4-39) 
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is  even  slower. 

These  results  are  extremely  counter  intuitive.  As  discussed  in  chapter  one,  the  high¬ 
est  degree  of  correlation  might  be  expected  in  the  specular  arrangement.  However, 
the  slower  decorrelations,  arising  because  of  the  slower  variations  in  the  Vg  terms, 
may  be  understood  by  considering  the  paths  the  light  travels  in  each  case.  For  the 
last  case  presented,  where  the  reference  is  in  the  back-scatter  direction,  the  angles 
between  the  various  paths  are  small  and,  since  it  is  the  angular  difference  which 
counts,  will  therefore  introduce  less  of  an  angular  decorrelation  than  in  the  other 
cases,  where  the  angles  between  the  paths  are  greater  (being  largest  in  the  specular 
case). 


In  plainer  terms,  the  directions  of  incidence  and  observation  are  all  very  similar  in 
the  back-scatter  case.  Changing  these  will  obviously  aifect  the  degree  of  similarity, 
but  not  to  a  major  extent.  Compare  that  to  where  observation  is  in  the  specular 
direction  —  the  directions  are  not  e\en  similar  in  the  first  place,  and  any  change 
makes  a  large  difference. 


4.2.2  Comparison  with  Experimental  Results 

Figure  4.8  is  a  plot  of  equation  (4.35)  with  0 i  =  10  degrees,  which  applies  to  correla¬ 
tions  taken  at  points,  compared  with  experimental  results,  which  necessarily  involve 
spatial  averaging  with  the  equipment  used.  It  shows  that  the  experimental  correla¬ 
tions  are  a  lot  lower  than  the  theory  predicts.  The  error  on  each  point  is  predicted 
to  be  around  ±5%,  being  the  statistical  uncertainty  on  a  sample  of  400  speckles, 
1/7(400). 
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Figure  4.8:  Equation  (4.35)  (solid  line)  plotted  with  experimental  correlations  (  +  ’s) 
taken  from  surface  #80  with  incident  light  at  10  degrees  and  s-polarised  and  observed 
speckle  patterns  in  the  specular  direction  —  also  s-polarised.  The  wavelength.  A  = 
0.633/m?.  The  error  on  each  point  is  ±5%. 

One  major  difference  between  the  theoretical  analysis  and  the  experimental  approach 
is  that  equation  (4.35)  produces  the  correlation  coefficient  at  a  point,  whereas  the 
experimental  situation  necessarily  deals  with  the  correlation  of  an  area  the  size  of 
the  camera’s  CCD  array. 

To  cope  with  this,  it  is  necessary  to  abandon  equation  (4.35)  and  turn  back  to 
equation  (4.33).  As  this  equation  provides  for  scattering  away  from  the  specular 
direction  and  out  of  plane  scattering,  it  is  possible  to  find  the  correlation  coefficient 
at  any  point  in  space,  not  just  at  the  point  in  the  specular  direction.  By  averaging 
equation  (4.33)  over  the  area  of  the  CCD  array,  a  result  similar  to  the  experimental 
observations  should  be  obtained. 

No  attempt  has  been  made  to  average  equation  (4.33)  analytically.  For  the  results 
to  be  shown  later,  the  equation  was  evaluated  for  each  pixel  in  the  CCD  array 
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contributing  to  the  portion  of  the  speckle  pattern  used  to  obtain  an  experimental 
result.  An  average  of  all  evaluations  gives  the  correlation  coefficient  to  be  compared 
to  experiment. 

Figure  4.9  shows  the  relationship  between  the  coordinates  of  a  point  on  the  CCD 
array,  (p,q),  and  the  angles  used  in  equation  (4.33).  This  relationship  is  explicitly 
stated  as  equations  (1.40). 


Figure  4.9:  Diagram  showing  the  relationship  between  a  point  on  the  CCD  array 
and  the  angles  ,0Z  and  03.  Note  that  the  centre  of  the  CCD  array  subtends  an 
angle  of  to  the  mean  surface  normal. 


where 


cos  0i 
sin  O3 


n 

9 

X 


(4.40) 


x 


d sin  0\  -f  pcos  0\ 
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z  =  dcosOl-psu\Ol 
R  =  yV  +  qi  A~2 


The  numerically  averaged  result  is  compared  against  experiment  in  figure  1.10.  This 
graph  is  identical  to  figure  4.8  in  all  other  respects.  It  can  be  seen  that  by  averaging 
equation  (4.33)  a  result  very  close  to  the  experimental  observation  is  obtained. 


Figure  4.10:  Equation  (4.33)  averaged  over  the  aiea  of  the  CCD  element  (solid 
line)  plotted  with  experimental  correlations  (+\s)  taken  from  surface  #80  with  the 
reference  at  10  degrees,  incidence  light  is  s- polarised  and  observed  speckle  patterns 
in  the  specular  direction  are  also  s-polarised.  A  =  0.633 /on. 

Having  now  determined  that  it  is  necessary  to  average  equation  (4.33)  to  gain  a 
similarity  with  experiment,  it  is  desirable  to  examine  the  effect  this  has  on  the 
various  arrangements  discussed  above  and  graphed  in  figures  4.5  to  4.7. 


Figures  4.11  to  4.13  are  the  equivalents  of  the  above  mentioned  graphs  but  now 
using  equation  (4.33)  averaged  over  an  area  around  the  observed  direction. 


As  expected  for  the  specular  case,  the  correlation  at  any  point  is  generally  less 
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Figure  4.12:  4>\,  V?)  (averaged)  plotted  with  the  reference  at 

grees,  and  0 2  =  0.  All  angles  satisfy  equation  (4.36). 


=  10  de- 


every  half  a  degree  —  if  a  finer  resolution  were  used  the  spike  might  be  expected  to 
have  an  angular  width  about  that  of  a  speckle. 


The  reason  for  the  reduction  in  the  correlations  is  that  the  terms  in  equation  (4.33) 
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Figure  4.13:  Cj[Oj,  02,  <Pi,  <p2)  (averaged)  plotted  for  the  reference  in  back-scatter 
(0X  —  10  degrees  and  02  =  —10  degrees).  All  angles  satisfy  equation  (1.36). 

which  varied  slowly  to  give  the  larger  correlations  in  the  previous  sect  ion,  in  partic¬ 
ular  VZJ,  vary  rapidly  with  02  while  0\  is  fixed,  as  is  the  case  when  averaging. 


It  is  interesting  to  note  that  figure  4.12  and  figure  4.13  are  now  in  line  with  the 
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intuitive  argument  presented  in  chapter  one,  justifying  the  use  of  the  specular  ar¬ 
rangement  when  considering  experimental  correlations. 


4.3  Results 


Figures  4.14  and  4.15  show  the  experimentally  observed  correlation  coefficients  for 
surface  #80,  plotted  with  the  resubs  of  averaging  equation  (4.33)  as  described  above, 
for  situations  where  incident  and  detected  light  are  both  s-polarised  and  p-polarised 
respectively.  Although  equation  (4.33)  does  not  take  any  polarisation  effects  into 
account,  it  is  obviously  adequate  for  predicting  the  observed  behaviour. 


Figure  4.14:  Experimental  intensity  correlations  for  surface  #80  with  incident  and 
detected  light  s-polarsed.  Curve  labels  represent  the  (incident)  angle  at  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  =  0.633/im.  The 
solid  curves  represent  the  results  of  averaging  equation  (4.33)  using  the  experimental 
parameters.  The  symbols  are  experimental  results. 

As  the  measurements  of  the  surface  parameters  obtained  had  some  margin  of  doubt 
associated  with  them,  figure  4.14  has  been  regraphcd  in  figure  4.16  with  the  aver- 
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Figure  4.15:  Experimental  intensity  correlations  for  surface  # SO  with  incident,  and 
detected  light  p-polarised.  Curve  labels  represent  the  (incident)  angle  at  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  —  0.633/m?.  The 
solid  curves  represent  the  results  of  averaging  equation  (4.33)  using  the  experimental 
parameters.  The  symbols  are  experimental  results. 

aged  equation  (4.33),  for  the  extremes  of  the  error  ranges  in  the  R.M.S  height  and 
correlation  length. 

For  the  sake  of  completeness,  the  average  intensity  plot  for  surface  # SO  is  given 
in  figure  4.17.  Both  the  incident  and  observed  light  are  s-polarised.  As  may  be 
expected,  as  the  angle  of  incidence  moves  further  away  from  being  normal  to  the 
surface,  the  amount  of  light  which  is  specularly  reflected  (as  opposed  t.o  any  other 
direction)  increases. 

As  a  matter  of  interest,  analytical  and  experimental  curves  are  ;uso  plotted  for 
surface  #239.  This  surface  is  predicted  to  show  double  scattering  more  evidently 
as  the  angle  of  incidence  increases.  For  the  co-polarised  situations,  figures  4. IS  and 
4.19,  theoretical  curves  as  above  are  plotted. 
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grey  area  represents  the  region  ol  uncertainty  in  the  measurement  of  the  suiface 
parameters. 

The  agreement  to  the  theoretical  curves  is  not  as  good  as  it  is  with  surface  #80.  This 


Figure  4.17:  Averaged  intensity  plot  for  surface  #S0  with  observation  in  the  specular 
direction.  The  illuminating  light  has  wavelength  0.633 /mi  and  it  and  the  detected 
light  are  s-polarised. 
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may  be  expected  as  surface  #239  is  bordering  on  the  margins  of  being  a  multiply 
scattering  surface.  It  must  be  stressed  that  due  to  large  gradients  involved,  the 
graphs  can  be  misleading:  to  see  the  discrepancy  between  the  experimental  and 
analytical  curves  in  its  true  light,  the  vertical  distance  between  an  experimental 
point  and  the  theoretical  curve  must  be  considered. 


Figure  4.18:  Experimental  intensity  correlations  foi  surface  #239  with  incident  and 
detected  light  s-polarised.  Curve  labels  represent  the  (incident)  angle  ai  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  =  0.63 ifj.ni.  The 
solid  curves  represent  the  results  of  averaging  equation  (4.33)  using  the  experimental 
parameters.  The  symbols  are  experimental  results. 


Now  that  some  double  scattering  is  possible,  cross-polarised  correlations  may  be 
obtained,  as  shown  in  figures  4.20  and  4.21. 

The  cross-polarised  component  is  very  much  less  intense  than  the  co-polarised  one,  so 
much  so  that  Kind-41!  reported  that  no  depolarisation  was  evident  in  scattering  from 
this  surface.  For  these  results,  long  exposure  times  were  required  to  produce  speckle 
images  with  pixel  values  in  the  range  specified  in  chapter  three.  The  exposure  time 
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Figure  4.19:  Experimental  intensity  correlations  for  surface  #239  with  incident  and 
detected  light  p-polarised.  Curve  labels  represent  the  (incident)  angle  at  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  =  0.633 ft  in.  The 
solid  curves  represent  the  results  of  averaging  equation  (4.33)  using  the  experimental 
parameters.  The  symbols  are  experimental  results. 

was  greater  than  30  seconds  for  most  of  the  data  shown,  coming  close  to  2  minutes 
for  the  higher  angles.  Although  the  large  exposure  times  allow  large  amounts  of 
ambient  light  to  enter  the  camera,  this  is  subtracted  when  the  speckle  images  are 
processed.  Further  problems  may  be  introduced  experimental  apparatus,  with  small 
amounts  of  detected  cross  polarised  light,  as  discussed  in  chapter  three. 

One  notable  feature  present  in  the  cross-polarised  curves  is  the  existence  of  extra 
humps  in  the  correlation  at  certain  angles  away  from  the  main  peak.  The  shapes  of 
these  curves  are  explained  in  chapter  six. 

The  behaviour  of  the  extra  correlation  humps  in  figure  4.20  is  interesting,  each 
individual  curve  has  peaks  at  around  three  degrees  either  side  of  the  main  peak,  but 
the  value  to  which  these  rise  varies,  becoming  lower  with  higher  reference  angles. 


Figure  4.20:  Experimental  intensity  correlations  for  surface  #239  with  incident  light 
s-po!arised  and  detected  light  p-polarised.  Curve  labels  represent  the  (incident) 
angle  at  which  the  reference  was  taken.  Observation  is  in  the  specular  direction  and 
A  =  0.633/t?n. 

Figure  4.21  shows  similar  behaviour,  but  with  the  extra  correlation  humps  being 
more  evident  at  the  larger  reference  angles.  At  the  lower  angles,  however,  the 
correlation  coefficients  seem  to  fluctuate  randomly  —  no  clearer  picture  is  given  if 
the  curves  are  examined  separately. 

Figure  4.22  and  figure  4.23  show  averaged  intensity  plots  for  surface  #239,  for  the 
case  of  observation  in  the  specular  direction,  for  s-polarised  light  incident  and  s-  and 
p-  polarised  light  detected  respectively.  The  relative  scaling  between  these  graphs 
has  not  been  preserved,  the  s-polarised  experiment  requiring  a  neutral  density  fdter 
in  the  system  to  attenuate  the  amount  of  light  incident  on  the  surface. 

As  can  be  seen  in  the  figures,  as  the  angle  of  incidence  moves  away  from  normal,  the 
amount  of  co-polarised  light  scattered  into  the  specular  direction  increases.  Con- 
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Figure  4.21:  Experimental  intensity  correlations  for  surface  #239  with  incident  light 
p-polarised  and  detected  light  s-polarised.  Curve  labels  represent  the  (incident) 
angle  at  which  the  reference  was  taken.  Observation  is  in  the  specular  direction  and 
A  =  0.633 fim. 

sequently,  the  amount  of  cross-polarised  light  (already  small)  decreases,  hence  the 


Figure  4.22:  Averaged  intensity  plot  for  surface  #239  with  observation  in  the  specu¬ 
lar  direction.  The  illuminating  light  has  wavelength  0.633/im  and  it  and  the  detected 
light  arc  s-polarised. 
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Figure  4.23:  Averaged  intensity  plot  for  surface  #239  with  observation  in  the  spec¬ 
ular  direction.  The  illuminating  light  has  wavelength  0.633/<?»  and  is  s-polarised. 
The  detected  light  is  p-polarised. 

long  exposure  times  needed  to  obtain  the  speckle  images. 


Since  there  are  no  obvious  side  lobes  in  figure  4.22  or  figure  4.23,  it  seems  that 
the  correlation  of  speckle  patterns  is  a  much  more  sensitive  probe  of  the  multiple 
scattering  regime  than  is  the  average  intensity.  This  final  piece  of  evidence  offered 
up  by  scattered  light  from  surface  #239  (that  the  surface  only  shows  slight  evidence 
of  double  scattering),  lends  a  great  deal  of  support  to  the  simple  model  near  the 
start  of  this  chapter  which  predicted  such  behaviour  on  the  basis  of  equation  (4.16). 
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Chapter  5 


Double  Scatter  Simulation 


5.1  Theory 


The  Kirchhoff  approximation  will  be  used  to  obtain  an  expression  for  the  double  scat¬ 
ter  component  of  the  scattered  light  in  the  far  field  by  applying  it  to  the  Helmholtz- 
KirchhofF  integral.  This  method  was  employed  by  Bruce^  for  one-dimensional  sur¬ 
faces  and  it  seemed  the  most  likely  candidate  for  an  extension  into  two  dimensions. 

Light  incident  on  the  surface  is  treated  as  being  scattered  into  two  states  —  one 
in  which  the  light  leaves  the  surface,  and  the  other  in  which  the  light  strikes  the 
surface  again  at  a  different  point.  Of  the  light  which  impinges  on  the  surface  again, 
some  will  scatter  and  leave  the  surface  and  some  will  strike  it  a  third  time.  Light 
continues  to  be  scattered  from  the  surface  until  it  has  all  escaped  (or  in  the  case  of 
a  non-perfect  conductor,  been  absorbed). 

Bruce  found  that  when  considering  perfectly  conducting  surfaces  for  which  the 
Kirchhoff  approximation  can  be  applied,  such  as  those  discussed  here,  over  97% 
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of  the  incident  light  had  left  the  surface  after  the  second  contact.!22!;  the  amount  of 
triple  (and  higher)  scattered  light  gave  a  negligible  contribution  to  the  final  scattered 
field. 

Only  the  single  and  double  scatter  components  of  the  scattered  field  will  be  dealt 
with  in  the  present  study.  The  small  power  loss  evidenced  with  this  approach  was 
deemed  an  acceptable  price  for  extending  the  calculation  to  deal  with  two  dimen¬ 
sional  surfaces. 


5.1.1  The  Single  Scatter  Contribution 

The  equation  to  be  used  for  the  single  scatter  contribution  is  equation  (1.10),  derived 
in  the  previous  chapter.  It  represents  the  contribution  to  the  electric  field  due  to  light 
which  only  interacts  with  the  surface  once.  Obviously,  for  surfaces  which  possess 
steep  sided  features,  light  will  scatter  from  more  than  one  point  before  it  finally 
escapes. 

5.1.2  The  Double  Scatter  Contribution 

Following  a  procedure  similar  to  that  used  to  gain  the  single  scatter  component,  an 
expression  will  be  derived  for  the  double  scatter  component  of  the  scat  tered  light,  in 
the  far  field. 

The  nomenclature  to  be  used  is  shown  in  figure  5.1,  the  subscripts  1  and  2  represent 
arbitrary  realisations  of  functions  at  points  1  and  2,  e.g.,  nj  is  the  surface  normal 
at  point  1,  and  is  a  function  of  the  coordinates  of  that  point. 


Figure  5.1:  Geometry  for  the  double  scatter  component  of  the  scattered  field. 
Given  an  incident  field  on  a  surface,  E\,  the  field  after  scattering  once  can  be  written. 


from  equation  (2. -IS),  as 


r  1  ff  n  tfV’ia  ,  dE i,  JO 

Ei  =  —  /  /  Eu~, - 4' i2— * — 

4  k  J  Jsi  a »i  a iii 


and  the  double  scattered  component  of  the  scattered  field  can  be  written  as 


r  _  1  f  f  TP  dt'va  dE-i,  ,c 

Ei  —  —  I  E-2 j-t - V’23— ; — ms2 

4?r  y  ctU2  <vn2 


where,  as  previously,  i/’  is  the  free  space  Green’s  function, 

-It, 

-  -nr 


The  fields  at  the  surface  and  their  normal  derivatives  in  equation  (5.1  and  equation 
5.2),  can  be  given  by  the  Kirchhoff  approximation,  such  that 


ElM  = 


(1  +  7vi)Ti 
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Els  —  (1  +  Tv  2  )  I' 2 

0EX,  _  OE, 

-  l'-Wl)aI,T 

^  -  o-^ 

U 112  ul\2 


(5,1) 


where  7£j  and  Tti  are  the  Fresnel  reflection  coefficients  of  the  surface  at  points  1  and 

2. 


Combining  equations  (5.1,  5.2  and  5.4)  loads  to  an  expression  for  the  double  scat¬ 
tered  field  depending  on  the  incident  field, 


(1  +K,)£\^2  -(1-R,K'i2~ 

C/11j  <711 1 


+  (1  -KiU'n  fd  +nl)Elr~i-  -(1  -  R.l)-:'U()L' 


ctoichij 


fhii  chii 


(I  Sul  Si  (5.5) 


Taking  the  incident  field  as  a  plane  wave  of  amplitude  Eo>  as  before 


Ei  =  E0CklVi 


t5.G) 


and  the  normal  derivative  at  the  surface  at  point  1  becomes 


Phi 

Uni 


=  7  Z:'o  11 ,  -k  1  f ,k|  r* 


(5.7) 


Also,  V'  12  may  be  written  from  equation  (5.3)  as 

Vll  = 


Ri 


(5.S) 


and  by  taking  the  double  scattered  component  of  light  out  to  the  far  field  and  so 
expressing  R3  as 

(r>-9) 


R.3  —  Ro  “  r2 
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then  (observing  that  ko  and  k3  are  identical), 


and 


d</’23 

On  2 


-  -m2.k3 


^ikoRo  \ 

~Ro  ) 


,~ik3.V2 
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The  remaining  terms  of  equation  (5.5)  are  derived  by  expanding  equation  (5.S)  into 
cartesian  coordinates  (i.e.,  Rj.  =  —  aq)2  -f  (1J2  —  J/i ) 2  +  (~2  —  ~i)2)  and  applying 

the  operator  (2.15).  Their  final  forms  are 

(C-iD)v  i2 

(-A  +  il3),!'l2  (5.12) 

-  ih 0)  +  'HAD  -  3 .-1C  +  /?/;]  cq, 

where 


eliii 

Rij'vi 
dn  2 
£012 


A 

ii 


112R2 

n2.k2 


C  = 


ni.Rj 

/?22 


£)  =  nj.ki 


(5.13) 


Thus,  for  a  given  arrangement  of  illumination  and  observation,  and  with  a  known 
surface  profile,  only  the  reduction  coefficients  remain  to  be  determined  in  equation 
(5.5).  With  a  perfectly  conducting  surface,  the  reflection  coefficients  can  be  set.  to 
+  1  or  —1  for  p-polarised  or  s-polarised  incident  light  respectively:  this  is  due  to  a 
perfect  conductor  being  a  perfect  reflector!39!.  For  the  one  dimensional  case,  where 
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no  change  in  polarisation  occurs  upon  scattering,  this  is  till  that  is  required;  but 
for  a  two  dimensional  surface  —  where  the  polarisation  can  change  upon  scattering 
out  of  the  incident  plane  —  changes  in  the  reflection  coefficients  must  be  taken  into 
account. 


5.1.3  Including  Polarisation 

Following  Mitzuer^  ami  Beckmann^,  the  electric  field  E  can  be  split  into  two 
orthogonal  components,  one  which  is  s-polarised  E~\  and  the  other  p-polarised 
E+.  The  scattering  geometry  is  shown  in  figure  5.2.  The  vectors  e~  and  e+  are 
unit  vectors  representing  the  s-polarised  and  p-polarised  electric  vector  alignments 
respectively.  x,y  and  z  arc  orthogonal  unit  vectors  describing  the  coordinate  system 
used.  The  scattered  field  can  be  expressed  in  terms  of  the  incident  field  as 


E~  1 

'IV  ri  l 

r+  r+ 

L  1  j-  1  ii  J 

.  E .+  . 

..here  E~  and  E*  is  the  incident  field,  E~  and  E+  the  scattered  field,  and  the  I's 
are  the  reflection  coefficients;  subscript  ||  indicates  the  component  of  the  reflection 
coefficient  which  does  not.  change  polarisation  am!  _L  indicates  the  component  of 
the  reflection  coefficient  which  maps  one  polarisation  state  to  the  other.  These 
reflection  coefficients  are  effectively  the  Fresnel  reflection  coefficients  (7v“  and  7v+) 
transformed  from  the  local  scattering  plane  into  the  world  frame.  Note  that  in  the 
world  frame  the  s-polarised  and  p-polarised  electric  vectors  lie  perpendicular  and 
parallel  to  the  plane  of  incidence  —  once  scattering  is  considered  from  a  local  point 
of  view  this  is  no  longer  true,  hence  the  coordinate  transforms.  These  new  reflection 
coefficients  are  derived  by  MitznerBT 

Fjj-  —  7L~  cos  cos  [3,  +  7l+  sin  /?,  sin  (3S 
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Figure  5.2:  Geometry  for  polarisation  showing  the  orientation  of  the  electric  field 
coordinate  systems  within  the  world  frame. 


rj|  =  7v"  sin  0,  sin  ft,  +  7v+'  cos  0t  cos  0a 

r:  =  -7Z~  sin 0,  cos  0a  +  7v+  cos  0t  sin  0a  (5.15) 

rt  =  —  R~  cos  /?,  sin  0a  +  7v+  sin  0,  cos  0a 

Pi  and  03  are  the  angles  between  the  local  scattering  plane  and  the  incident  and 
scattered  polarisation  coord  in  a1  ?  systems  respectively,  as  shown  in  figure  5.3,  and 
can  be  defined  as 

cos 0X -=  tx.e~  (5  16) 

or 


* 


si  nfix  =  q^.ej 


(5.17) 
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Figure  5.3:  Relative  orientation  of  the  electric  field  coordinate  systems  with  respect 
to  the  local  scattering  plane. 


for  x  —  i,s,  where 


n  x  kj 


(5.  IS) 


J-  i  |  |  \  1  ^  / 

|n  x  kx| 

is  a  unit  vector  along  the  intersection  of  the  local  scattering  plane  with  the  wave 
fronts,  as  shown  (along  with  qx)  in  figure  5.4,  and  e~  is  the  direction  of  the  s- 
polarised  component  of  the  electric  field.  Expressed  in  world  coordinates. 


The  process  of  keeping  track  of  the  polarisation  can  be  thought  o’  as  first  a  trans¬ 
formation  of  the  (e~,ef)  system  into  one  which  lies  in  the  (t,,  q,)  plane.  This  is 
followed  by  scattering  the  light  according  to  equation  (5,1)  or  (5.2).  Finally  the  co¬ 
ordinate  system  (which  now  lies  in  the  (ts,q,)  plane)  is  transformed  to  the  (ej\e+) 
arrangement.  The  overall  effect  is  a  change  in  polarisation. 

A  limitation  must  now  be  imposed:  the  Fresnel  reflection  coefficients  only  apply  to 
a  specular  reflection  from  a  flat  plane,  and  so  the  above  procedure  will  not  work  if 
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Figure  5.4:  Detail  of  the  relationship  between  the  electric  field  coordinate  system 
and  the  local  scattering  plane. 


k,,  k,  and  n  do  not  lie  within  the  same  plane.  When  they  do,  t,  =  =  t,  and  the 

polarisation  is  correct.  Having  imposed  this  limitation,  the  polarisation  change  no 
longer  depends  explicitly  on  the  surface  normal,  and  t  may  be  rewritten  in  terms  of 
the  incident  and  scattered  wave  vectors: 


k,  x  k, 
|k,  x  k,| 


(5.21) 


A  model  is  now  needed  to  cope  with  light  scattered,  out  of  the  local  plane  of  inci¬ 
dence  (a  diffraction  effect  supported  by  the  Kirchhoff  method  used).  If  any  light  is 
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scattered  in  some  direction  ks,  from  a  scattering  element,  with  incident  light  in  the 
direction  k,,  then  the  existence  of  a  virtual  plane,  oriented  such  that  the  scattered 
light  is  specularly  reflected,  may  be  considered.  The  polarisation  of  the  light  is  then 
given,  as  above,  by  applying  the  reflection  coefficients  on  this  plane  with  t.,  and  ts  as 
in  equation  (5.21).  The  actual  orientation  of  the  virtual  plane  is  never  needed  as  k, 
and  k,  are  already  determined.  This  approximation  to  whatever  polarisation  change 
is  occurring  is  justified  since  the  surfaces  to  be  studied  here  must  conform  to  t.lu* 
Kirchhoff  approximation,  i.e.,  most  of  the  light  scattered  from  any  point  will  be  as 
if  from  a  plane  mirror.  A  qualitative  argument  successfully  explaining  polarisation 
changes  in  enhanced  back-scatter  has  been  based  on  a  similar  assumption^11!. 

Clearly  it  is  necessary,  from  equation  (5.1 1),  to  work  out  both  the  s-  and  p-polarisod 
components  of  £'2  (cf.  equation  (5.1)),  before  a  single  realisation  of  £3  (equation 
(5.2)),  be  it  s-  or  p-polarised,  can  be  determined. 


5.2  Implementation 


Given  that  the  single  and  double  scatter  components  of  the  electric  field  are  now 
realised,  the  total  scattered  field  is  simply  the  sum  of  the  two  components. 

The  angular  intensity  correlation,  the  quantity  sought  here,  depends  on  two  realisa¬ 
tions  of  the  intensity,  if  =  /(ki.ka)  and  /2  =  /(kjjkg),  taken  Irom  N  realisations 
of  the  surface,  and  is  obtained  using  equation  (2.15)  from  chapter  two: 


7!  5  h 

N 


(5.22) 
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The  total  electric  field  can  be  expressed  as  ail  intensity  (defined  as  the  energy  crossing 
a  unit  area  in  unit  time)  by: 

I  =  EE ‘  (5.23) 

4?r 

Examining  equation  (5.22),  it  can  be  seen  that  any  constant  factors  in  /,  and  hence 
(from  equation  (5.23))  in  the  electric  field  E,  may  be  discarded  when  only  considering 
the  correlation.  Note  that  this  only  applies  to  factors  common  to  both  the  single 
and  double  scatter  components  of  E  and  I  —  their  relative  scaling  must  not  be 
disturbed.  This  removal  ot  all  terms  constant  with  respect  to  angle  means,  especially, 
that  the  correlation  does  not  depend  upon  /?0,  a  term  which  (as  it  is  the  position  of 
observation  in  the  far  field)  tended  towards  infinity. 

5,2.1  Discretisation  for  Computer  Modeling 

The  integration  over  the  surface  in  equations  (4.4)  and  (5.5)  shall  be  performed 
simply  by  multiplying  the  value  calculated  at  the  centre  of  a  unit  by  the  area  of  that 
unit,  as  shown  in  figure  5.5.  This  makes  discretisation  of  the  equations  simple;  it  is 
achieved  merely  by  changing  the  integral  signs  to  sums  over  the  ,r  and  y  directions 
and  leplacing  dS  with  A5„v,  the  area  of  the  unit.  Obviously  this  is  a  very  simplistic 
approximation  to  the  actual  integral  —  a  more  accurate  calculation  might  extend 
Simpson’s  rule  into  two-climensions. 

It  is  generally  deemed  necessary,  while  considering  discretisation  of  equation  (5.5) 
and  all  associated  parts,  that  the  phase  terms  present,  namely  e-,kl  r> ,  and 

must  be  viewed  with  respect  to  their  variation.  In  order  to  maintain  a  realistic 
idea  of  this  phase  over  a  discretised  surface,  its  variation  must  be  much  smaller  than 
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Figure  5.5:  Numerical  integration  performed  by  taking  the  mid-point  ot  two  reali¬ 
sations  of  a  function  and  multiplying  it  by  their  separation. 

27T  between  two  consecutive  surface  points.  Taking  just  the  k0R2  term,  as  this  varies 
fastest,  the  condition  can  be  written  as 

A(k0R2)  «  2k  (5.24) 

where  A  indicates  the  change  over  consecutive  points.  Replacing  k0  by  2ir/\.  the 
condition  becomes 

A R2  «  A  (5.25) 

The  largest  change  in  R2  will  be  of  the  order  of  the  distance  between  points  on 
the  discretised  surface,  so  to  track  the  phase  reliably  across  a  discretised  surface 
there  must  be  many  sample  points  per  wavelength.  There  is,  however,  a  problem 
with  this  approach:  the  discretised  surface  will  have  so  many  points  on  it  that 
the  calculation  would  be  unrealistically  slow;  also,  the  fast  shading  algorithm  (see 
below),  would  require  so  much  memory  that  it  would  not  be  possible  to  implement 
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it.  The  integration  will  therefore  be  performed  with  a  grid  separation  controlled  by 
the  correlation  length  of  the  gradient  of  the  surface,  this  being  the  c"1  half  width  of 
the  correlation  function  of  the  gradient  (from  equation  (2.30)): 

CS'S‘[t)  =  (5.26) 

When  this  is  equal  to  c-1,  the  value  of  t  =  Ts',  the  correlation  length  of  the  gradients 
(found  by  Newton-Raphson  iteration),  is  approximately  given  by 

t$i  ss  C.51104r  (5.27) 

i.e.,  for  a  surface  with  a  correlation  length  r  =  4.31  A,  surface  points  will  be  separated 
by  4.31  A  =+=  0.51 104  =  2.20A.  If  2048  points  are  available,  then  the  surface  side  should 
be  y/204S  x  (2.20 A)*  =  99. 74 A. 

Surface  segments  spaced  by  this  distance  will  still  have  similar  gradients,  and  hence 
each  segment  may  be  considered  to  be  locally  flat.  The  distances  between  points  on 
a  surface,  will  now  constitute  a  random  distribution  (cf.  Monte  Carlo  method).  It 
can  be  seen  from  figuie  5.6  that  the  phase  spread  over  an  ensemble  average  of  100 
surfaces  is  relatively  uniform,  implying  that  it  does  not  need  to  be  tracked  on  a  finer 
scale.  The  randomness  of  the  surface  points  is  enhanced  by  adding  a  small  random 
offset  to  each  x  and  y  coordinate  (which  is  linear  and  up  to  a  maximum  of  A/ 10) 
and  interpolating  to  find  the  correct  s  coordinate  on  the  surface;  i.e.,  the  surface 
realisations  are  not  confined  to  a  fixed  grid  which  may  impose  its  own  structure  onto 
any  results. 

To  perform  the  integration,  a  sum  over  flat  surface  segments  must  be  conducted  (a 
requirement  of  the  Kirchhoff  approximation).  Somehow,  these  flat  facets  must  be 
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Figure  5.6:  The  distribution  of  tin.-  phase  differences  between  two  points  binned  50 
times  over  the  wavelength,  on  a  100  x  100  wavelength  surface  consisting  of  20-SS 
points,  with  parameters  a  =  1.42A  and  r  =  4. 31  A,  over  an  ensemble  average  of  100 
surfaces.  Point  to  point  shading  has  been  taken  into  account. 


obtained  from  the  data  available,  namely  a  square  grid  of  surface  heights  (given  the 
nature  of  the  surface  generation  technique  which  relies  on  fast  Fourier  transforms  — 
see  below).  In  actual  fact,  the  grid  must  be  disturbed  to  create  a  random  distribution 
of  points,  but  for  the  purposes  of  this  explanation  a  square  grid  will  be  considered. 

Since  it  is  unlikely  that  four  independent  points  will  all  lie  in  a  flat  plane,  each 
‘square1  of  the  surface  (defined  by  the  four  corner  points),  must  be  split  into  two 
triangles,  both  of  which  are  uniquely  defined  planes,  as  shown  in  figure  5.7.  This 
complicates  the  summation  procedure  slightly  as  it  is  now  necessary  to  cater  for  the 
two  sub-triangles  within  each  unit. 
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The  A SXy  required  for  a  triangular  unit  can  be  defined  as 


a c 

A5^  "  ~2n7~ 


(5.2S) 


the  area  of  a  surface  facet  with  a  surface  normal  n  (n.  is  the  component  which  lies 
along  the  z  direction)  on  a  triangular  base  of  area  AxAy/2.  The  positions  of  the 
mid-point  of  each  facet  (see  figure  5.8)  are  given  by  averaging  the  positions  of  the 
three  corners,  for  facet  A  the  mid-point  is 


At 
Ay 
A , 


3 

Piv+Plv  +Piv 

3 


=  Pi,  +  3  A-r 

=  Pi.  + 


Pj.+Pj.+P^ 

3 
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Figure  5.8:  Geometry  for  calculating  the  positions  of  surface  facets. 


and  for  facet  B 


5 1  +fjj  ~t  P*  i 
3 

f\ ,,  +>Ji,j  +!\ ,, 


3 


-  Pu  +  §A.r 
=  P\v  +  ^Ay 


(5.30) 


where  the  points  Pj,  P2,  P3  and  P4  are  four  realisations  at  the  corners  of  the  square 
unit-  The  surface  normal  for  each  facet  can  be  derived  by  taking  the  vector  cross 
product  of  two  sides  of  a  facet: 

nA  =  P2Pi  x  P2P3  =  -At /(P3t  -  P2l)x  +  A x{Pu  -  P2i) y  +  AxAyz  (5.31) 
nB  =  P4P3  x  P,P,  =  At/(Pi_.  -  P4,  )x  -  A.x(P3l  -  Piz  )y  +  AxAiyz  (5.32) 

Of  course,  these  vectors  still  need  to  be  normalised  by  dividing  each  component  by 
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As  mentioned  above,  the  randomness  of  the  distribution  of  points  on  the  surface 
needs  to  be  increased,  this  is  achieved  by  shifting  the  corners  of  each  square  by  a 
small  random  amount.  The  surface  normal  may  still  be  defined  by  a  cross  product 
of  two  vectors  describing  the  sides, 

n.-i  =  a.4  x  b.4  =  P2P\  x  P2P3 

nB  =  as  x  bB  =  P4P3  x  PAP{  (5.33) 

although  the  solution  is  not  so  trivial  as  before,  it  is  nevertheless  a  simple  calculation. 


The  most  significant  complication  introduced  by  randomising  the  surface  nodes  is 
that  the  surface  area  (IS  of  each  facet  needs  explicit  calculation  (see  figure  5.9),  not 
particularly  difficult  in  itself,  but  meaning  one  more  quantity  to  keep  track  of.  The 
area  of  a  triangle  described  by  two  vectors  joined  at  their  bases  is. 


dS  =■  - be  =  ^6|a  —  ob| 


where 


o  = 


a.b 

ab 


(5.3-1) 


(5.35) 


and  a  and  b  are  two  vectors  describing  the  sides  of  the  facet  as  in  equation  (5.33), 
assigned  such  that  6  >  a.  Note  that  c  may  be  calculated  via  the  dot  product  of  a 
and  the  normal  to  b,  but  this  method  requires  a  normalisation  and  hence  a  square 
root.  Bearing  in  mind  that  calculation  speed  is  of  the  essence,  the  above  method 
was  chosen  as  it  avoids  the  slow  square  root  and  uses  a  number  of  faster  operations. 


By  this  point  in  the  development  of  the  theory,  it  is  not  practical  to  write  the 
calculations  as  a  single  equation  or  even  just  a  few  equations.  However,  this  is  of  no 
concern  when  considering  that  it  will  be  implemented  on  a  computer  —  in  fact  it  is 
an  advantage  to  have  expressed  the  theory  as  many  small  parts. 
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Figure  5.9:  Geometry  for  calculation  of  area  between  two  vectors. 

5.2.2  Calculating  Shaded  Areas 


So  far,  no  allowance  has  been  made  for  the  effect  the  surface  height  variation  has  in 
terms  of  casting  shadows;  all  points  on  the  surface  have  been  assumed  illuminated 
at  all  times.  Clearly,  for  steep  sided  features  of  the  surface,  especially  at  high  (away 
from  normal)  angles  of  incidence,  there  are  going  to  be  portions  of  the  surface  which 
are  in  shadow.  The  inverse  ot  this  also  applies:  points  on  the  surface  which  cannot 
be  seen  from  the  observation  direction  cannot  directly  contribute  to  the  scattered 
light  in  that  direction. 

In  the  double  scatter  case,  the  effect  of  point  to  point  shading  must  also  be  taken 
into  account.  Light  cannot,  reach  a  point  on  one  side  of  a  hill  from  a  point  on  the 
other  side  —  this  must  be  accounted  for.  Note  that  if  the  surface  is  discretised 
into  N  points,  then  the  point  to  point  calculation  will  have  to  be  performed  from  1 
point  to  N  points  N  times,  i.e.,  N 2  point  to  point  calculations.  It  is  vital  that  this 
calculation  is  computationally  as  fast  as  possible.  A  direct  ray  trace,  possibly  the 
most  obvious  method  for  determining  shaded  regions,  is  too  slow  to  be  practiced. 
Fortunately,  it  is  not  the  only  solution  available. 
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Considering,  for  now,  a  one  dimensional  surface,  figure  5.10  shows  t  he  shaded  regions 
of  that  surface  when  light  is  incident  on  it.  The  vector  s  is  defined  as  that  which  lies 


light. 

in  the  direction  of  the  line  joining  the  last  illuminated  point  with  the  current  point. 
Starting  at.  the  left  side  of  the  discretised  surface,  as  each  point  is  stepped  through, 
s  is  redefined.  If  a  point  lies  in  an  illuminated  area,  then  s  will  simply  he  the  vector 
joining  two  consecutive  points  (cf,  sj.  If,  on  the  other  hand,  the  point  lies  in  a 
shaded  region,  then  s  will  be  the  vector  joining  the  start,  of  that  shaded  region  with 
the  current  point  (cf.  S2).  To  determine  if  the  current  point  is  illuminated  or  shaded, 
the  dot  product  of  s  and  the  normal  to  the  wave  vector  of  the  illuminating  light,  m 
(as  shown  in  figure  5.10),  is  taken.  If 


s.m  >  0 


(5.36) 


then  the  current  point  is  illuminated  and  the  tail  of  s  is  redefined  as  this  point, 
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othenvise  the  point  is  shaded  and  the  tail  of  s  is  left  alone. 

For  illumination  to  the  light  of  the  surface,  the  direction  the  surface  is  processed  in 
must  be  reversed,  i.e.,  stepping  starts  at  tn  ght  hand  side  and  the  definition  of  m 
is  such  that  it  still  points  ‘up’.  To  calculate  observed  points,  ‘illumination’  can.  be 
considered  to  be  from  the  observation  direction  —  ‘shaded’  regions  are  then  areas 
which  are  net  observable,  and  ‘illuminated’  regions,  areas  which  are.  Note  that  this 
algorithm  only  touches  upon  each  surface  point  once. 

As  long  as  the  observed  direction  is  in  the  plane  of  incidence,  a  tsvo-dimensioual 
surface  can  be  thought  of  (for  the  purposes  of  the  above  shadow  calculations)  as 
many  one-dimensional  surfaces  lying  side  by  side. 

A  similar  method  to  that  above  can  be  used  for  point  to  point  shadow  calculations 
(sec  figure  5.11),  the  only  difference  being  that  the  wave  vector,  and  hence  it's 
normal  m,  arc  no  longer  constant.  However,  the  point  to  point  calculation  cannot 
be  extended  to  that  for  a  two-dimensional  surface  so  trivially. 

The  solution  to  the  problem  of  calculating  point  to  point  shadows  is  to  work  with  a 
cylindrical  system  —  cartesian  vectors  are  no  longer  the  optimum  tool.  Fortunately, 
the  condition  s.m  >  0  can  be  rewritten  in  term;  of  gradients.  For  this  case;  as  you 
step  out  along  the  surface,  if  the  gradient  of  the  line  joining  the  current  point  to 
the  starting  point  is  greater  than  the  maximum  gradient  so  far,  then  the  point  is 
illuminated.  The  gradient  concept  can  easily  be  extended  into  the  two-dimensional 
case,  as  no  speo.ucation  of  direction  has  been  made. 


The  final  problem  of  determining  which  order  to  process  the  points  on  the  surface 
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Figure  5.11:  Shaded  regions  of  a  one-dimensional  surface  due  to  point  to  point 
scattering. 


is  still  not  simple  —  the  most  computationally  effective  being  to  spiral  out  from 
the  ‘source’  point  and  calculating,  only  once  for  each  point,  whether  that  point  is 
shaded  or  not.  The  difficulty  lies  in  working  out  which  points  lie  along  the  line  of 
sight  between  the  current  point  and  the  source  position.  Mere,  the  most  sensible 
course  of  action  is  to  precalculate  the  point  to  point  dependencies;  i.e.,  for  each 
point,  create  a  list  of  points  which  will  be  directly  affected  by  a  change  in  it.  Figure 
5.12  shows  the  concept  schematically.  These  lists  can  be  created  in  such  a  way  as  to 
form  a  template  which  can  be  overlayed  on  the  surface  at  any  position.  The  template 
must  be  twice  the  area  of  the  surface;  thus  ensuring  that  the  centre  of  the  template 
can  be  positioned  at  one  corner  of  the  surface  while  the  whole  surface  still  remains 
under  its  influence,  even  though  only  one  quarter  of  it  is  used.  If  the  template  were 
the  same  size  as  the  surface,  then  only  when  the  template  is  positioned  above  the 
surface  centre  would  the  whole  surface  be  processed. 


All  that  need  be  done  now  is  to  spiral  out,  at  each  point  working  out  the  gradient 
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Figure  5.12:  Creation  of  a  template  for  determining  point  to  point  shadows  on  a 
two-dimensional  surface.  The  dependencies  are  created  by  tracing  a  straight  line 
between  a  point  and  the  centre,  inserting  a  reference  to  that  point  into  the  list  of  all 
points  the  line  passes  through.  Point,  47  has  been  added  to  the  list  of  every  point 
that  has  an  effect  on  it;  a  completed  list  for  point  47  is  also  shown, 

of  the  slope  from  the  centre  to  the  current  point.  This  is  compared  with  the  value 
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of  the  ‘maximum1  gradient  associated  with  that  point;  if  it  is  the  new  maximum 
the  point  is  marked  as  illuminated,  otherwise  as  shaded.  The  calculated  gradient 
must  still  replace  the  maximum  gradient  of  each  point  in  the  dependency  list  if  it  is 
greater,  regardless  of  whether  the  original  point  was  illuminated  or  not. 

This  last  statement  may  seem  unreasonable:  how  can  a  point  which  is  not  illumi¬ 
nated  cast  a  shadow?  The  problem  here  is  that  each  ‘point'  actually  has  a  finite 
size  —  the  whole  defined  as  shaded  or  illuminated  by  the  state  of  the  centre  of  that 
point.  It  is  possible  though,  that  one  corner  of  the  point  might  be  illuminated,  the 
centre  being  shaded;  this  could  cast  a  shadow,  which  if  ignored,  would  propagate  an 
error  though  to  all  dependencies.  The  cjuestion  of  whether  the  whole  point  should 
be  set  as  shaded  or  illuminated  is  unresolved;  a  better  model  could  give  each  point 
a  more  analogue  value  depending  on  the  percentage  of  its  area  shaded.  However,  it 
should  be  noted  that  as  the  surface  varies  smoothly,  only  those  point  s  on  a  boundary 
between  light  and  dark  will  not  be  totally  one  state  or  the  other  (a  small  percentage 
of  the  total  number  of  points).  This  implies  that  for  a  point  which  is  on  the  bound¬ 
ary,  over  half  of  its  area  will  be  of  the  same  state  as  its  centre.  Finally,  if  the  surface 
has  been  discretiscd  to  a  sufficient  degree,  adjacent  points  will  not  have  a  noticeably 
different  effect  on  the  light  which  scatters  from  them,  therefore  some  averaging  of 
the  shadcd/illuminat.ed  state  will  occur. 

The  speed  of  this  process  can  be  improved  by  watching  for  when  the  gradient  and 
distance  from  the  initial  position  take  the  light  above  the  surface  —  at.  this  point, 
all  dependents  in  the  point’s  list  can  be  marked  as  shaded. 


Figure  5.13  shows  a  two-dimensional  surface  profile  (left,  image)  and  the  clfect  of  the 
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Figure  5.13:  Representations  of  a  computationally  generated  surface  profile  showing 
the  results  of  the  shading  techniques.  The  surface  height  at  a  point  is  given  by  tin- 
brightness  (light  is  high).  The  ldack  regions  indicate  shaded  areas. 


two  shading  techniques,  one  of  illumination  from  a  certain  angle  (middle  image), 
and  the  other  of  point  to  point  shading  (right,  image).  The  dark  line  down  the  right 
hand  side  of  the  middle  image  is  present  because  the  images  actually  represent  point 
heights,  but  the  shading  calculations  work  on  the  connections  between  them;  i.e.  a 
surface  with  a  side  of  32  points  will  have  31  facets  along  it. 

Extending  the  two-dimensional  shadow  calculations  from  a  rectangular  grid  into 
one  in  which  each  rectangle  is  subdivided  into  two  triangles  is  fairly  simple.  For 
the  illumination  and  observation  calculations,  each  ‘row’  of  the  surface  can  still 
be  considered,  except  now  a  row  has  twice  as  many  points.  The  point  to  point 
calculation  requires  slightly  more  thought,  two  templates  are  required,  one  for  each 
of  the  two  possible  ‘source’  triangles.  More  care  must  also  be  taken  when  spiraling 
out  —  the  same  process  outlhied  above  can  be  used,  but  the  order  in  which  the  sub- 
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triangles  are  processed  must  be  intelligently  defined:  the  one  closer  to  the  centre 
first. 

The  application  of  the  methods  described  above  is  not  limited  to  calculating  shaded 
regions  for  numeiical  calculations.  The  author  believes  that  these  methods  are  the 
fastest  possible  for  determining  lines  of  sight,  with  potential  use  in  various  areas. 


5.2.3  Numerical  Random  Surfaces 

A  numerical  random  surface  with  a  standard  height  deviation  of  a  and  correlation 
length  r  can  be  created  simply  by  convolving  Gaussian  white  noise,  of  standard 
deviation  a.  with  a  Gaussian  correlation  function,  of  1/e  half  width  t.  This  is 
shown  pictorially  in  figure  5.14.  Computationally,  the  convolution  is  performed 
by  a  multiplication  in  Fourier  space.  Zero  mean  Gaussian  white  noise  of  standard 
deviation  cr  has  (from  Parseval’s  theorem)  a  Fourier  transform  which  is  also  zero 
mean  Gaussian  white  nuise  of  standard  deviation  a.  The  Fourier  transform  of  the 
correlation  function  is  obtained  by  applying  a  fast  Fourier  transform  algorithm  to 
an  array  containing  a  Gaussian  envelope  function  with  a  1/e  half  width  of  r.  This 
transform  can  be  remembered  and  need  only  be  calculated  once  for  each  execution 
of  the  program,  no  matter  how  many  surfaces  may  need  to  be  generated. 

The  actual  process  of  creating  a  random  surface  consists  of  obtaining  a  random 
number  from  a  Gaussianly  distributed  setl4'l,  multiplying  this  number  by  the  rele¬ 
vant  entry  in  the  correlation  transform  array  and  storing  it  in  a  result  array.  After 
every  element  of  the  result  array  has  been  allocated,  an  inverse  fast  Fourier  trans¬ 
form  is  performed  upon  it.  The  final  result  is  an  array  of  Gaussianly  distributed 
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Figure  5.14:  Convolution  of  Gaussian  white  noise  with  a  Gaussian  correlation  func¬ 
tion.  The  method  is  equally  valid  for  two-dimensional  surface  creation. 


m 


numbers  with  a  standard  deviation  a  and  a  Gaussian  autocorrelation  function  with 
correlation  length  r. 


5.3  Results 

5.3.1  Single  Scatter  Component  of  Simulation 

Since  the  single  scatter  component  and  the  analytical  equation  derived  in  the  previ¬ 
ous  chapter  both  have  a  common  origin,  some  degree  of  similarity  between  the  two 
should  be  expected.  Figure  5.15  shows  the  analytical  and  computational  correla¬ 
tions  obtained  for  the  specular  geometry  with  the  reference  at  eight  degrees.  The 
simulation  was  run  for  a  surface  with  parameters  identical  to  those  of  surface  #80, 
consisting  of  2048  points  and  with  a  side  length  of  730A.  A  square  region  with  24 
points  per  side  was  illuminated  (shading  was  incorporated)  and  the  ensemble  aver- 


# 
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age  consisted  of  600  surfaces.  The  analytical  curve  was  plotted  using  the  expression 
for  specular  observation  and  with  no  averaging,  equation  (4.35)  in  chapter  four. 


Figure  5.15:  Analytical  and  single  scatter  computational  intensity  correlations  for 
observation  in  the  specular  direction  with  the  reference  at  8  degrees,  a  ~  2.27 ;.nn, 
t  —  20.7/<777.  and  A  =  0.633/4777.  Both  curves  have  data  points  every  half  degree. 


The  differences  in  the  two  curves  may  be  attributed  to  two  factors;  the  approxima¬ 
tions  made  in  the  analytical  derivation  and  the  non-infinite,  discretised  surface  in 
the  simulation.  The  relatively  small  number  of  sampling  points  per  surface  might 
be  expected  to  reduce  the  computational  correlations,  but  as  a  reduction  is  not 
evident,  the  arguments  presented  previously  in  favour  of  large  distances  between 
surface  points  appear  to  be  justified.  However,  although  this  appears  unlikely,  there 
may  still  be  a  netting  off  of  two  opposite  effects,  one  increasing  the  correlation  and 
one  decreasing  it,  which  could  mask  an  error  introduced  by  any  particular  approxi¬ 
mation.  Shadowing  is  not  likely  to  have  much  effect  for  such  a  low  sloped  surface. 
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Figure  5.16  shows  the  effect  of  including  polarisation  tracking  in  the  simulation,  the 
curves  show  the  correlations  of  the  polarisation  components.  The  s-polarised  curve 
coincides  with  the  the  calculation  run  with  no  polarisation  tracking.  In  all  cases  the 
incident  light  was  s-polarised,  all  simulation  details  are  as  above. 


Figure  5.16:  Single  scatter  intensity  correlations  with  polarisation  for  observation 
in  the  specular  direction  with  the  reference  at  8  degrees.  The  intensity  of  the  p- 
polarised  component  was  about  1  O'4  times  smaller  than  the  s-polarised  one.  a  — 
2.27/tm,  t  =  20.7/rm  and  A  =  0.633/tm. 


An  interesting  point  to  note  here  is  that  the  simulation  predicts  a  small  amount  of 
depolarisation  in  the  specular  direction.  It  is  generally  accepted  that  single  scatter 
calculations  do  not  predict  depolarisation  in  this  direction,  and  any  cross  polarised 
component  noticed  in  experimental  situations  has  been  attributed  to  some  degree 
of  multiple  scattering^. 


The  polarisation  change  shown  here  is  an  artifact  of  equation  (4.10)  and  the  method 
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chosen  to  control  polarisation,  the  reflection  coefficients.  The  cross-polarising  reflec¬ 
tion  coefficient  becomes  zero  for  any  specular  scattering  and  on  the  face  of  things 
this  will  mean  no  depolarisation.  However,  when  considering  that  a  small  amount 
of  light  may  be  scattered  into  the  specular  direction  from  slightly  tilted  surface 
facets,  the  last  term  of  equation  (4.10),  not  multiplied  by  the  reflection  coefficient, 
is  non-zero,  and  makes  a  contribution  to  the  cross-polarised  light. 

This  result  requires  a  more  in  depth  study  of  polarisation  changes  from  scattering  ele¬ 
ments.  However,  since  the  predicted  depolarisation  is  very  small  (the  cross-polarised 
component  is  about  four  orders  of  magnitude  smaller  than  the  co-polarised  one),  it 
will  not  be  noticed  by  the  experimental  equipment,  and  may  be  ignored  for  the  pur¬ 
poses  of  this  work.  It  will  make  an  insignificant  contribution  to  any  cross-polarised 
component  due  to  multiple  scattering.  It  may  be  that  the  result  is  merely  an  error, 
due  to  the  specular  approximation  used  to  calculate  the  reflection  coefficients  or 
the  finite  size  of  the  surface  facets.  Alternatively,  it  may  be  that  the  technique  to 
control  polarisation  using  reflection  coefficients  can  not  be  applied  to  light  which  is 
not  locally  specularly  reflected,  or  that  the  presumption  that  no  depolarisation  is 
present  in  the  specular  component  of  light  scattered  from  a  single  scattering  rough 
surface  is  invalid.  Future  work  might  like  to  consider  this. 

5.3.2  Experimental  Results  from  a  Multiply  Scattering  Sur¬ 
face 

Figures  5.17,  5.18,  5.20  and  5.21  show,  for  incident-observed  polarisation  slates  of 
s-s,  p-p,  s-p  and  p-s  respectively,  experimental  correlation  results  for  surface  #83. 
Within  each  figure,  curves  are  plotted  for  reference  angles  of  2,  10,  20,  30,  40  and 
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50  degrees. 


Figure  5.17:  Experimental  intensity  correlations  for  surface  #83  with  incident  and 
detected  light  s-polarised.  Curve  labels  represent  the  (incident)  angle  at  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  =  0.633 /mi. 


Figure  5.18:  Experimental  intensity  correlations  for  surface  #S3  with  incident  and 
detected  light  p-polarised.  Curve  labels  represent  the  (incident)  angle  at  which  the 
reference  was  taken.  Observation  is  in  the  specular  direction  and  A  =■  0.633/<rn. 
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As  can  be  seen  in  figure  5.17  and  figure  5. IS  —  data  from  co-polariscd  arrangements 
—  the  narrowing  of  correlations  at  higher  angles,  or  rather  the  widening  at  lower 
angles,  observed  in  the  single  scattering  case  (cf.  chapter  four)  is  less  evident.  As 
expected  then,  this  surface  does  not  conform  to  the  analytical  theory  in  chapter 
four.  Figure  5.19  emphasises  the  failure  by  showing  the  results  of  the  analytical 
expression  compared  to  the  results  in  figure  5.17. 
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Figure  5.19:  Experimental  data  as  in  figure  5.17  —  surface  #$'■!,  incident  and  de¬ 
tected  light  s-polarised  —  with  averaged  analytical  curves  (dashed  lines)  as  described 
in  chapter  four  overlayed. 

Figure  5.20  and  5.21  are  more  interesting;  these  are  plots  in  the  cross-polarised  case, 
showing  the  double  scatter  component  of  the  scattered  field.  They  are  similar  to 
the  co-polarised  graphs,  figure  5.17  and  figure  5.18,  in  as  much  as  the  widths  of 
the  peaks  are  the  same  and  uniform,  but  at  the  lower  angles  an  extra  bump  in 
the  correlation  coefficients  has  appeared  at  the  sides.  When  considering  that  the 
experimental  equipment  correlates  two  separate  speckle  patterns  (not  just  the  point 
intensities,  as  do  the  calculations),  this  suprising  result  means  the  speckle  produced 
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by  one  angular  arrangement,  in  the  double  scattering  case,  is  remarkably  similar  to 
that  produced  by  another  totally  different  arrangement.  Chapter  six  discusses  why 
this  could  be  so. 


Figure  5.20:  Experimental  intensity  correlations  for  surface  #83  with  incident  light 
s-polarised  and  detected  light,  p-polariscd.  Curve  labels  represent  the  (incident) 
angle  at  which  the  reference'  was  taken.  Observation  is  in  the  specular  direction  and 
A  =  0.633/(7?i. 

To  show  the  extent  of  the  double  scatter  extra  correlation  peak  more  clearly,  a  surface 
plot  of  the  correlation  coefficients  is  given,  in  figure  5.22,  over  all  combinations  of  the 
two  correlated  angular  arrangements,  for  surface  #83,  in  the  case  where  s-polarised 
light  is  incident  and  p-polarised  light  detected.  Although  a  peak  is  evident  either 
side  of  the  main  01  =  02  ridge,  this  can  be  considered  as  simply  a  repetition  of  data 
—  everything  on  one  side  of  the  ridge  should  be  visible  on  the  other.  The  peak  of  the 
extra  correlation  occurs  when  the  reference  image  is  at  two  degrees  and  correlated 
with  one  at  eight  and  a  half,  or  vice  versa  due  to  the  symmetry.  The  data  was  taken 
with  an  eight  bit  Pulnix  CCD  camera  which  was  used  before  the  low  noise  Princeton 
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Figure  5.21:  Experimental  intensity  correlations  for  surface  #83  with  incident  light 
p-polarised  and  detected  light  s-polarised.  Curve  labels  represent-  the  (incident) 
angle  at  which  the  reference  was  taken.  Observation  is  in  the  specular  direction  and 
A  =  0  633/mi  . 

Instruments  camera  was  available. 

To  help  compare  these  results  to  the  discussion  of  intensity  results  previously  re¬ 
ported  with  this  surface^11!,  figure  5.23  shows  the  intensity  plot,  for  surface  #83, 
where  as  usual,  the  detector  is  in  the  specular  direction  and  the  illuminating  light 
has  a  wavelength  of  0.633 fim.  The  figure  shows  two  graphs  for  the  average  inten¬ 
sity,  one  with  s-polarised  light  detected,  the  other  p-polarised.  Both  curves  were 
obtained  with  s-polarised  light  incident. 

As  with  previous  average  intensity  results  presented,  in  the  co-polarised  situation, 
the  light  level  detected  tends  to  increase  as  the  angle  of  incidence  moves  away  from 
normal.  The  cross  polarised  case,  L  vvever,  shows  some  interesting  behaviour;  the 
intensity  drops  and  then  rises  back  up  again  before  tailing  off  as  with  surface  #239. 


m 
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Figure  5.22:  Surface  plot  of  correlation  coefficients  for  surface  #83  with  inc  lent 
light  s-polarised  and  observed  speckle  patterns  p-polarised. 

The  rise  in  intensity  is  one  of  the  ‘side  lobes’  noticed  by  Mendez  and  O’Donnolb"! 
when  scattering  from  this  surface,  occurring  here  at  a  larger  angle  than  they  observed 
due  to  the  geometrical  arrangement  used.  No  back-scatter  peak  is  seen  at  normal 
incidence  because  the  equipment  was  not  set  up  to  observe  in  that  direction  —  the 
detector’s  view  of  the  surface  is  obstructed  by  a  mirror  used  in  the  beam  delivery. 

Figure  5.24  shows  a  close-up  of  the  p-polarised  component  of  figure  5.23.  Two  peaks 
are  clearly  evident  in  the  intensity  profile.  The  first  peak  is  the  side-lobe  previously 
mentioned.  The  second  peak,  not  observed  in  the  work  referenced  abo/e  —  most 
likely  because  of  the  geometrical  arrangement  used,  is  discussed  in  chapter  six.  It  is 
shown  to  be  a  second  order  effect  (where  the  back-scatter  peak  is  zeroth  order  and 
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Figure  5.23:  Average  intensity  plot  tor  surface  #33  with  observation  in  the  specular 
direction  and  with  s-polarised  (squares)  and  p-polarised  light  (circles)  detected.  The 
illuminating  wavelength,  is  0.633/<m  and  s-polarised. 


the  side-lobe  first  order). 


Figure  5.24 :  Average  intensity  plot  for  surface  #83.  This  graph  shows  an  expanded 
view  of  the  p-polarised  component  of  figure  5.23 
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5.3.3  Computational  Results 

Figure  5.25  shows  the  correlation  curve  produced  by  the  computational  method 
detailed  above  for  a  surface  with  rms  height  and  correlation  length  as  for  surface 
#83.  The  graph  shows  the  single  and  double  scatter  components  for  s-polarisod  and 
p-polarised  light  observed  (with  the  incident  field  s-polarised),  and  also  the  combi¬ 
nations  of  the  single  and  double  components,  the  total  curves,  for  each  polarisation 
state. 

As  discussed  in  chapter  four,  there  is  no  point  trying  to  directly  compare  this  tesult 
with  an  experimental  one  as  the  latter  is  spatially  averaged,  reducing  the  correlation 
values.  A  superficial  examination  reveals  that  the  total  computational  values  are, 
as  expected,  higher  than  the  experimental  ones  at  any  given  angle.  It  may  be 
seen,  however,  that  the  distinctive  bumps  present  in  the  experimental  curves  are 
noticeably  lacking  in  figure  5.25. 

Figure  5.26  is  a  curve  produced  by  the  method  when  an  ensemble  average  is  taken 
over  only  100  surface  realisations.  As  can  be  seen,  extra  ecu  relation  peaks  arc- 
now  present,  at  the  lower  angles.  The  reason  they  are  not  observed  in  figure  5.P5 
is  because  they  have  been  averaged  out.  This  computational  method  gives  extra 
correlation  peaks  of  different  heights  and  at  different  positions  for  each  run  of  100 
frames.  Averaging  over  greater  numbers  of  frames  washes  the  peaks  out. 

The  inconsistent  manner  in  which  the  extra  peaks  appear  may  be  explained  by 
the  following  points.  A  possible  mechanism  to  describe  the  origin  of  the  extra 
correlation  peaks,  given  in  chapter  six,  relies  heavily  on  the  phase  change  of  the 
light  scattered  between  sides  of  a  well  in  the  surface.  The  main  approximations 
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Figure  5.25:  Computational  correlation  curves,  for  a  surface  with  parameters  identi¬ 
cal  to  $urface#S3,  with  the  reference  image  at  10  degrees  and  showing  the  different 
components.  The  incident  field  has  a  wavelength  of  0.633/' m  and  is  s-polarised. 
The  ensemble  average  was  taken  over  S00  frames  of  a  201S  point  surface  with  a  side- 
length  of  GOA. 


necessarily  introduced  into  the  computational  method  in  order  to  get  some  results, 
namely  the  large  spacing  between  surface  points,  will  mean  that  this  very  critical 
distance  (around  seven  wavelengths  for  surface  #S3)  is  constrained  by  the  spacing 
and  size  of  the  surface  facets.  There  are  less  than  fi  e  surface  facets  over  this  critical 
distance.  The  random  additions  to  the  coordinates  of  the  surface  points  to  give 
a  uniform  phase  distribution  of  the  light  scattered  from  all  points  of  the  surface 
will  not  help  in  this  localised  situation  —  the  critical  distance,  over  any  five  facets 
contributing  to  the  extra  correlation  peak,  being  randomly  adjusted. 


Nevertheless,  although  it  is  not  possible  to  compare  figure  5.25  and  figuie  5. 26 
directly  to  the  experimental  results,  they  do  show  clearly  which  components  of  the 
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Figure  5.26:  Computational  correlation  curves,  the  variable  parameters  are  identical 
to  that  of  figure  5.25  except  that  the  ensemble  average  was  taken  over  only  100 
frames. 


correlation  curves  are  a  result  of  the  single  scatter  and  double  scatter  processes. 
It  is  clearly  evident  that  the  double  scatter  process  is  responsible  for  introducing 
side  peaks  into  the  correlation  curves.  It  may  also  be  seen  that  the  cross  polarised, 
double  scatter  component  of  the  simulation  and  the  total  curve  are  the  same;  proving 
that  the  cross-polarised,  single  scatter  component  (discussed  above)  has  a  negligible 
effect  on  the  overall  result. 

Figure  5.27  and  figure  5.28  show  graphs  of  the  correlation  coefficients  for,  respec¬ 
tively,  s-polarised  and  p-polarised  components  of  the  total  (i.e.,  single  plus  double) 
field  scattered  from  the  computational  statistical  equivalent  of  surface  #83  using 
the  method  above.  All  the  results  were  taken  oven  a  2048  point  surface  with  a 
side  iength  of  60  wavelengths  and  a  wavelength  of  Q. 633/m?.  Each  graph  shows  four 
curves,  with  the  reference  angles  set  at  10,  20,  30  and  40  degrees. 
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Figure  5.27:  Total  (single  +  double)  computational  correlation  curves  for  surface 
#83  for  various  reference  images  (at  10,  20,  30  and  40  degrees).  The  incident  field 
has  a  wavelength  of  0.633/mt  and  is  s-polarised.  The  correlations  consider  only  the 
s-polarised  scattered  component.  The  ensemble  average  was  taken  over  600  frames 
of  a  2048  point,  surface  with  a  side  length  of  60A. 

It  is  immediately  noticeable  that  both  of  these  figures  show  roughly  the  expected 
shapes  for  the  correlation  curves.  The  curves  in  each  diagram  are  all  about  the  same 
width  irrespective  of  angle.  The  finding  differs  from  the  single  scatter  prediction  that 
the  widths  of  the  peaks  should  reduce  as  the  angle  of  incidence  increases,  but  is  in 
accordance  with  the  experimental  results  (cf.  figure  5.17  and  figure  5.20). 

The  above  results,  derived  from  the  computational  model,  go  to  show  that  the 
departure  of  the  single  scattering  Beckmann  theory  from  the  experimental  results 
arises  soley  because  of  its  failure  to  cope  with  multiple  scattering.  In  fact,  it  has 
been  shown  that  consideration  of  just  the  s:ngle  and  double  scatter  components  of 
the  scattered  field  is  sufficient  to  bring  theoretical  calculations  into  line  with  the 
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Figure  5.28:  Total  (single  +  double)  computational  correlation  curves  for  surface 
#83  for  various  reference  images  (at  10,  20,  30  and  40  degrees).  The  incident  field 
has  a  wavelength  of  0.633/im  and  is  s-polarised.  The  correlations  consider  only  the 
p-polarised  scattered  component.  The  ensemble  average  was  taken  over  600  frames 
of  a  2048  point  surface  with  a  side  length  of  GOA. 

experimental  results  (from  surface  #83). 


As  shown  in  chapter  four,  surface  #83  was  expected  to  fall  within  the  regime  of 
the  Kirchhoff  approximation,  but  not  to  conform  with  the  Beckmann  solution.  The 
similarity  between  experimental  results  and  results  from  the  computational  model 
lacing  the  same  constraints  shows  that  the  conceptual  basis  of  the  statistical  analysis 
of  the  surface  parameters  is  trustworthy.  An  arguement  is  presented  in  the  next 
chapter  which  is  based  upon  similar  concepts  and  can  explain  the  shapes  of  various 
scattering  results  from  doubly  scattering  surfaces. 
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Chapter  6 


Discussion  and  Conclusion 


6.1  Single  Scatter 


It  has  been  shown  that  for  a  Gaussian  random  rough  surface  with  its  nils  surface 
height  and  correlation  length  within  certain  ranges,  such  that  the  surface  may  be 
considered  to  be  single  scattering,  Beckmann  theory  adequately  explains  all  exper¬ 
imental  results  presented  once  allowances  for  the  particular  experimental  situation 
are  incorporated.  As  an  existing,  well  established  theory  copes  well  with  the  situa¬ 
tion,  no  further  explanation  is  deemed  necessary  here. 


If  any  future  work  is  proposed  in  this  area,  it  should  be  to  quantitatively  examine 
the  effect  that  spatially  averaging  the  speckle  around  the  observed  direction  has  on 
the  correlation  coefficients. 
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6.2  Double  Scatter 

For  Gaussian  random  rough  surfaces  with  rms  heights  and  correlation  lengths  such 
that  they  may  be  considered  as  double  scattering,  but  still  conforming  to  the  re¬ 
strictions  on  the  Kirchhoff  approximation,  experimental  results  were  found  to  depart 
from  the  Beckmann  theory  (which  only  considers  a  single  scattering  mechanism)  to 
a  greater  degree  as  the  surface  .  ameters  fall  further  from  the  ranges  that  indicate 
single  scattering.  Polarisation  changes  were  also  noticed  which  are  not  addressed  by 
the  Beckmann  approach. 

As  there  is  no  analytical  theory  known  to  the  author  which  can  account  for  multiple 
scattering  from  two  dimensional  surfaces,  a  computational  model  has  had  to  provide 
the  theoretical  comparison  for  the  results  obtained.  Vast  demands  on  computer 
resources  limited  the  extent  of  this  comparison.  However,  a  very  simple  model  is 
now  proposed  which  explains  some  features  of  the  experimental  results,  namely  the 
peaks  in  the  correlation  curves  and  their  relation  to  the  side  lobe  in  the  intensity 
plots. 

The  first  part  of  this  analysis  will  concentrate  on  an  examination  of  the  scattered 
intensity  envelope. 

6.2.1  Scattered  Intensity  Envelope  Considerations 

The  following  is  based  on  an  intuitive  argument  that  some  observable  macro-effect 
(e.g.,  a  side  lobe)  is  simply  related  to  micro-causes  (e.g.,  phase  differences  due  to 
a  collection  of  individual  scattering  paths).  This  is  not  a  new  approach;  various 
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authors  (mentioned  below)  have  explained  the  double  scatter  effect  in  this  way,  but 
only  in  a  qualitative  manner. 

Consider  figure  6.1,  which  shows  one  valley  on  a  surface  and  one  possible  pair  of 
interfering  paths  ace  and  deb ,  both  of  which  travel  along  the  same  line  between  the 
two  points  of  contact  with  the  surface.  All  incident  light  is  at  an  angle  0  with  respect 
to  the  mean  surface  normal,  and  only  the  light  scattered  in  the  specular  direction  is 
considered.  As  vector  quantities  are  used,  the  representation  can  be  applied  to  two 
dimensional  surfaces. 


Figure  6.1:  A  pair  of  the  many  possible  interfering  paths  within  a  valley  of  a  rough 
surface. 


Previous  arguments  along  these  linesMM  have  concluded  by  stating  that  the  back- 
scatter  peak  observed  when  illuminating  multiply  scattering  surfaces  is  due  to  con¬ 
structive  interference  between  two  equal  and  opposite  paths,  one  the  time  reversal 
of  the  other  (a  less  general  situation  than  with  paths  ace.  and  deb).  These  paths 
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arc  guaranteed  to  constructively  interfere,  all  other  paths  resulting  in  no  net.  inter¬ 
ference.  This  argument  may  be  extended  to  include  the  so  called  ‘side  lobes’  by 
simply  considering  the  first  order  interference  fringe  between  more  general  paths!3'1'. 
Figure  6.2  shows  the  equivalent  of  the  side  lobes  evident  when  observation  is  in  the 
specular  direction. 
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Figure  6.2:  A  graph  showing  the  cross-polarised  intensity  scattering  envelope  from 
surface  #83.  This  is  a  duplicate  of  figure  5.24  in  chapter  five.  Observation  is  in  the 
specular  direction  and  the  incident  light  has  a  wavelength  of  0.633//m. 


Examining  the  arrangement  more  closely,  the  phase  difference  between  the  paths, 
<f>,  is  easily  determined: 

<t>  =  (k,.r!  -  k2.r2) -- (ki.r2 -k2.ri) 

=  (ki  +  k2).(r!  -  r2)  (6.1) 


For  the  specular  case  considered. 


k,  = 


2tt 

T 


(sin  Cx  —  cos  Oz) 


(6.2) 


m 


e 


6.2.  DOUBLE  SCATTER 


175 


k2 


-(sin$x  +  cos  0  z) 


giving  the  phase  difference 

47T 

(/)  =  -—  sin  0Ax  (6.3) 

A 

where  Ax  is  the  distance  between  the  two  scattering  points  along  the  x  axis.  As  no 
concessions  to  a  one  dimensional  surface  have  been  made,  this  is  true  even  for  a  two 
dimensional  surface  —  it  is  a  consequence  of  the  in-plane  scattering  geometry  used. 


For  constru.  live  interference,  <p  must  be  an  integer  number  of  2tr,  i.e.. 


2nx  =  sin  0Ax 

A 


(6.4) 


where  n  is  an  integer,  or  equivalently, 


sin  0  — 


nX 

2A^ 


(6.5) 


The  question  is  now  posed  as  to  what  A.r  is  in  a  macro  sense.  I3v  averaging  equation 
(6.5),  over  an  ensemble  of  wells,  the  surface  may  be  considered  as  a  whole,  such  that. 

sinfl  = -  (6.6) 

2  »/a* 

where  the  quantity  i)±d.  has  simply  replaced  A.r,  and  may  be  considered  as  the  mean 
well  width  (across  the  direction  of  the  incident  plane). 

Equation  (6.6)  can  now  qualitatively  explain  the  features  of  figure  6.2.  The  back- 
scatter  peak  is  the  zeroth  order  fringe  n  —  0,  i.e.,  when  0  =  0]  and  subsequent  peaks 
are  the  higher  orders,  n  =  1,2...,  their  actual  position  determined  by  Ax. 


When  0  =  0,  n  =  0  regardless  of  the  fluctuations  in  Ax.  This  is  the  backscat- 
tering  condition  —  all  contributions  are  in  phase.  However,  higher  n  values  will 
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not  preclude  the  variance  of  A.r.  The  prominence  of  any  fringe  at  these  values  will 
depend  upon  the  variance  in  the  phase.  The  standard  deviation  in  the  phase  may 
be  determined  from  equation  (6.3), 

4  7T  . 

a<t>  ~  Sln  (<>-7) 

The  larger  o>  is,  the  smaller  the  intensity  of  any  higher  order  fringes.  Obviously,  to 
determine  the  exact  condition  the  general  form  of  the  probability  density  function 
of  Ax  must  be  found.  The  values  of  //Ar  and  <?Ar  may  also  then  be  obtainable. 

Antibackscatterl43!  is  also  allowed  by  this  model,  since  it  has  now  been  established 
that,  time  reversed  paths  are  not  required  for  a  peak  in  an  intensity  envelope. 

One  final  feature  of  the  intensity  plots,  which  has  not  been  previously  addressed  and 
which  this  model  can  explain,  is  why  the  co  polarised  and  cross-polarised  compo¬ 
nents  of  the  detected  light  show  peaks  at  different  angles.  This  is  evident  in  figure 
5.23  in  chapter  five,  and  also  for  the  side  lobes  in  other  scattered  intensity  results 
from  two  dimensional  surfacW1’^"*. 

Figure  6.3  shows  intensity  results  similar  to  those  shown  in  chapter  five,  but  was 
obtained  by  summing  up  all  of  the  pixels  in  a  speckle  image  taken  with  a  Pulnix 
CCD  camera^'8).  This  camera  was  used  before  the  lower  noise  Princeton  Instruments 
camera  was  available.  The  general  shape  of  these  curves  should  be  ignored  —  they 
result  from  an  automatic  gain  which  the  camera  possessed  —  but  the  angles  at. 
which  the  peaks  appear  may  be  relied  upon.  The  images  which  make  up  the  cross 
polarised  result  include  some  of  those  which  went  to  make  up  the  surface  plot  in 
chapter  five,  figure  5.22. 


Figure  6.3  is  considered  now  as  it  is  known  to  results  from  an  area  of  surface  #83 
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from  which  correlation  results  have  been  obtained  and  it  clearly  shows  the  positions 
of  the  peaks.  It  should  be  stressed  that  the  following  will  still  apply  to  figure  5.23 
in  chapter  five. 


Figure  6.3:  Intensity  plots  for  s  polarised  and  p  polarised  light  detected  from  sur¬ 
face  #83,  showing  the  angles  of  the  peaks  observed  for  observation  in  the  specular 
direction  when  s-polarised  light  is  incident  with  wavelength  0.633////?.  Data  points 
were  taken  ever)  half  degree. 


It  is  now  necessary  to  consider  what  is  happening  from  a  two  dimensional  point  of 
view  and  introduce  the  concept  of  a  mean  well  (i.e.,  the  average  dip  in  the  surface). 
Consider  one  such  well,  which,  upon  double  scattering,  may  give  rise  to  both  a 
co- polarised  and  cross-polarised  component,  as  shown  in  figure  6.1.  It  is  clearly 
evident  in  the  diagram,  that  the  mean  well  width  from  the  point  of  view  of  the 
cross-polarised  component  is  shorter  than  that  fer  the  co-polarised  component,  i.e., 
for  s-polari'-ed  light  incident. 


7,P  '  iis 

7lAx  ’/Ai 


(6.S) 


where  the  superscript:-  represent  the  polarisation  state  of  the  observed  scattered 


Figure  6.4:  Schematic  diagram  showing  how  a  well  in  a  sunaoe  doubly  scatters 
polarised  light  to  give  both  a  co-polarised  and  cross-polarised  cop  ponent. 


Considering  equation  (0.6),  it  may  he  seen  that  for  a  part-ion, .-l-  n  (greater  than 
zero), 

0,  >  0,  (0.9) 

i.e.,  as  observed,  the  angle  at  which  a  peak  of  a  particular  order  appears,  \s  larger 
for  the  cross-polarised  component  titan  for  the  co-polarised  one. 

Continuing  with  the  model  of  a  circular  (the  average  shape)  mean  well,  a  more 
definite  relationship  between  //^r  and  r/^r  may  ire  established,  purely  due  to  the 
circular  geometry.  From  figure  6.4, 


?/L  =  '/I*  cos  45 


(6.10) 


'Fhis  lead.*'  to  a  re'atioi;  between  the  angles  at  which  the  peaks  occur;  for  a  particular 


*  ■’!  "  C 


1^1 


.•  » 
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* 
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sin  0p  —  r/dsin  0 , 


(6.1 1 ) 
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Taking  the  angle  at  which  the  s-polarised  side  lobe  appears  in  figure  6.3,  03  = 
6.5  ±  0.25  degrees,  and  applying  equation  (G.ll),  the  angle  the  p-polarised  side  lobe 
should  be  evident  is,  0p  =  9.2  ±  0.4  degrees.  This  compares  well  with  the  angle 
measured  from  figure  6.3  of  9.25  ±  0.25,  giving  a  great  deal  of  weight  to  equation 
(6.11)  and  the  simple  mean  well  model. 

6.2.2  Correlation  Considerations 

Following  an  argument  similar  to  that  given  at  the  beginning  of  the  previous  section, 
a  high  degree  of  correlat  ion  may  be  said  to  exist  when,  on  average,  the  difference  in 
the  phase  differences  between  any  two  pairs  of  interfering  paths  is  close  to  zero.  Note 
that  this  is  regardless  of  what  those  phase  differences  are.  This  is  built,  on  the  basis 
that  a  high  degree  of  correlation  occurs  when  speckles  appear  in  the  same  place, 
•vii.  respect  to  other  speckles,  that  they  do  at  another  angle.  The  presumption  is 
then  made  that  individual  wells  in  the  surface  must  be  making  the  same  (or  similar) 
contributions  to  the  scattered  field  at  both  angles. 

Although  it  may  seem  unlikely  that  the  multiple  contributions,  from  the  various 
wells,  should  all  maintain  the  same  phase  relationships  with  respect  to  each  other, 
without  keeping  the  same  phase,  there  is  no  reason  for  imposing  such  a  restriction 
and  therefore  the  absolute  pha.'e  change  will  be  left  undefined. 

Starting  with  the  phase  difference  between  a  pair  of  paths  on  the  surface',  equation 
(6.3),  the  difference  between  that  intrcduccd  by  one  pair  and  that  introduced  by 
another  may  be  obtained.  In  particular,  the  difference  introduced  by  one  angle  ol 
incidence  compared  to  another  may  be  determined. 
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Figure  6.5:  A  pair  of  the  many  possible  interfering  paths  within  a  valley  of  a  rough 
surface  shown  at  different  angles  of  incidence. 


Comparing  figure  6.1  with  figure  6.5,  it  may  be  seen  that  the  latter  consists  of  two 
of  the  interfering  path  pairs  of  the  former,  incident  at  different  angles.  For  now,  the 
path  the  light  takes  between  scattering  points  on  the  surface  is  assumed  the  same 
for  all  angles  for  the  purposes  of  comparison. 

The  phase  differences  within  each  pair  are 


0  1  = 

4/r 

—  sin  Vj  A.r 

A 

(6.12) 

02  = 

Aw  . 

—  sir.  x 

A 

(6.15) 

where  the  subscripts  1  and  2  denote  the  two  different  pairs. 


The  difference  in  the  phases  between  one  angle  of  incidence  and  the  other  will 
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therefore  be 

4  7T 

A<i>  =  (p2  —  4>\  —  "T- (sin 02  -  sin Oi)Ax  (6.14) 

A 

In  equation  (6.14)  Ax  may  by  replaced,  following  the  same  argument  as  made  pre¬ 
viously,  by  7/^.c  —  a  surface  wide  parameter. 

As  stated  above,  a  high  correlation  will  be  expected  when  the  difference  in  the  two 
phases  is  zero  or  an  integer  multiple  of  2n.  Introducing  this  into  equation  (6.14) 
gives 

u  -\ 

sin  02  =  sm  6:i  +  -  (C.15) 

When  n  =  U,  0\  —  02,  a  correlation  between  two  of  the  same  thing  —  obviously  I 
—  the  main  correlation  peak  in  the  experimental  results  in  chapter  five. 

However,  higher  order  interactions  may  still  give  a  degree  of  correlation,  albeit  at¬ 
tenuated  by  the  fluctuations  in  Ax,  explaining  the  extra  peaks  observed  in  the 
correlation  of  speckle  patterns  bom  doubly  scattering  surfaces. 

Unfortunately,  experimental  work  conducted  did  not  produce  more  than  one  extra 
correlation  peal;  (the  first  order  one),  but  the  computational  results,  shown  in  figure 
5.26  in  chapter  five,  clearly  show'  two  extra  peaks  to  the  side  of  the  main  n  =  0 
peak.  This  evidence  affords  some  confidence  in  the  above  argument. 

G.2.3  Further  Analysis  of  the  Simple  Mode! 

The  model,  as  it  stands,  is  not  perfect;  for  instance,  it  fails  to  predict  the  relative 
placement  of  any  higher  order  intensity  peaks  and  is  likely  to  fail  in  the  correlation 


1S2 
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case  as  well.  Given  that  the  first  order  peak  is  at  0 i,  the  second  order  peak  should 
be  expected,  from  equation  (6.6),  at 

sin  #2  =  2  sin  Oy  (6.16) 

In  figure  6.3,  with  the  first  order  peak  at  6.5  ±  0.25  degrees  tor  the  co-polarised 
situation,  the  second  order  peak  should  be  expected  at  13.1  ±0.5  degrees.  There  is 
clearly  a  dip  around  this  region  in  the  graph. 

It  is,  however,  possible  to  explain  this  failure.  One  of  the  main  assumptions  the 
model  was  based  upon  was  that  the  mean  well  width  is  fixed,  and  that  the  interfering 
light  traveled  along  equal  but  opposite  paths  between  the  two  scattering  points  of 
surface. 

It  is  now  proposed  that  if  light  is  incident  on  a  well  at  a  large  angle,  it  will  travel 
a  different  distance  between  points  than  that  im  ident  at  a  small  angle.  Tins  can 
be  seen  in  figure  6.6.  Tracing  the  concept  through,  it  is  realised  that  if  a  larger 
angle  of  incidence  decreases  the  width  between  scattering  points,  then  the  angle  the 
peak  appears  at  will  correspondingly  increase.  Releasing  the  restriction  of  fixing  y±r 
explains  why  a  second  order  peak  is  observed  at  around  15.5  degrees,  in  figure  6.3, 
instead  of  the  predicted  position  of  around  13  degrees. 

The  possible  scattering  paths  that  might  exist  higher  up  the  well,  traversing  across 
larger  distances,  may  be  discounted  by  realising  that,  in  a  large  proportion  of  these 
cases,  the  other  side  of  the  well  is  not  likely  to  rise  high  enough  to  catch  the  light 
for  the  second  scatter. 

Note  that  now  should  be  written  either,  7/&.-(n),  as  a  function  of  n,  or  more 
properly,  r/Ar(0),  as  a  function  of  the  angle  of  incidence.  Sonic  idea  of  the  form  of 
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Figure  6.6:  Diagram  showing  how  an  increase  in  the  angle  oi'  incidence  leads  to  a 
decrease  in  the  distance  traversed  across  the  well. 

this  function  may  be  gained  bv  considering  the  gradients  of  the  sides  of  a  mean  well. 

It  is  first  necessary  to  make  an  assumption  that  in  double  scattering,  for  scattering 
in  the  specular  direction,  the  paths  between  scattering  points  are  horizontal  (by 
which  it  is  meant  that  they  are  parallel  to  the  mean  surface  plane).  This  may  be 
justified  by  realising  that  if  a  mean  well  is  being  considered,  this  will  be  symmetrical 
and  the  light  will  therefore  travel  along  a  horizontal  path.  As  discussed  in  chapter 
four,  and  shown  in  figure  4.6,  this  will  lead  to  a  gradient  at  the  surface  of 


where  0  is  constrained  to  be  positive,  and  m  is  negative  if  0  is  the  angle  of  incidence, 
and  positive  if  0  is  the  angle  of  observation  (for  incident  light  from  the  left  in  the 
above  diagrams;.  The  change  in  the  gradient  between  the  two  scattering  points  is 
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therefore  2m,  the  difference  between  —in  and  in. 

In  chapter  two,  an  expression  was  obtained  for  the  rate  of  change  of  the  gradient 
with  respect  to  the  distance  across  a  Gaussian  surface,  namely  its  mis  curvature, 
equation  (2.33). 


As  the  portion  of  the  surface  considered  here  is  at  the  bottom  of  a  well,  i.e.,  in  the 
region  of  a  high  point  in  the  curvature,  the  approximation  to  the  mis  curvature  may 
be  on  the  low  side.  To  allow  for  this  a  quantity  p  (which  gives  a  result  of  less  than 
1)  is  introduced, 


dm 

dr 


2y/'3cr 
PT 2 


(6. 18) 


<t  and  r  are,  respectively,  the  rms  height  and  correlation  length  of  the  surface. 
Strictly  speaking,  p  will  be  a  function  of  0 ,  decreasing  as  0  increases  —  for  now  it 
will  be  considered  constant. 


Dividing  the  change  in  the  gradient  by  this  rate  of  change  will  give  an  expression 
for,  7 the  distance  between  the  two  scattering  points. 


= 


pt 


r  0 
tan  |  —  —  - 


(6. lb) 


where  0  is  the  angle  of  incidence.  Obviously  this  expression  will  apply  to  the  co¬ 
polarised  case  —  it  is  related  to  the  cross-polarised  case  by  equation  (6.10). 


Substituting  equation  (6.10)  into  equation  (6.6),  an  absolute  expression  for  the  an¬ 
gle:  ,  0 ,  at  which  the  side  lobes  appear  may  be  obtained, 


sin  0  tan  1 - -  1  = 

•>  9 


2pr2 


(6.20) 


which  may  be  transformed  with  suitable  trigonometric  identities  to  give  (for  terms 


6.2.  DOUBLE  SCATTER 


185 


up  to  order  tan2(d/2)), 


(o\  (2  -  0  -  v'i'2  -  O' 

tan  I  —  ss - 1 - 

V2 )  2(C  +  2) 


-4C(C  +  2) 


(6.21) 


where 


n\  \/3o 
‘Ipr1 


(6.22) 


Applying  this  to  surface  #S3,  p  may  be  determined  by  looking  at  the  n  =  1  peak  in 
figure  6.3  where  0  =  6.5  ±  0.25.  It  is  found  that  p  =  0.66  ±  0.1  -1  (the  large  error  is 
due  to  the  inclusion  of  the  uncertainty  in  the  measured  surface;  parameters,  <r  and 

r). 


For  n  =  2,  equation  (6.20)  predicts  0  —  15  ±  5  degrees  (which  will  be  slightly  on 
the  low  side  due  to  the  assumption  that,  p  is  fixed).  It  may  be  seen  that  the  second 
peak  in  figure  6.3  is  near  the  centre  of  this  range  at  15.5  ±  0.25  degrees, 

Incidentally,  the  mean  well  widths  which  give  rise  to  these  intensity  peaks  are  ap¬ 
proximately  4.3 pin  for  //  =  1  and  3.6//)/!  tor  n  —  2.  That  is  just  under  7  and  6 
wavelengths  respectively. 

The  calculated  value  of  p  for  figure  6.2  is  0.94  ±  0.20.  The  difference  between  this 
value  and  that  found  from  figure  6.3  can  be  attributed  to  the  different  areas  of  the 
surface  used  in  each  case,  which  implies  that  the  surface  statistics  are  not  uniform 
over  the  whole  surface  on  such  a  small  scale  (and  was  found  to  be  the  case^'J). 

The  relaxation  of  the  constraint  on  ij&x  also  carries  through  to  the  correlation  model. 
The  assumption  made  at  the  start  of  this  study,  that  the  two  scattering  paths 
traverse  the  same  line  between  identical  scattering  points,  (refer  back  to  figure  6.5 
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to  see  this),  is  no  longer  consistent  with  the  understanding  of  the  processes  at  this 
stage. 

As  previously  mentioned,  for  a  high  degree  of  correlation,  all  that  is  required  is  for 
any  two  paths  to  give  rise  to  the  same  phase  change.  It  is  not  necessary  for  these 
paths  to  travel  the  same  route. 

Rewriting  equation  (G.l-1)  to  reflect  this,  one  obtains 

Nq  =  02  —  Oi  =  —  (sinfl2A.r-2  —  sin  01  A.iq )  (6.2:1) 

where  Ax  is  not  fixed,  but.  depends  upon  the  angle  of  incidence.  The  subscripts  1 
and  '2  denote  the  two  different  correlated  arrangements. 

Shifting  again  into  using  the  mean  width  (which  is  now  a  function  of  angle)  and 
confining  the  two  phases  to  be  integer  numbers  of  2jt  apart, 

2 

Cn\02i),u  (C2)  -  ~  +  sinfhy^dth)  (6  21) 

J>  A 


Substituting  equation  (6.19)  (modified  for  the  cross-polarised  case)  into  this  gives 


x  Oi 


y/bnXa 


sin  02  tan  -  -  --  =  -  +  sin  0l  tan  I  -  -  - 


2/7  r  2 


7T  9X 


(6.2o) 


which  may  be  solved  for  02  using  equation  (6.21),  where 

(H 


\Z(miXg 
(,  =  — — — - h  sin  0 1  tan 


2  pr2 


(6.26) 


Unfortunately,  as  discussed  above,  there  is  no  reason  to  expect  the  phase  differences 
in  the  correlation  peaks  to  be  related  to  the  intensity  envelope  peaks.  In  fact,  the 
paths  which  combine  to  give  a  high  degree  of  correlation  (i.e,,  speckles  in  similar 
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positions)  are  unlikely  to  be  the  same  paths  which  produce  the  high  points  in  t.he 
intensity  envelopes.  However,  the  differences  introduced  into  the  paths  by  illumi¬ 
nating  at  different  angles  are  adequately  covered  by  equation  (6.2b)  if  p  is  set  to 
allow  for  the  curvatures  involved. 


Examining  figure  5.20,  it  may  be  seen  that  for  the  reference,  0 at  2  degrees,  a  peak 
in  the  correlation  curve  appears  when  Oj  —  S.5  ±  0.25  degrees.  This  gives  a  value 
for  p  in  the  correlation  arrangement  of  1.0  ±  0.3  (including  all  errors  due  to  the 
uncertainty  in  measurement  of  the  surface  parameters).  It  has  been  assumed  that 
the  observed  peak  is  the  first,  order  interference  one,  i.e.,  for  n  —  1. 


With  the  reference  at  0 1  =  10  degrees,  there  is  an  observed  peak  at  n  —  —1.  Its 
predicted  position,  at  0}  —  3  ±  1.5  degrees,  compares  favourably  with  t.lte  measured 
angle  of  2.5  ±  0.25  degrees.  However,  note  that  the  angular  differences  between  the 
two  peaks  and  their  reference  angles  are  too  similar  to  pose  a  strenuous  test  for  the 
theory.  The  fluctuations  in  A.r,  coupled  with  the  spatial  averaging  necessarily  in¬ 
volved  with  the  experimental  technique  used,  prevent  any  higher  ordei  peaks,  which 


could  rnoie  rigorously  test  the  predictions,  from  being  observed  in  the  correlation 
experiments. 


It  may  be  realised  that  the  larger  value  of  p  for  this  situation,  means  that  the 
curvature,  averaged  around  the  area  of  a  well  which  contributes  to  a  high  degree 
of  correlation,  is  less  than  that  in  the  intensity  ease;  i.e.,  the  scattering  processes 
which  contribute  to  the  t  wa  correlation  in  surface  #83  occur  at  higher  positions 
in  the  well  than  do  those  which  contribute  to  a  side  lobe  in  the  intensity  envelope. 


As  mentioned  above,  p  was  assumed  to  be  constant.  The  assumption  will  have  a 
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much  larger  effect  on  the  correlation  results,  than  it  did  in  the  intensity  situation, 
as  the  curvature  is  "likely  to  be  varying  more  rapidly  in  the  region  of  the  surface 
contributing  in  the  former  case.  This  statement  is  based  on  the  previous  finding, 
that  the  portion  of  the  surface  contributing  to  the  peaks,  is  further  away  from  the 
bottom  of  a  well  for  the  correlation  case  as  compared  to  the  intensity  case.  As  this 
is  not  allowed  for,  the  predictions  of  the  positions  of  any  correlation  peaks  should 
not  be  expected  to  match  exactly  with  the  observed  positions. 

Further  evidence  of  the  change  in  /»  with  tin*  angle  of  incidence  is  seen  in  the  corre¬ 
lation  coefficients  lor  scattering  from  surface  #23t).  Considering  figure  1.20,  values 
for  />  are  found  to  be  about  0.6  around  the  10  degree  point,  increasing  to  about  2.0 
by  the  30  degree  mark.  These  movements  may  be  explained  if  at  higher  angles  of 
incidence,  the  areas  of  a  particular  well  contributing  to  extra  correlation  peaks  are 
higher  up  the  sides  of  that  well;  i.e.,  the  average  curvature  between  the  two  points 
is  reduced.  The  movement  of  the  contributing  paths  to  higher  positions  in  the  well 
also  explains  why  the  extra  correlation  peaks  are  lower  at  the  higher  angles  —  the 
fluctuations  in  A.v  will  be  greater  at  higher  angles. 

6.3  Concluding  Remarks 

The  first  three  chapters  of  this  thesis  were  concerned  with  building  up  the  back¬ 
ground  and  the  concept  of  multiple  scattering  and  correlation  coefficients.  An  intro¬ 
duction  was  given  into  the  statistical  notions  and  electromagnetic  theories  used  later 
in  the  work;  and  the  methods  and  controls  required  to  produce  results  presented  in 
the  following  chapters  were  detailed. 
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Perhaps  the  only  real  surprise  contained  within  these  three  chapters  was  at  the  end 
of  chapter  two,  where  a  statistical  study  gave  an  analytical  tool  for  determining 
whether  a  Gaussian  random  rough  surface  falls  within  the  regime  of  the  Kirch- 
hoff  approximation.  A  distinction  was  made  between  the  validity  of  the  KirchhofF 
approximation,  as  applicable  to  a  point  on  the  surface,  and  the  single  scattering 
Beckmann  theory.  The  three  surfaces  used  throughout  the  work  were  all  seen  to  fall 
within  the  bounds  of  the  Kirchhoff  approximation  (at  least  at  normal  incidence). 

The  fourth  chapter  saw  the  start  of  more  original  work;  it  examined  the  angular 
correlation  of  speckle  patterns  produced  by  low  slope  surfaces.  A  study  similar  to 
that  used  to  determine  which  surfaces  fall  within  the  Kirchhoff  approximation  was 
employed,  to  give  some  idea  of  whether  a  surface  should  be  expected,  based  on  its 
measured  parameters  a  and  r,  purely  to  be  a  single  scattering  one  or  not.  It  is  shown 
that  of  the  three  surfaces  under  examination,  one  was  purely  single  scattering,  one 
double  scattering,  and  the  third  lay  on  the  border  between  the  two  regions. 

An  analytical  expression  was  derived,  following  the  work  of  Beckmann'1'*  and  Leger 
and  Perrin-31',  for  the  intensity  correlation  coefficients  expected  in  a  general  two 
dimensional  scenario.  It  was  then  tailored  to  suit  in  plane  scattering  and  some  of  its 
possibilities  were  explored,  these  departed  from  the  results  one  might  have  expected 
intuitively.  The  more  general  expression  was  averaged  around  the  direction  of  ob¬ 
servation,  to  represent  more  accurately  the  relatively  wide  view  considered  in  the 
experimental  case.  Experimental  results  from  a  single  scattering  surface  were  found 
to  be  in  good  agreement,  and  reexamining  the  possibilities  explored  earlier  with 
this  new  averaged  technique  gave  results  in  line  with  the  intuitive  understanding. 
Finally,  results  were  presented  for  a  surface  which  departs  slightly  from  the  single 
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scattering  regime. 

In  chapter  five,  the  double  scatter  Kirchhoff  theory  was  derived  for  two  dimen¬ 
sional  scattering.  It  had  to  be  solved  computationally  and  most  of  this  chapter  was 
concerned  with  explaining  the  methods  used  to  make  it  possible.  Polarisation  was 
incorporated  into  the  calculations,  and  surface  self  shadowing  was  allowed  for.  A 
technique  capable  of  calculating  which  points  were  visible  from  another,  over  tiu* 
whole  surface  in  just  one  pass,  was  presented.  Experimental  results  and  tin'  results 
from  the  computational  model  were  shown. 

In  the  previous  sections  of  this  final  chapter,  an  accepted  model  used  to  explain 
intensity  and  correlation  results,  which  succeeded  in  predicting  the  occurrence  of 
certain  features,  was  shown  to  be  insufficient  at  a  less  superficial  level.  The  model 
was  advanced  and  shown  to  adequately  explain  and  predict  various  features  evident 
in  intensity  and  correlation  situations,  when  observing  in  the  specular  direction, 
more  accurately. 

It  was  shown  that  the  positioning  of  the  wells  within  the  surface  does  not  matter, 
only  their  widths  have  any  effect,  upon  the  angular  positioning  of  multiple  scatteiing 
features.  Also,  it  was  suggested  that  it  was  a  region  near  the  bottom  of  a  well  which 
gave  rise  to  the  intensity  envelope  side  lobes,  and  that  less  restriction  is  placed  upon 
the  region  giving  rise  to  the  extra  correlations. 

This  may  be  explained  by  considering  that  for  an  overall  intensity  gain,  i.e.,  a 
side  lobe,  a  vast  majority  of  the  wells  on  a  surface  must  be  introducing  similar 
phase  changes  (so  that  their  contributions  to  the  scattered  field  will  constructively 
combine).  It  seems  reasonal.de  to  expect  that  any  similarities  in  the  various  wells  on 
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the  surface,  will  Ire  towards  their  bottoms  -  every  well  has  a  bottom,  but  how  far 
up  its  sides  go  is  less  delinite. 

For  a  high  degree  of  correlation,  it  is  only  required  that  a  single  well  should  give  the 
same  phase  contribution  to  the  scattered  field  at  one  angle  as  it  did  at  the  other 
angle  being  considered.  It  should  be  expected  that  this  will  occur  surface  wide  (on 
an  average  basis)  at  particular  angles  only,  but  not  necessarily  at  ihose  angles  which 
give  rise  to  intensity  side  lobes. 

The  main  successes  of  this  model  were  the  relationships  formed  between  peaks  of 
different  orders  in  both  the  intensity  envelope  and  correlation  curves,  and  between 
the  co-  and  cross-polarised  peaks  (of  a  particular  order)  observed  in  the  intensity 
results. 

No  attempt  was  made  to  reformulate  the  equations  derived,  such  that  they  may  be 
applied  in  a  more  general  situation  than  in  the  specular  direction.  However,  the 
basic  principles  are  not  affected  by  the  geometry;  it  should  he  a  relatively  simple 
matter  to  extend  the  model  beyond  the  presented  specular  scenario. 

As  a  final  point,  it  would  be  desirable  to  have  the  model  rebuilt  on  a  more  ngorous 
basis,  using,  say,  the  moments  of  A.r  or  even  finding  a  function  for  it.  It  should  then 
be  possible  to  exactly  predict  the  angular  position  of  every  feature  in  tin*  correlation 
and  intensity  plots,  and  maybe  even  gain  some  idea  of  their  relative  contrasts. 
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